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 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,639,116 C→A 88.1% A523A (GCG→GCT)  pbpC ← penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,639,1160CA88.1% 324.5 / 46.4 144A523A (GCG→GCT) pbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base A (67/60);  total (75/69)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 144 total aligned reads displayed.

CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCC  >  CP009273/2638972‑2639261
                                                                                                                                                |                                                                                                                                                 
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                              <  2:3238796/148‑1 (MQ=255)
    ttCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAgg                                                                                                                                            >  1:3231331/1‑148 (MQ=255)
          cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                      >  1:1520976/1‑148 (MQ=255)
            gcgcGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGa                                                                                                                                    >  2:488222/1‑148 (MQ=255)
            gcgcGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGa                                                                                                                                    <  2:1311098/148‑1 (MQ=255)
             cgcgACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGAt                                                                                                                                   <  1:1724211/148‑1 (MQ=255)
               cgACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAg                                                                                                                                 <  2:3238513/148‑1 (MQ=255)
                gACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGc                                                                                                                                <  2:3099793/148‑1 (MQ=255)
                 aCAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGcc                                                                                                                               <  2:1049675/148‑1 (MQ=255)
                  cAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCa                                                                                                                              <  2:1092678/148‑1 (MQ=255)
                    gTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATa                                                                                                                            >  2:2721259/1‑148 (MQ=255)
                       aGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCt                                                                                                                         >  2:3152325/1‑148 (MQ=255)
                        gATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTg                                                                                                                        >  1:566834/1‑148 (MQ=255)
                        gATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTg                                                                                                                        >  1:566845/1‑148 (MQ=255)
                            ttattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGc                                                                                                                    >  2:1547670/1‑148 (MQ=255)
                             tattGACCTGATTCAACAATGGTATGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGcc                                                                                                                   >  1:588514/1‑148 (MQ=255)
                             tattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGcc                                                                                                                   <  1:1960023/148‑1 (MQ=255)
                             tattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGcc                                                                                                                   >  2:3270686/1‑148 (MQ=255)
                             tattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGcc                                                                                                                   <  2:3126158/148‑1 (MQ=255)
                             tattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGcc                                                                                                                   >  2:1475112/1‑148 (MQ=255)
                              attGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGccc                                                                                                                  >  1:103107/1‑148 (MQ=255)
                               ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                 >  1:2279821/1‑148 (MQ=255)
                                  aCCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGttt                                                                                                              >  2:439727/1‑148 (MQ=255)
                                  aCCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGttt                                                                                                              >  2:439547/1‑148 (MQ=255)
                                  aCCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGttt                                                                                                              >  1:481278/1‑148 (MQ=255)
                                   ccTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGtttt                                                                                                             >  2:53116/1‑148 (MQ=255)
                                    cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATCGCCATAGCTGGTGCCCGTTTTc                                                                                                            <  2:1366185/148‑1 (MQ=255)
                                    cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                            <  2:2389663/148‑1 (MQ=255)
                                    cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                            <  2:733027/148‑1 (MQ=255)
                                    cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                            <  2:2768146/148‑1 (MQ=255)
                                     tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCTAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                           >  2:2399351/1‑148 (MQ=255)
                                     tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                           >  1:3126603/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                       >  2:1530787/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                       >  1:1502718/1‑148 (MQ=255)
                                          caacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGcc                                                                                                      >  2:530092/1‑148 (MQ=255)
                                           aacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCa                                                                                                     <  1:1820956/148‑1 (MQ=255)
                                            acaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCATTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAg                                                                                                    >  1:2375075/1‑148 (MQ=255)
                                            acaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAg                                                                                                    <  2:1172627/148‑1 (MQ=255)
                                             caaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGt                                                                                                   <  2:2886002/148‑1 (MQ=255)
                                             caaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGt                                                                                                   <  1:3270686/148‑1 (MQ=255)
                                                      ggCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcg                                                                                          <  1:2669585/148‑1 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                        >  2:2267380/1‑148 (MQ=255)
                                                            ggCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGc                                                                                    >  1:1735562/1‑148 (MQ=255)
                                                             gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                                                                   >  2:2481274/1‑148 (MQ=255)
                                                             gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                                                                   <  2:1090025/148‑1 (MQ=255)
                                                             gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                                                                   <  2:93800/148‑1 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                  <  2:588464/148‑1 (MQ=255)
                                                                  gCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCa                                                                              <  1:2675917/148‑1 (MQ=255)
                                                                  gCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCa                                                                              >  2:3291089/1‑148 (MQ=255)
                                                                  gCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCa                                                                              >  2:3291080/1‑148 (MQ=255)
                                                                   ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                             >  1:770504/1‑148 (MQ=255)
                                                                   ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                             <  2:3080510/148‑1 (MQ=255)
                                                                   ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                             <  2:3080139/148‑1 (MQ=255)
                                                                      aaCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAt                                                                          <  1:2483693/148‑1 (MQ=255)
                                                                       aCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATc                                                                         >  2:2867702/1‑148 (MQ=255)
                                                                        cTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATcc                                                                        >  2:2244191/1‑148 (MQ=255)
                                                                              aacaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAg                                                                  <  2:2370490/148‑1 (MQ=255)
                                                                              aacaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAg                                                                  >  2:623638/1‑148 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                 >  2:2595706/1‑148 (MQ=255)
                                                                                  acGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGtt                                                              >  1:1044375/1‑148 (MQ=255)
                                                                                      gCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcgc                                                          <  1:1036590/148‑1 (MQ=255)
                                                                                         tGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtc                                                       <  2:858534/148‑1 (MQ=255)
                                                                                          gCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtca                                                      <  2:2770766/148‑1 (MQ=255)
                                                                                           ccGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcat                                                     >  2:493989/1‑148 (MQ=255)
                                                                                             gTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtcatca                                                   <  2:2151859/148‑1 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                 <  2:1478977/148‑1 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                 >  1:2607716/1‑148 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                 >  1:2607710/1‑148 (MQ=255)
                                                                                                  gTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCAt                                              <  2:1502758/148‑1 (MQ=255)
                                                                                                  gTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCAt                                              >  1:2267633/1‑148 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                            <  2:2279821/148‑1 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                            >  2:1096451/1‑148 (MQ=255)
                                                                                                       ccAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATcc                                         >  2:806852/1‑148 (MQ=255)
                                                                                                         aGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCgg                                       <  1:1980781/148‑1 (MQ=255)
                                                                                                            ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                    <  2:2678041/148‑1 (MQ=255)
                                                                                                             cAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCAGGCGaa                                   >  1:2567001/1‑148 (MQ=255)
                                                                                                               gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                 <  1:488211/148‑1 (MQ=255)
                                                                                                               gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                 >  1:1131171/1‑148 (MQ=255)
                                                                                                               gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                 >  1:1131175/1‑148 (MQ=255)
                                                                                                               gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                 >  1:1370625/1‑148 (MQ=255)
                                                                                                               gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                 >  1:1370815/1‑148 (MQ=255)
                                                                                                                  cccAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATc                              >  2:938161/1‑148 (MQ=255)
                                                                                                                    cAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCa                            <  1:513175/148‑1 (MQ=255)
                                                                                                                    cAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCa                            <  1:513202/148‑1 (MQ=255)
                                                                                                                      aTGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcg                          <  2:1034260/148‑1 (MQ=255)
                                                                                                                       tGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgc                         >  1:2366900/1‑148 (MQ=255)
                                                                                                                          cATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgccgc                      >  2:3154858/1‑148 (MQ=255)
                                                                                                                            tAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccc                    <  2:3106515/148‑1 (MQ=255)
                                                                                                                            tAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccc                    >  2:507519/1‑148 (MQ=255)
                                                                                                                            tAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccc                    >  2:507535/1‑148 (MQ=255)
                                                                                                                             agcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                   >  2:3032842/1‑148 (MQ=255)
                                                                                                                             agcgagcgTTAACCCCAATAGCCCAGGCGTCACGAAAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGGTCCACGCCGcccc                   >  2:1387533/1‑148 (MQ=255)
                                                                                                                              gcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccccc                  <  2:2989564/148‑1 (MQ=255)
                                                                                                                                   cgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacga             <  2:2375075/148‑1 (MQ=255)
                                                                                                                                     ttAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACa           >  1:1518974/1‑148 (MQ=255)
                                                                                                                                     ttAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACa           >  2:2861278/1‑148 (MQ=255)
                                                                                                                                            cAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg    >  1:1466748/1‑148 (MQ=255)
                                                                                                                                            cAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg    <  2:2602104/148‑1 (MQ=255)
                                                                                                                                            cAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg    >  1:1466761/1‑148 (MQ=255)
                                                                                                                                              aTCGCCCAGGCGACACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGcc  >  1:512709/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCC  >  CP009273/2638972‑2639261

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
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  >  CP009273/2638972‑2639264
                                                                                                                                                |                                                                                                                                                    
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCA                                                                                                                                                 <  SRR19973070.1786965/150‑1 (MQ=60)
  GATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGG                                                                                                                                               >  SRR19973070.3252196/1‑150 (MQ=60)
      CGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTC                                                                                                                                           >  SRR19973070.2411305/1‑150 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                         >  SRR19973070.1531651/1‑150 (MQ=60)
             CGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAG                                                                                                                                    <  SRR19973070.1736148/150‑1 (MQ=60)
                   AGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAG                                                                                                                              >  SRR19973070.1406469/1‑150 (MQ=60)
                   AGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAG                                                                                                                              <  SRR19973070.3317064/150‑1 (MQ=60)
                      AAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTG                                                                                                                           >  SRR19973070.570642/1‑150 (MQ=60)
                      AAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTG                                                                                                                           >  SRR19973070.570653/1‑150 (MQ=60)
                           ATTATTGACCTGATTCAACAATGGTATGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCC                                                                                                                      >  SRR19973070.592530/1‑150 (MQ=60)
                           ATTATTGACCTGATTCAACAATGGTATGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCC                                                                                                                      >  SRR19973070.592580/1‑150 (MQ=60)
                            TTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCC                                                                                                                     >  SRR19973070.104087/1‑150 (MQ=60)
                             TATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCG                                                                                                                    <  SRR19973070.1973231/150‑1 (MQ=60)
                             TATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCG                                                                                                                    <  SRR19973070.1973228/150‑1 (MQ=60)
                             TATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCG                                                                                                                    >  SRR19973070.2294966/1‑150 (MQ=60)
                                TGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTT                                                                                                                 >  SRR19973070.484673/1‑150 (MQ=60)
                                   CCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCC                                                                                                              >  SRR19973070.3146787/1‑150 (MQ=60)
                                       ATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGC                                                                                                          >  SRR19973070.1513281/1‑150 (MQ=60)
                                       ATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGC                                                                                                          >  SRR19973070.1513241/1‑150 (MQ=60)
                                         TCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCA                                                                                                        <  SRR19973070.3158778/150‑1 (MQ=60)
                                          CAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCATTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAG                                                                                                       >  SRR19973070.2390774/1‑150 (MQ=60)
                                           AACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGT                                                                                                      <  SRR19973070.1833496/150‑1 (MQ=60)
                                             CAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGG                                                                                                    <  SRR19973070.3291719/150‑1 (MQ=60)
                                                      GGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCG                                                                                           <  SRR19973070.2686926/150‑1 (MQ=60)
                                                      GGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCG                                                                                           <  SRR19973070.550213/150‑1 (MQ=60)
                                                          CTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGC                                                                                       >  SRR19973070.1747591/1‑150 (MQ=60)
                                                           TGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCA                                                                                      >  SRR19973070.974903/1‑150 (MQ=60)
                                                                 AGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAC                                                                                >  SRR19973070.775908/1‑150 (MQ=60)
                                                                  GCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACT                                                                               <  SRR19973070.2693291/150‑1 (MQ=60)
                                                                    CAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAT                                                                             <  SRR19973070.2686545/150‑1 (MQ=60)
                                                                      AACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCC                                                                           <  SRR19973070.2499842/150‑1 (MQ=60)
                                                                             CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                    >  SRR19973070.113929/1‑150 (MQ=60)
                                                                                CAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTT                                                                 >  SRR19973070.1051738/1‑150 (MQ=60)
                                                                                      GCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTT                                                           <  SRR19973070.1043913/150‑1 (MQ=60)
                                                                                             GTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGC                                                    >  SRR19973070.2624733/1‑150 (MQ=60)
                                                                                             GTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGC                                                    >  SRR19973070.2624739/1‑150 (MQ=60)
                                                                                                CGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCAT                                                 >  SRR19973070.2282691/1‑150 (MQ=60)
                                                                                                         AGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCG                                        <  SRR19973070.1994180/150‑1 (MQ=60)
                                                                                                           TCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCAGGCGAA                                      >  SRR19973070.2583778/1‑150 (MQ=60)
                                                                                                            CCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAAT                                     <  SRR19973070.2306213/150‑1 (MQ=60)
                                                                                                             CAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATG                                    >  SRR19973070.1138978/1‑150 (MQ=60)
                                                                                                             CAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATG                                    >  SRR19973070.1138982/1‑150 (MQ=60)
                                                                                                             CAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATG                                    >  SRR19973070.1380140/1‑150 (MQ=60)
                                                                                                             CAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATG                                    >  SRR19973070.1379949/1‑150 (MQ=60)
                                                                                                               GATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGAT                                  <  SRR19973070.491638/150‑1 (MQ=60)
                                                                                                               GATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGAT                                  <  SRR19973070.491649/150‑1 (MQ=60)
                                                                                                                ATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTACCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATC                                 <  SRR19973070.1245244/150‑1 (MQ=60)
                                                                                                                    CAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACG                             <  SRR19973070.516724/150‑1 (MQ=60)
                                                                                                                    CAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACG                             <  SRR19973070.516697/150‑1 (MQ=60)
                                                                                                                     AATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGC                            >  SRR19973070.2382553/1‑150 (MQ=60)
                                                                                                                                   CGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACA              >  SRR19973070.1529625/1‑150 (MQ=60)
                                                                                                                                       AACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAA          >  SRR19973070.2583988/1‑150 (MQ=60)
                                                                                                                                          CCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGG       >  SRR19973070.1476998/1‑150 (MQ=60)
                                                                                                                                          CCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGG       >  SRR19973070.1477011/1‑150 (MQ=60)
                                                                                                                                          CCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAAACCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGG       >  SRR19973070.1655951/1‑150 (MQ=60)
                                                                                                                                            CAATCGCCCAGGCGACACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCC     >  SRR19973070.516228/1‑150 (MQ=60)
                                                                                                                                               TAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  <  SRR19973070.3252454/150‑1 (MQ=60)
                                                                                                                                               TAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  <  SRR19973070.3252409/150‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                                    
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  >  CP009273/2638972‑2639264

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: ‑