Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1624 250.6 12620048 98.0% 12367647 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,170 +C intergenic (+172/+279) hemN → / ← glnG O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,1701.C100.0% 647.9 / NA 231intergenic (+172/+279)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (140/91);  total (140/91)

Only 100 of 251 total aligned reads displayed.

TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACA  >  NC_000913/4051077‑4051256
                                                                                              |                                                                                      
tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                                                         >  2:3707697/1‑93 (MQ=255)
tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcg                                                                                         <  1:3065786/93‑1 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                                        >  2:5554043/1‑92 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                                        >  2:4879700/1‑92 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                                        <  2:364044/93‑2 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                                        <  1:4987257/93‑2 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATgcggcgc                                                                                        <  2:6247635/93‑2 (MQ=255)
      acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                                                   >  1:5908137/1‑93 (MQ=25)
      acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                                                   >  1:4588911/1‑93 (MQ=25)
        aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                 >  2:5376000/1‑93 (MQ=25)
        aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                 >  2:202863/1‑93 (MQ=25)
         aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                                                >  2:506437/1‑93 (MQ=21)
         aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                                                >  2:1622416/1‑93 (MQ=25)
         aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTAAGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                                                >  1:4090216/1‑93 (MQ=25)
           aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                                              >  1:5299381/1‑93 (MQ=21)
           aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                                              >  2:1958348/1‑93 (MQ=21)
            ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                                             >  1:139338/1‑93 (MQ=21)
            ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                                             <  2:2499436/93‑1 (MQ=25)
             cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                            >  1:5947526/1‑93 (MQ=21)
                 tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                        <  2:5429924/93‑1 (MQ=21)
                 tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                        <  2:4961131/93‑1 (MQ=21)
                  gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                                       >  1:4562755/1‑93 (MQ=18)
                   aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                                                                      >  1:2786110/1‑93 (MQ=18)
                    aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                     <  1:1598191/93‑1 (MQ=18)
                      ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                   <  2:1811236/93‑1 (MQ=18)
                       tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                  >  2:3801757/1‑93 (MQ=18)
                       tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                  >  1:1166656/1‑93 (MQ=18)
                        attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                                 >  2:1133889/1‑93 (MQ=17)
                        attGCTTTTGTAGGCCGGATAAGGCGCTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                                                 >  1:333243/1‑93 (MQ=17)
                         ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTGCg                                                                >  1:797131/1‑93 (MQ=17)
                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                               >  1:4000136/1‑93 (MQ=18)
                           gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                                              <  2:794349/93‑1 (MQ=17)
                            cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                                             <  1:3039328/93‑1 (MQ=17)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            >  1:1686347/1‑93 (MQ=17)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            >  2:273918/1‑93 (MQ=18)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                            >  2:2596419/1‑93 (MQ=17)
                              tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                                           <  2:5997534/93‑1 (MQ=18)
                               ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                                          <  1:1457813/93‑1 (MQ=18)
                                 gTAGGCCGGATAAGGCATTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTTATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                        <  2:292950/93‑1 (MQ=18)
                                  tAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                       <  1:4199203/93‑1 (MQ=18)
                                   aGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                                      >  2:5590266/1‑93 (MQ=18)
                                    ggCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                                                     >  2:4170849/1‑93 (MQ=18)
                                      ccGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                                   <  2:2272238/93‑1 (MQ=18)
                                      ccGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                                   <  2:551816/93‑1 (MQ=18)
                                       cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                                  <  1:5580421/93‑1 (MQ=18)
                                       cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                                  >  1:1985624/1‑93 (MQ=18)
                                       cGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                                                  <  2:818619/93‑1 (MQ=18)
                                        ggATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                                 >  1:3523968/1‑93 (MQ=18)
                                        ggATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                                 >  1:2699067/1‑93 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                <  2:934254/93‑1 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                <  1:4673747/93‑1 (MQ=18)
                                          aTAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                                               >  2:3322386/1‑93 (MQ=21)
                                           tAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                              >  1:4511702/1‑93 (MQ=21)
                                             aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                                            >  2:4367844/1‑93 (MQ=21)
                                             aGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                                            >  2:4270329/1‑93 (MQ=21)
                                              ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                                           <  2:5867019/93‑1 (MQ=21)
                                              ggCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                                           >  2:4509932/1‑93 (MQ=21)
                                               gCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                                          >  1:1039951/1‑93 (MQ=21)
                                                cGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa                                         <  2:5940901/93‑1 (MQ=21)
                                                cGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa                                         <  1:5145540/93‑1 (MQ=21)
                                                 gTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                        >  1:4367858/1‑93 (MQ=25)
                                                  tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                       >  1:5098779/1‑93 (MQ=21)
                                                  tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                       >  1:5393244/1‑93 (MQ=21)
                                                  tttACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                       <  2:137866/93‑1 (MQ=21)
                                                        ccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                 >  2:2893259/1‑93 (MQ=25)
                                                        ccgcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                 <  2:2996652/93‑1 (MQ=25)
                                                          gcATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                               >  2:2775761/1‑93 (MQ=37)
                                                           cATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                              >  1:4834191/1‑93 (MQ=37)
                                                            aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                             <  1:2602011/93‑1 (MQ=37)
                                                            aTCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                             <  1:6079792/93‑1 (MQ=37)
                                                             tCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                            >  2:6134764/1‑93 (MQ=37)
                                                                aggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                          >  2:3100512/4‑93 (MQ=255)
                                                                 ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                         <  2:377373/91‑1 (MQ=255)
                                                                  gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaa                        >  1:1729237/2‑93 (MQ=255)
                                                                  gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaa                        >  2:3245908/2‑93 (MQ=255)
                                                                  gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaa                        >  1:2877995/2‑93 (MQ=255)
                                                                   cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa                       <  1:64779/93‑1 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                     >  1:4410995/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                     >  1:3295986/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                     <  1:4706399/93‑1 (MQ=255)
                                                                      cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAg                    >  1:1899970/1‑93 (MQ=255)
                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                   >  1:2518427/1‑93 (MQ=255)
                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                   >  1:3972183/1‑93 (MQ=255)
                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                   >  1:3395061/1‑93 (MQ=255)
                                                                         aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                 >  1:3797038/1‑93 (MQ=255)
                                                                         aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                 <  1:2719035/93‑1 (MQ=255)
                                                                          gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCa                >  1:220810/1‑93 (MQ=255)
                                                                           cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa               <  2:2089220/93‑1 (MQ=255)
                                                                           cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa               >  1:3733589/1‑93 (MQ=255)
                                                                            atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt              >  1:785817/1‑93 (MQ=255)
                                                                            atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt              >  1:5598555/1‑93 (MQ=255)
                                                                             tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc             >  2:4795074/1‑93 (MQ=255)
                                                                              gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCa            >  2:4496854/1‑93 (MQ=255)
                                                                              gatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCa            <  2:656058/93‑1 (MQ=255)
                                                                               atgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAt           <  2:4704944/93‑1 (MQ=255)
                                                                                tgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt          <  2:3565460/93‑1 (MQ=255)
                                                                                  ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc        >  2:4093523/1‑93 (MQ=255)
                                                                                    tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCgg      >  1:4861140/1‑93 (MQ=255)
                                                                                      aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGca    <  2:3734911/93‑1 (MQ=255)
                                                                                       tGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac   <  2:1037917/93‑1 (MQ=255)
                                                                                              |                                                                                      
TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACA  >  NC_000913/4051077‑4051256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  >  NC_000913/4051076‑4051263
                                                                                              |                                                                                              
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCG                                                                                                >  SRR4302533.5650432/1‑94 (MQ=60)
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCG                                                                                                >  SRR4302533.2789551/1‑94 (MQ=60)
 TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCG‑C                                                                                              <  SRR4302533.3084616/94‑1 (MQ=60)
  GGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGct                                                                                              >  SRR4302533.3978150/1‑92 (MQ=60)
  GGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGct                                                                                              <  SRR4302533.5017495/94‑3 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          >  SRR4302533.747810/1‑94 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          >  SRR4302533.5943825/1‑94 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          >  SRR4302533.4616708/1‑94 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          <  SRR4302533.3015288/94‑1 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCAGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          >  SRR4302533.1313114/1‑94 (MQ=60)
        CAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGC                                                                                        <  SRR4302533.5132857/94‑1 (MQ=60)
         AAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCG                                                                                       <  SRR4302533.1995143/94‑1 (MQ=60)
         AAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCG                                                                                       <  SRR4302533.3631397/94‑1 (MQ=60)
         AAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCGGGTCGCG                                                                                       >  SRR4302533.1174583/1‑94 (MQ=60)
         AAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTAAGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCG                                                                                       >  SRR4302533.4115117/1‑94 (MQ=60)
          AAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGT                                                                                      >  SRR4302533.2118051/1‑94 (MQ=60)
           AACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTC                                                                                     >  SRR4302533.5331128/1‑94 (MQ=60)
            ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                                                    >  SRR4302533.140225/1‑94 (MQ=60)
            ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                                                    <  SRR4302533.6116099/94‑1 (MQ=60)
             CCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTT                                                                                   >  SRR4302533.5983499/1‑94 (MQ=60)
              CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTA                                                                                  >  SRR4302533.3168882/1‑94 (MQ=60)
              CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTA                                                                                  >  SRR4302533.3872964/1‑94 (MQ=60)
               TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAT                                                                                 <  SRR4302533.4731566/94‑1 (MQ=60)
               TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAT                                                                                 >  SRR4302533.4415618/1‑94 (MQ=60)
               TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAT                                                                                 >  SRR4302533.4448201/1‑94 (MQ=60)
                GTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATC                                                                                >  SRR4302533.5961439/1‑94 (MQ=60)
                  TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                                              >  SRR4302533.2460835/1‑94 (MQ=60)
                  TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                                              >  SRR4302533.3028026/1‑94 (MQ=60)
                  TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                                              >  SRR4302533.4590424/1‑94 (MQ=60)
                   GAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGG                                                                             <  SRR4302533.2701261/94‑1 (MQ=60)
                   GAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGG                                                                             >  SRR4302533.2803227/1‑94 (MQ=60)
                     AATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCC                                                                           <  SRR4302533.1607866/94‑1 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          >  SRR4302533.1564115/1‑94 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          >  SRR4302533.1707442/1‑94 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          <  SRR4302533.748532/94‑1 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          <  SRR4302533.1408315/94‑1 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          <  SRR4302533.2543461/94‑1 (MQ=60)
                       TTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTA                                                                         >  SRR4302533.5250475/1‑94 (MQ=60)
                       TTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTA                                                                         >  SRR4302533.1173655/1‑94 (MQ=60)
                       TTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTA                                                                         <  SRR4302533.4851986/94‑1 (MQ=60)
                        TATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAC                                                                        <  SRR4302533.5698357/94‑1 (MQ=60)
                        TATTGCTTTTGTAGGCCGGATAAGGCGCTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAC                                                                        >  SRR4302533.335282/1‑94 (MQ=60)
                         ATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTGCG                                                                       >  SRR4302533.801823/1‑94 (MQ=60)
                          TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                                                      >  SRR4302533.4024493/1‑94 (MQ=60)
                            GCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACT                                                                    >  SRR4302533.4050198/1‑94 (MQ=60)
                             CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                                                   <  SRR4302533.3057929/94‑1 (MQ=60)
                             CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                                                   >  SRR4302533.5441150/1‑94 (MQ=60)
                             CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                                                   >  SRR4302533.1696561/1‑94 (MQ=60)
                              TTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                                                  <  SRR4302533.959561/94‑1 (MQ=60)
                                TTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCT                                                                <  SRR4302533.1466604/94‑1 (MQ=60)
                                 TGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTC                                                               >  SRR4302533.284571/1‑94 (MQ=60)
                                 TGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTC                                                               >  SRR4302533.4886882/1‑94 (MQ=60)
                                 TGTAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTC                                                               >  SRR4302533.4886881/1‑94 (MQ=60)
                                   TAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                                             <  SRR4302533.3773952/94‑1 (MQ=60)
                                   TAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                                             >  SRR4302533.4198055/1‑94 (MQ=60)
                                   TAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                                             <  SRR4302533.4224615/94‑1 (MQ=60)
                                   TAGGCCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                                             >  SRR4302533.1660694/1‑94 (MQ=60)
                                       CCGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATC                                                         >  SRR4302533.1997771/1‑94 (MQ=60)
                                        CGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCTGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCC                                                        <  SRR4302533.3406422/94‑1 (MQ=60)
                                        CGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCC                                                        <  SRR4302533.5614144/94‑1 (MQ=60)
                                        CGGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCC                                                        >  SRR4302533.3545686/1‑94 (MQ=60)
                                        CGGATAAGGCGTTTACGACGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCC                                                        >  SRR4302533.2715622/1‑94 (MQ=60)
                                         GGATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCG                                                       >  SRR4302533.5015830/1‑94 (MQ=60)
                                          GATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGG                                                      <  SRR4302533.4702067/94‑1 (MQ=60)
                                          GATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGG                                                      <  SRR4302533.2993870/94‑1 (MQ=60)
                                          GATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGG                                                      <  SRR4302533.4607486/94‑1 (MQ=60)
                                           ATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC                                                     >  SRR4302533.4539023/1‑94 (MQ=60)
                                           ATAAGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC                                                     >  SRR4302533.791053/1‑94 (MQ=60)
                                              AGGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATA                                                  <  SRR4302533.4565127/94‑1 (MQ=60)
                                               GGCGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAA                                                 >  SRR4302533.1046064/1‑94 (MQ=60)
                                                 CGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                                               <  SRR4302533.5176639/94‑1 (MQ=60)
                                                 CGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                                               >  SRR4302533.4394225/1‑94 (MQ=60)
                                                 CGTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                                               >  SRR4302533.5626222/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              <  SRR4302533.2755866/94‑1 (MQ=60)
                                                  GTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302533.5129583/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302533.5425646/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              <  SRR4302533.2730852/94‑1 (MQ=60)
                                                    TTACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                                            >  SRR4302533.1880284/1‑94 (MQ=60)
                                                     TACGCCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATACACAACG                                           >  SRR4302533.5045511/1‑94 (MQ=60)
                                                         CCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCT                                       >  SRR4302533.3116427/1‑94 (MQ=60)
                                                         CCGCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCT                                       >  SRR4302533.1586402/1‑94 (MQ=60)
                                                           GCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTC                                     <  SRR4302533.4729486/94‑1 (MQ=60)
                                                           GCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTC                                     >  SRR4302533.4863492/1‑94 (MQ=60)
                                                           GCATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTC                                     >  SRR4302533.4332007/1‑94 (MQ=60)
                                                             ATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG                                   <  SRR4302533.6116562/94‑1 (MQ=60)
                                                             ATCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG                                   <  SRR4302533.2618000/94‑1 (MQ=60)
                                                              TCAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGC                                  >  SRR4302533.2713453/1‑94 (MQ=60)
                                                               CAGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGCGCTTCAGCT                                 >  SRR4302533.6026768/1‑94 (MQ=60)
                                                                AGG‑CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTA                                <  SRR4302533.5305919/94‑1 (MQ=60)
                                                                  GGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAA                               >  SRR4302533.1739722/1‑94 (MQ=60)
                                                                  GGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAA                               >  SRR4302533.3066413/1‑94 (MQ=60)
                                                                  GGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAA                               >  SRR4302533.2895704/1‑94 (MQ=60)
                                                                   GCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAA                              <  SRR4302533.3347753/94‑1 (MQ=60)
                                                                    CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAC                             <  SRR4302533.65199/94‑1 (MQ=60)
                                                                     ATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACA                            >  SRR4302533.3316331/1‑94 (MQ=60)
                                                                     ATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACA                            <  SRR4302533.4060169/94‑1 (MQ=60)
                                                                     ATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACA                            >  SRR4302533.4437636/1‑94 (MQ=60)
                                                                      TCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAG                           <  SRR4302533.4023961/94‑1 (MQ=60)
                                                                      TCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAG                           <  SRR4302533.4734923/94‑1 (MQ=60)
                                                                      TCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAG                           >  SRR4302533.1911575/1‑94 (MQ=60)
                                                                       CTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC                          >  SRR4302533.816177/1‑94 (MQ=60)
                                                                       CTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC                          >  SRR4302533.3416029/1‑94 (MQ=60)
                                                                       CTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC                          >  SRR4302533.3996373/1‑94 (MQ=60)
                                                                       CTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC                          >  SRR4302533.2533942/1‑94 (MQ=60)
                                                                        TGAGCATGATGCCTGAGGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCC                         <  SRR4302533.5095290/94‑1 (MQ=60)
                                                                         TAGCATGATGCTTGATGAGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCC                        <  SRR4302533.1141211/94‑1 (MQ=60)
                                                                         GAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCC                        >  SRR4302533.3820282/1‑94 (MQ=60)
                                                                          AGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCA                       <  SRR4302533.4637146/94‑1 (MQ=60)
                                                                          AGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCA                       >  SRR4302533.222195/1‑94 (MQ=60)
                                                                          AGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCA                       <  SRR4302533.2735723/94‑1 (MQ=60)
                                                                          AGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCACA                       >  SRR4302533.4313697/1‑94 (MQ=60)
                                                                           GCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAA                      >  SRR4302533.3756491/1‑94 (MQ=60)
                                                                            CATGATGCCTGATGCGGCGCTGGTCTCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAT                     >  SRR4302533.3538121/1‑94 (MQ=60)
                                                                            CATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAT                     >  SRR4302533.5632380/1‑94 (MQ=60)
                                                                            CATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAT                     >  SRR4302533.790437/1‑94 (MQ=60)
                                                                            CATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAT                     >  SRR4302533.1714226/1‑94 (MQ=60)
                                                                             ATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATC                    <  SRR4302533.3405726/94‑1 (MQ=60)
                                                                             ATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATC                    >  SRR4302533.1932871/1‑94 (MQ=60)
                                                                             ATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATC                    >  SRR4302533.2373623/1‑94 (MQ=60)
                                                                                ATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATT                 <  SRR4302533.5607444/94‑1 (MQ=60)
                                                                                ATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATT                 <  SRR4302533.534102/94‑1 (MQ=60)
                                                                                   CCTGATGCGGCGCTGGTCGCGTCTTATCAGTCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              <  SRR4302533.447883/94‑1 (MQ=60)
                                                                                   CCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              <  SRR4302533.45948/94‑1 (MQ=60)
                                                                                    CTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGG             <  SRR4302533.6243974/94‑1 (MQ=60)
                                                                                    CTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGG             >  SRR4302533.4890601/1‑94 (MQ=60)
                                                                                        TGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACA         >  SRR4302533.3606829/1‑94 (MQ=60)
                                                                                         GCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTGCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACAC        >  SRR4302533.840031/1‑94 (MQ=60)
                                                                                          CGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACA       <  SRR4302533.5972133/94‑1 (MQ=60)
                                                                                          CGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACA       >  SRR4302533.5144449/1‑94 (MQ=60)
                                                                                           GGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAG      >  SRR4302533.2478425/1‑94 (MQ=60)
                                                                                            GCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGC     >  SRR4302533.569525/1‑94 (MQ=60)
                                                                                              C‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  >  SRR4302533.134660/1‑94 (MQ=60)
                                                                                              C‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  <  SRR4302533.3788693/94‑1 (MQ=60)
                                                                                              |                                                                                              
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  >  NC_000913/4051076‑4051263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑