Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 900,246 C→T S66L (TCA→TTA)  SAUSA300_RS04420 → cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793900,2460CT100.0% 283.1 / NA 95S66L (TCA→TTA) SAUSA300_RS04420cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (90/5);  total (90/5)

CGACAAGTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATCATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACACGCGGAACAACTGC  >  NC_007793/900150‑900345
                                                                                                |                                                                                                   
cGACAAGTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGattatt                                                                                                  >  1:4954580/1‑100 (MQ=255)
 gACAAGTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTg                                                                                                 >  1:2844058/1‑100 (MQ=255)
    aaGTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCa                                                                                              <  1:2853339/100‑1 (MQ=255)
    aaGTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCa                                                                                              <  1:2419160/100‑1 (MQ=255)
      gTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAAc                                                                                            <  1:3114862/100‑1 (MQ=255)
      gTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAAc                                                                                            >  1:2799152/1‑100 (MQ=255)
        cAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTg                                                                                          >  1:3989875/1‑100 (MQ=255)
        cAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTg                                                                                          >  1:5012832/1‑100 (MQ=255)
        cAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGACAACTg                                                                                          >  1:4306342/1‑100 (MQ=255)
         taaCGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGTGTTATAATTCAAACGTTCATCGTGGTGTCCATACATTAGGATTATTGGCAACTGa                                                                                         >  1:565266/2‑100 (MQ=255)
           aggCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATg                                                                                       >  1:2702601/3‑100 (MQ=255)
           aCGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATg                                                                                       >  1:3362232/1‑100 (MQ=255)
           aCGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATg                                                                                       >  1:2620494/1‑100 (MQ=255)
            tgCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATgg                                                                                      >  1:1709964/2‑100 (MQ=255)
              ccTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACATGCGTTATAATTCAAACGTTCATCGTGGTGTTAATACATTAGGATTATTGGCAACTGATGGtt                                                                                    >  1:4351774/1‑100 (MQ=255)
              ccTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGtt                                                                                    >  1:648877/1‑100 (MQ=255)
              ccTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGtt                                                                                    >  1:20639/1‑100 (MQ=255)
              ccTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGtt                                                                                    >  1:5704740/1‑100 (MQ=255)
               cTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTa                                                                                   >  1:1832055/1‑100 (MQ=255)
               cTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTa                                                                                   >  1:772581/1‑100 (MQ=255)
               cTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACAGAtcggaa                                                                                   >  1:1817578/1‑95 (MQ=255)
                tgtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTAt                                                                                  >  1:4614789/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACtagatcggaag                                                                                 >  1:4795161/1‑90 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGacgagatcg                                                                                 >  1:1473598/1‑92 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:1921288/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:5539784/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:2849508/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:2840364/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:4719025/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:4615267/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:335287/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:3828668/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:2660116/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:4544586/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:430976/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:3519734/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:3168124/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATg                                                                                 >  1:3685521/1‑100 (MQ=255)
                 gtgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGAGTTCATACATTAGAATTATTGGCAACTGATGGTTATg                                                                                 >  1:2660564/1‑100 (MQ=255)
                 gggCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAAGATTATTGGCAACTGATGGTTATg                                                                                 >  1:2663655/3‑100 (MQ=255)
                  cgCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGGGTTCATACATTAGGATTATTGGCAACTGATGGTTATGa                                                                                >  1:396811/2‑100 (MQ=255)
                   gCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGGTTATTGGCAACTGATGGTTATGaa                                                                               >  1:3589860/1‑100 (MQ=255)
                   gCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGaa                                                                               >  1:5735031/1‑100 (MQ=255)
                   gCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGaa                                                                               >  1:937195/1‑100 (MQ=255)
                   gCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGaa                                                                               >  1:969797/1‑100 (MQ=255)
                     aaGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGaaaa                                                                             >  1:5750836/1‑100 (MQ=255)
                         gtTAAATGTTTTAGAAGATTACTACAAGCGTCATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGcc                                                                         >  1:4473023/1‑100 (MQ=255)
                            aaaTGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGCGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGt                                                                      >  1:610713/1‑100 (MQ=255)
                                     agaagaTTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGtt                                                             >  1:131775/1‑100 (MQ=255)
                                           ttACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTcgtcgt                                                       >  1:791576/1‑100 (MQ=255)
                                             actacAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGttt                                                     >  1:907926/1‑100 (MQ=255)
                                             actacAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGttg                                                     <  1:2318194/100‑2 (MQ=255)
                                               tacAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTtta                                                   >  1:1546974/1‑100 (MQ=255)
                                               tacAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTtta                                                   >  1:388762/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:5468034/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:420247/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:4607127/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:394879/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:1288798/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:1367927/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:3055017/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGCTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:1544770/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGTGCTCATACATTAGGATTATTTGCAACTGATGGTTATGAAAATGCCCGTGAAACTGGTCGTCGTTttat                                                  >  1:219543/1‑100 (MQ=255)
                                                acAAGCGTTATAATTCAAACGTTCATCGTGGGGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTttat                                                  >  1:1431866/1‑100 (MQ=255)
                                                        tataATTCAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGa                                          <  1:5339422/100‑1 (MQ=255)
                                                               cAAACGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAACGCGAAgatcgg                                   >  1:4239822/1‑95 (MQ=255)
                                                                  aggTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTga                                >  1:4070478/3‑100 (MQ=255)
                                                                   cGTTCATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAgatcggaaga                               >  1:3525185/1‑91 (MQ=255)
                                                                       cATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGaaa                           >  1:152953/1‑100 (MQ=255)
                                                                       cATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGaaa                           >  1:1393647/1‑100 (MQ=255)
                                                                       cATCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGaaa                           >  1:3074553/1‑100 (MQ=255)
                                                                         tCGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATc                         >  1:2132079/1‑100 (MQ=255)
                                                                          ggtggtgTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCa                        >  1:5586975/2‑100 (MQ=255)
                                                                          cGTGGTGTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCa                        >  1:2193742/1‑100 (MQ=255)
                                                                           gtggtgTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCAt                       >  1:1472357/1‑100 (MQ=255)
                                                                           gcggtgTTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCAt                       >  1:2671546/3‑100 (MQ=255)
                                                                             ggtgTTCATACATTAGGATTATTGGTAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGATGTACTTTGATGAAATCAttt                     >  1:51892/1‑100 (MQ=255)
                                                                              tgcttCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCAtttt                    >  1:588838/4‑100 (MQ=255)
                                                                                gtTCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTa                  >  1:1021164/1‑100 (MQ=255)
                                                                                 gtCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTac                 >  1:3055378/2‑100 (MQ=255)
                                                                                  tCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:5156183/1‑100 (MQ=255)
                                                                                  tCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:5533420/1‑100 (MQ=255)
                                                                                  tCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:3982386/1‑100 (MQ=255)
                                                                                  tCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:2200727/1‑100 (MQ=255)
                                                                                  tCATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:2266061/1‑100 (MQ=255)
                                                                                  ccATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:2957693/2‑100 (MQ=255)
                                                                                  ccATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGCCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTaca                >  1:5110112/2‑100 (MQ=255)
                                                                                   gaTACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTacac               >  1:1043860/2‑100 (MQ=255)
                                                                                   cATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTacac               >  1:5553080/1‑100 (MQ=255)
                                                                                   cATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTacac               >  1:3373303/1‑100 (MQ=255)
                                                                                   cATACATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTacac               >  1:1307689/1‑100 (MQ=255)
                                                                                     cccATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACAcgc             >  1:1916675/3‑100 (MQ=255)
                                                                                      aCATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACAcgcg            >  1:4289985/1‑100 (MQ=255)
                                                                                      aCATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACAcgcg            >  1:240614/1‑100 (MQ=255)
                                                                                       cATTAGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACACGCgg           >  1:3799961/1‑100 (MQ=255)
                                                                                           aGGATTATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACACGCGGagat       >  1:1774555/1‑97 (MQ=255)
                                                                                                cATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACACGCGGAACAACTGc  >  1:3177365/1‑100 (MQ=255)
                                                                                                |                                                                                                   
CGACAAGTCAAACGCCTGTGCAAGTGTTAAATGTTTTAGAAGATTACTACAAGCGTTATAATTCAAACGTTCATCGTGGTGTTCATACATTAGGATCATTGGCAACTGATGGTTATGAAAATGCCCGTGAAACCGTTCGTCGTTTTATTAATGCGAAGTATTTTGAAGAAATCATTTTTACACGCGGAACAACTGC  >  NC_007793/900150‑900345

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑