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 F78 I0 R1 626 241.3 10389380 77.7% 8072548 149.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,354,405 C→T 92.7% intergenic (+430/‑130) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,354,4050CT92.7% 561.3 / 23.6 177intergenic (+430/‑130)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (77/87);  total (83/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 179 total aligned reads displayed.

AATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAA  >  CP009273/3354259‑3354551
                                                                                                                                                  |                                                                                                                                                  
aaTATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCAtt                                                                                                                                                 <  2:524585/150‑1 (MQ=255)
aaTATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCAtt                                                                                                                                                 <  2:524565/150‑1 (MQ=255)
   atGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATtaaa                                                                                                                                              >  2:2416498/1‑150 (MQ=255)
    tGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATtaaat                                                                                                                                             >  1:1849097/1‑150 (MQ=255)
    tGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATtaaat                                                                                                                                             >  1:2026046/1‑150 (MQ=255)
       atatAAAAAACCAATATATTATTCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATtaaataaa                                                                                                                                          <  1:103113/150‑1 (MQ=255)
        tataAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTaaataaat                                                                                                                                         <  1:1053422/150‑1 (MQ=255)
           aaaaaaCCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACa                                                                                                                                      <  2:1248683/150‑1 (MQ=255)
           aaaaaaCCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACa                                                                                                                                      <  2:1763248/150‑1 (MQ=255)
            aaaaaCCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACAt                                                                                                                                     >  1:2486052/1‑150 (MQ=255)
              aaaCCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATcc                                                                                                                                   >  1:1636691/1‑150 (MQ=255)
                aCCAATATATTATCCTTAATTCTCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGAAACTATAAAATTTATTAAATAAATACATACCa                                                                                                                                 >  2:3938756/1‑150 (MQ=255)
                aCCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCa                                                                                                                                 >  2:1406588/1‑150 (MQ=255)
                 ccAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAg                                                                                                                                <  2:4563107/150‑1 (MQ=255)
                      atatTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGt                                                                                                                           <  1:203678/150‑1 (MQ=255)
                      atatTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGt                                                                                                                           <  1:2806870/150‑1 (MQ=255)
                       tattatCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGtt                                                                                                                          <  1:3960790/150‑1 (MQ=255)
                         ttatCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTgg                                                                                                                        >  1:2907851/1‑150 (MQ=255)
                         ttatCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTgg                                                                                                                        <  2:3868263/150‑1 (MQ=255)
                            tCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTa                                                                                                                     <  2:3114675/150‑1 (MQ=255)
                            tCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTa                                                                                                                     <  2:2877604/150‑1 (MQ=255)
                             ccTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTAt                                                                                                                    <  1:3936888/150‑1 (MQ=255)
                             ccTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTAt                                                                                                                    >  2:3850344/1‑150 (MQ=255)
                                tAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTc                                                                                                                 <  2:3425699/150‑1 (MQ=255)
                                 aaTTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCt                                                                                                                <  2:1511544/150‑1 (MQ=255)
                                  aTTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTg                                                                                                               <  1:3010825/150‑1 (MQ=255)
                                  aTTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTg                                                                                                               <  2:2911164/150‑1 (MQ=255)
                                     aTCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAAc                                                                                                            <  1:4525776/150‑1 (MQ=255)
                                     aTCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAAc                                                                                                            <  1:4062687/150‑1 (MQ=255)
                                      tCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACg                                                                                                           <  2:3725277/150‑1 (MQ=255)
                                       cTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGa                                                                                                          >  1:712190/1‑150 (MQ=255)
                                       cTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGa                                                                                                          >  1:1550990/1‑150 (MQ=255)
                                          tatatTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGtt                                                                                                       <  1:35453/150‑1 (MQ=255)
                                            tatTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCACGTTGGCTATTTCTGAACGAGttat                                                                                                     >  2:4815271/1‑150 (MQ=255)
                                             atTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCACGTTGGCTATTTCTGAACGAGttatt                                                                                                    >  2:3554189/1‑150 (MQ=255)
                                               tttCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGttattta                                                                                                  <  1:849125/150‑1 (MQ=255)
                                               tttCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGttattta                                                                                                  <  1:4843510/150‑1 (MQ=255)
                                                     aTTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACa                                                                                            >  2:4268072/1‑150 (MQ=255)
                                                     aTTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACa                                                                                            >  2:343741/1‑150 (MQ=255)
                                                         ttCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTca                                                                                        <  1:3736148/150‑1 (MQ=255)
                                                             gcgcAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACaaa                                                                                    >  2:776986/1‑150 (MQ=255)
                                                              cgcAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAAc                                                                                   >  2:3675378/1‑150 (MQ=255)
                                                               gcAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACa                                                                                  <  1:3378974/150‑1 (MQ=255)
                                                               gcAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACa                                                                                  <  1:774799/150‑1 (MQ=255)
                                                               gcAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACa                                                                                  >  2:1749927/1‑150 (MQ=255)
                                                                  aaCATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGat                                                                               <  2:2843678/150‑1 (MQ=255)
                                                                    cATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGatat                                                                             <  2:4934202/150‑1 (MQ=255)
                                                                    cATTCATGTAATGAATAATTATCACCACAAACAACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGatat                                                                             <  2:4437486/150‑1 (MQ=255)
                                                                      ttCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGatatat                                                                           <  2:1192766/150‑1 (MQ=255)
                                                                          tGTAATGAATAATTATCACCACACACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATAc                                                                       <  1:3926464/150‑1 (MQ=255)
                                                                           gTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACt                                                                      <  2:1732067/150‑1 (MQ=255)
                                                                           gTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACt                                                                      >  2:3156133/1‑150 (MQ=255)
                                                                            tAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTg                                                                     >  1:3681485/1‑150 (MQ=255)
                                                                               tGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCtt                                                                  <  1:1968510/150‑1 (MQ=255)
                                                                                gAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTc                                                                 <  2:2519815/150‑1 (MQ=255)
                                                                                gAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTc                                                                 >  1:582962/1‑150 (MQ=255)
                                                                                gAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTc                                                                 >  1:583138/1‑150 (MQ=255)
                                                                                 aataatTATCACCACAAACCACCTCACCATCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCa                                                                <  1:2525122/150‑1 (MQ=255)
                                                                                  ataatTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAg                                                               <  1:3929515/150‑1 (MQ=255)
                                                                                    aatTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCt                                                             <  1:2761778/150‑1 (MQ=255)
                                                                                      ttATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCttt                                                           <  2:2425407/150‑1 (MQ=255)
                                                                                       tATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCtttt                                                          <  1:4712032/150‑1 (MQ=255)
                                                                                          caccacAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttat                                                       >  2:2200213/1‑150 (MQ=255)
                                                                                          caccacAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttat                                                       >  2:2200177/1‑150 (MQ=255)
                                                                                           accacAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttatt                                                      >  2:1280819/1‑150 (MQ=255)
                                                                                           accacAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttatt                                                      <  1:149741/150‑1 (MQ=255)
                                                                                            ccacAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttattt                                                     <  1:216143/150‑1 (MQ=255)
                                                                                             cacaAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTttttatttt                                                    >  1:368294/1‑150 (MQ=255)
                                                                                               caAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaa                                                  >  2:4938176/1‑150 (MQ=255)
                                                                                               caAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaa                                                  >  1:1195771/1‑150 (MQ=255)
                                                                                                  accaccTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaa                                               >  2:709221/1‑150 (MQ=255)
                                                                                                  accaccTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTaaaaa                                               >  2:5035388/1‑150 (MQ=255)
                                                                                                     accTCACCAGCACCTTATTATTCCCATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaata                                            <  1:3596814/150‑1 (MQ=255)
                                                                                                      ccTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataa                                           >  1:350792/1‑150 (MQ=255)
                                                                                                       cTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAaaataaa                                          >  1:4858070/1‑150 (MQ=255)
                                                                                                         cACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAAtt                                        <  1:506984/150‑1 (MQ=255)
                                                                                                         cACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAAtt                                        <  1:4815271/150‑1 (MQ=255)
                                                                                                          aCCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAAttt                                       >  1:867803/1‑150 (MQ=255)
                                                                                                           ccAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTg                                      >  1:2831794/1‑150 (MQ=255)
                                                                                                            cAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGAGATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGt                                     >  1:3184361/1‑150 (MQ=255)
                                                                                                             aGCACCGTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCACGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGtt                                    <  2:3130381/150‑1 (MQ=255)
                                                                                                                 ccTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGttatta                                >  1:1961685/1‑150 (MQ=255)
                                                                                                                  cTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTaa                               <  1:4736874/150‑1 (MQ=255)
                                                                                                                   ttattaTTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAg                              <  2:4340504/150‑1 (MQ=255)
                                                                                                                        attCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGatat                         <  2:354071/150‑1 (MQ=255)
                                                                                                                             caTGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                    <  1:895584/150‑1 (MQ=255)
                                                                                                                             caTGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                    >  1:2191622/1‑150 (MQ=255)
                                                                                                                             caTGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCa                    <  1:2454363/150‑1 (MQ=255)
                                                                                                                               tGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCata                  >  2:1284132/1‑150 (MQ=255)
                                                                                                                                   cGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtt              >  2:2028086/1‑150 (MQ=255)
                                                                                                                                   cGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtt              >  2:150728/1‑150 (MQ=255)
                                                                                                                                   cGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtt              <  2:4869381/150‑1 (MQ=255)
                                                                                                                                   cGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGtt              <  2:2486052/150‑1 (MQ=255)
                                                                                                                                    gTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGttt             <  2:4441844/150‑1 (MQ=255)
                                                                                                                                      aaCTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGttttt           <  1:4432543/150‑1 (MQ=255)
                                                                                                                                      aaCTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTGTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGttttt           <  1:1833789/150‑1 (MQ=255)
                                                                                                                                        cTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTAGTAAGGATATGTTCATATGTTTTTCa         >  2:2616766/1‑150 (MQ=255)
                                                                                                                                           taAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaa      >  2:3955340/1‑150 (MQ=255)
                                                                                                                                            aaaaTTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCaaaaa     >  2:4905838/1‑150 (MQ=255)
                                                                                                                                               aTTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGaa  <  2:624509/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTTCAAAAAGAA  >  CP009273/3354259‑3354551

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTT  >  CP009273/3354256‑3354542
                                                                                                                                                     |                                                                                                                                         
CAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATTATCGTAACTATAAAATTT                                                                                                                                           >  SRR19973083.424487/1‑150 (MQ=60)
    ATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTA                                                                                                                                       >  SRR19973083.734221/1‑150 (MQ=60)
    ATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTA                                                                                                                                       >  SRR19973083.2084507/1‑150 (MQ=60)
      ATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAA                                                                                                                                     >  SRR19973083.2981302/1‑150 (MQ=60)
      ATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAA                                                                                                                                     >  SRR19973083.355752/1‑150 (MQ=60)
       TGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAAT                                                                                                                                    >  SRR19973083.1858422/1‑150 (MQ=60)
       TGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAAT                                                                                                                                    >  SRR19973083.2036465/1‑150 (MQ=60)
       TGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAAT                                                                                                                                    >  SRR19973083.2036393/1‑150 (MQ=60)
          ATATAAAAAACCAATATATTATTCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAA                                                                                                                                 <  SRR19973083.103796/150‑1 (MQ=60)
           TATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAAT                                                                                                                                <  SRR19973083.1058727/150‑1 (MQ=60)
             TAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATAC                                                                                                                              <  SRR19973083.323099/150‑1 (MQ=60)
               AAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACAT                                                                                                                            >  SRR19973083.2498750/1‑150 (MQ=60)
                 AAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCC                                                                                                                          >  SRR19973083.1644841/1‑150 (MQ=60)
                   ACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCA                                                                                                                        >  SRR19973083.3835893/1‑150 (MQ=60)
                    CCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAG                                                                                                                       <  SRR19973083.2786843/150‑1 (MQ=60)
                    CCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAG                                                                                                                       >  SRR19973083.4457928/1‑150 (MQ=60)
                         ATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGT                                                                                                                  <  SRR19973083.204815/150‑1 (MQ=60)
                         ATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGT                                                                                                                  <  SRR19973083.2820792/150‑1 (MQ=60)
                          TATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTT                                                                                                                 <  SRR19973083.3978428/150‑1 (MQ=60)
                            TTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGG                                                                                                               >  SRR19973083.2922086/1‑150 (MQ=60)
                                CCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTAT                                                                                                           <  SRR19973083.3954345/150‑1 (MQ=60)
                                CCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTAT                                                                                                           <  SRR19973083.3954413/150‑1 (MQ=60)
                                CCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTAT                                                                                                           <  SRR19973083.3954418/150‑1 (MQ=60)
                                     ATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTG                                                                                                      <  SRR19973083.3025347/150‑1 (MQ=60)
                                        ATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAAC                                                                                                   <  SRR19973083.4080718/150‑1 (MQ=60)
                                        ATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAAC                                                                                                   <  SRR19973083.4545409/150‑1 (MQ=60)
                                          CTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGA                                                                                                 >  SRR19973083.1558791/1‑150 (MQ=60)
                                          CTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGA                                                                                                 >  SRR19973083.715728/1‑150 (MQ=60)
                                             TATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTT                                                                                              <  SRR19973083.35727/150‑1 (MQ=60)
                                                  TTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTA                                                                                         <  SRR19973083.853132/150‑1 (MQ=60)
                                                  TTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTA                                                                                         <  SRR19973083.4864301/150‑1 (MQ=60)
                                                  TTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTA                                                                                         <  SRR19973083.4864341/150‑1 (MQ=60)
                                                       TATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTAC                                                                                    >  SRR19973083.874363/1‑150 (MQ=60)
                                                            TTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCA                                                                               <  SRR19973083.3753102/150‑1 (MQ=60)
                                                            TTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCA                                                                               >  SRR19973083.1140404/1‑150 (MQ=60)
                                                                 CGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAAC                                                                          <  SRR19973083.1658918/150‑1 (MQ=60)
                                                                  GCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACA                                                                         <  SRR19973083.778577/150‑1 (MQ=60)
                                                                  GCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACA                                                                         <  SRR19973083.3394647/150‑1 (MQ=60)
                                                                     AACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGAT                                                                      >  SRR19973083.3656217/1‑150 (MQ=60)
                                                                       CATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATAT                                                                    >  SRR19973083.2833284/1‑150 (MQ=60)
                                                                       CATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATAT                                                                    >  SRR19973083.3573216/1‑150 (MQ=60)
                                                                             TGTAATGAATAATTATCACCACACACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATAC                                                              <  SRR19973083.3943944/150‑1 (MQ=60)
                                                                             TGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATAC                                                              <  SRR19973083.1845053/150‑1 (MQ=60)
                                                                               TAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTG                                                            >  SRR19973083.3698278/1‑150 (MQ=60)
                                                                               TAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTG                                                            >  SRR19973083.3698332/1‑150 (MQ=60)
                                                                                  TGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTT                                                         <  SRR19973083.1978576/150‑1 (MQ=60)
                                                                                   GAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTC                                                        >  SRR19973083.585890/1‑150 (MQ=60)
                                                                                   GAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTC                                                        >  SRR19973083.586055/1‑150 (MQ=60)
                                                                                   GAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTC                                                        >  SRR19973083.586066/1‑150 (MQ=60)
                                                                                    AATAATTATCACCACAAACCACCTCACCATCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCA                                                       <  SRR19973083.2537943/150‑1 (MQ=60)
                                                                                     ATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAG                                                      <  SRR19973083.3946901/150‑1 (MQ=60)
                                                                                     ATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAG                                                      <  SRR19973083.3947003/150‑1 (MQ=60)
                                                                                       AATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCT                                                    <  SRR19973083.2775489/150‑1 (MQ=60)
                                                                                       AATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCT                                                    <  SRR19973083.2775519/150‑1 (MQ=60)
                                                                                        ATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTT                                                   >  SRR19973083.793833/1‑150 (MQ=60)
                                                                                          TATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTT                                                 <  SRR19973083.4732196/150‑1 (MQ=60)
                                                                                            TCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTA                                               >  SRR19973083.3146552/1‑150 (MQ=60)
                                                                                              ACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATT                                             <  SRR19973083.150633/150‑1 (MQ=60)
                                                                                               CCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTT                                            <  SRR19973083.1543308/150‑1 (MQ=60)
                                                                                               CCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTT                                            <  SRR19973083.217323/150‑1 (MQ=60)
                                                                                                CACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTT                                           >  SRR19973083.370022/1‑150 (MQ=60)
                                                                                                  CAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAA                                         >  SRR19973083.1201682/1‑150 (MQ=60)
                                                                                                        ACCTCACCAGCACCTTATTATTCCCATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATA                                   <  SRR19973083.3613291/150‑1 (MQ=60)
                                                                                                         CCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAA                                  >  SRR19973083.352465/1‑150 (MQ=60)
                                                                                                          CTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAA                                 >  SRR19973083.4878925/1‑150 (MQ=60)
                                                                                                            CACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATT                               <  SRR19973083.509442/150‑1 (MQ=60)
                                                                                                            CACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATT                               <  SRR19973083.4835912/150‑1 (MQ=60)
                                                                                                             ACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTT                              >  SRR19973083.871854/1‑150 (MQ=60)
                                                                                                              CCAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTG                             >  SRR19973083.2845806/1‑150 (MQ=60)
                                                                                                               CAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGT                            >  SRR19973083.4669029/1‑150 (MQ=60)
                                                                                                               CAGCACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGAGATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGT                            >  SRR19973083.3199524/1‑150 (MQ=60)
                                                                                                                  CACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTAT                         >  SRR19973083.2217668/1‑150 (MQ=60)
                                                                                                                  CACCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTAT                         >  SRR19973083.508911/1‑150 (MQ=60)
                                                                                                                    CCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTA                       >  SRR19973083.1971708/1‑150 (MQ=60)
                                                                                                                    CCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTA                       >  SRR19973083.1971515/1‑150 (MQ=60)
                                                                                                                    CCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTTTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTA                       >  SRR19973083.1971723/1‑150 (MQ=60)
                                                                                                                    CCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTA                       <  SRR19973083.697168/150‑1 (MQ=60)
                                                                                                                    CCTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAAGACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTA                       <  SRR19973083.1395482/150‑1 (MQ=60)
                                                                                                                     CTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAA                      >  SRR19973083.2199304/1‑150 (MQ=60)
                                                                                                                     CTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAA                      <  SRR19973083.4757105/150‑1 (MQ=60)
                                                                                                                     CTTATTATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAA                      <  SRR19973083.4756301/150‑1 (MQ=60)
                                                                                                                           ATTCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATAT                <  SRR19973083.4286680/150‑1 (MQ=60)
                                                                                                                             TCACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGT              <  SRR19973083.2537302/150‑1 (MQ=60)
                                                                                                                              CACATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTT             >  SRR19973083.4859948/1‑150 (MQ=60)
                                                                                                                                CATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCA           <  SRR19973083.899833/150‑1 (MQ=60)
                                                                                                                                CATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCA           >  SRR19973083.2202609/1‑150 (MQ=60)
                                                                                                                                CATGATCGTAACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCA           <  SRR19973083.2466942/150‑1 (MQ=60)
                                                                                                                                         AACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTT  <  SRR19973083.4451828/150‑1 (MQ=60)
                                                                                                                                         AACTATAAAATTTATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTGTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTT  <  SRR19973083.1843028/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                         
CAAAATATGAATATAAAAAACCAATATATTATCCTTAATTATCTATATATTTTCCTATTTTTCAGCGCAAACATTCATGTAATGAATAATTATCACCACAAACCACCTCACCAGCACCTTATTATTCACATGATCGTAACTATAAAATTCATTAAATAAATACATCCCAGCAAGTTGGCTATTTCTGAACGAGTTATTTATTTACATTCACAAACAGATATATATACTGCTTCAGCTTTTTATTTTAAAAAATAAATTTGTTATTAAGGATATGTTCATATGTTTTT  >  CP009273/3354256‑3354542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: