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 F342 I0 R1 1331 225.9 7641038 95.8% 7320114 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,049,740 C→A 100% L47L (CTC→CTA)  fau → 5‑formyltetrahydrofolate cyclo‑ligase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,049,7400CA98.8% 493.8 / ‑1.7 163L47L (CTC→CTA) fau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (77/84);  total (77/86)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 165 total aligned reads displayed.

CCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAT  >  CP009273/3049593‑3049884
                                                                                                                                                   |                                                                                                                                                
ccGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTGATGGCGAACTc                                                                                                                                                  >  2:113282/1‑148 (MQ=255)
     aaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGTTGGCGAACTAGacac                                                                                                                                             >  1:1917797/1‑148 (MQ=255)
      aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAcc                                                                                                                                            <  1:1486609/148‑1 (MQ=255)
        tGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTACGCTATTTTTATGGCGAAAAAGACACCCa                                                                                                                                          >  1:3679429/1‑148 (MQ=255)
          aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCTCTCTTTTGATGGCGAACTAGACACCCAGc                                                                                                                                        >  2:1154015/1‑148 (MQ=255)
          aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGc                                                                                                                                        >  2:1153647/1‑148 (MQ=255)
          aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGc                                                                                                                                        <  1:2461136/148‑1 (MQ=255)
              ggCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACt                                                                                                                                    <  2:791265/148‑1 (MQ=255)
              ggCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGACGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACt                                                                                                                                    <  2:791266/148‑1 (MQ=255)
                    gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGa                                                                                                                              <  1:3407756/148‑1 (MQ=255)
                    gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGa                                                                                                                              <  2:3634107/148‑1 (MQ=255)
                     tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTCCTCTCTTTTTTTTGCGAACTAGACACCCACCcacacaaataa                                                                                                                             >  1:1811178/1‑143 (MQ=255)
                      cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaac                                                                                                                            >  2:2046116/1‑148 (MQ=255)
                        tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaacaa                                                                                                                          <  2:3014708/148‑1 (MQ=255)
                        tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaacaa                                                                                                                          <  2:3014712/148‑1 (MQ=255)
                        tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaacaa                                                                                                                          <  2:3304/148‑1 (MQ=255)
                            gCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTGCTCTATTTTTTTGGCGTACTAGAAACCCACCCACACAAATaaaaacaca                                                                                                                      >  1:1297090/1‑145 (MQ=255)
                              gTTACCGCCGGAACAACAGCAGGAATTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTgg                                                                                                                    <  1:3210638/148‑1 (MQ=255)
                                  aCGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcgcg                                                                                                                <  1:3672412/148‑1 (MQ=255)
                                   cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATTGCGAACTAGACACCCAGCCCCACATAGAACAACACAGgcgcgc                                                                                                               >  1:3028558/1‑148 (MQ=255)
                                   cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACAAGACACCCACCCACACAAAGAACAACTCTGgcgcgc                                                                                                               >  1:1892264/1‑148 (MQ=255)
                                        gAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATGGCGCACAAGACACCCCCCCACACCAAGAACAACTCTGGCGCGCCGGTa                                                                                                          >  1:2351328/1‑148 (MQ=255)
                                         aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTaa                                                                                                         >  2:1261212/1‑148 (MQ=255)
                                             acaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcgc                                                                                                     <  1:1084935/148‑1 (MQ=255)
                                                 cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCAAATACGGTCGCTTTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCAATCATAGAACAACTCTGGCGCGCCGGTAAGCGCGtat                                                                                                 <  1:885192/148‑1 (MQ=255)
                                                    gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATc                                                                                              >  2:1940933/1‑148 (MQ=255)
                                                    gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATc                                                                                              <  2:3727777/148‑1 (MQ=255)
                                                       aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTTGGGAACTAGACACCCCGCCACACATAGAACAACTCGGGCGCGCCGGTAAGCGCGGATATCtta                                                                                           >  1:1998925/1‑147 (MQ=255)
                                                       aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTTATTGCGAACTAGACACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGGATATAtta                                                                                           >  1:2367958/1‑147 (MQ=255)
                                                       aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTGCTCACTTTTGAAGGCGAACTAGACACCCAGCCACACATAGAACAACACTGGCCCCCCCCAAAGCGCGAATATCtaa                                                                                           >  2:1438474/1‑146 (MQ=255)
                                                        tGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTcc                                                                                          <  2:205756/148‑1 (MQ=255)
                                                        tGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTTATTGCGAACTAGAAACCCAGCCACAAATAGAACAACAATTGCGCGCCGGTAAGCGCCTATATAttac                                                                                          >  1:1539537/1‑146 (MQ=255)
                                                        tGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGATTGCGAACTAGACACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGTATATAttgc                                                                                          >  1:1539582/1‑146 (MQ=255)
                                                         gggTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATcatgca                                                                                         >  2:1786522/1‑143 (MQ=255)
                                                           gTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGt                                                                                       <  2:3200909/148‑1 (MQ=255)
                                                            tCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGtt                                                                                      <  2:2119231/148‑1 (MQ=255)
                                                               acaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATATTCCAGTTTTg                                                                                   >  2:864466/1‑148 (MQ=255)
                                                                 aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTACACTTTTGcc                                                                                 >  2:2415500/1‑147 (MQ=255)
                                                                   gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCCGCCACACCTAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCCGTTTTGCCTc                                                                               >  2:223738/1‑148 (MQ=255)
                                                                   gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTTCAta                                                                               >  2:2053784/1‑147 (MQ=255)
                                                                      gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCCGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                                            <  1:345347/148‑1 (MQ=255)
                                                                      gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                                            <  2:255587/148‑1 (MQ=255)
                                                                       cTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCGCTCTTTTGATGGCGAACTAGACAACCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTACAGTTTTGCCTCCCtt                                                                           >  2:1127127/1‑148 (MQ=255)
                                                                           ccGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCACTCATTTGATGGCGACCTAGACAACCCCCCACACAAAGAACAACTCTGGCGCGCCCGTAACCGCCGATATCTTACAGATTTTCAACCCtgaagt                                                                       >  2:3370729/1‑143 (MQ=255)
                                                                            cGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTg                                                                      <  1:20123/148‑1 (MQ=255)
                                                                            cGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTg                                                                      <  1:20429/148‑1 (MQ=255)
                                                                             ggATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGc                                                                     <  1:485725/148‑1 (MQ=255)
                                                                             ggATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGc                                                                     <  2:851155/148‑1 (MQ=255)
                                                                              gatgatGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTTTTGGCGAACTAGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGcc                                                                    >  2:1946759/1‑148 (MQ=255)
                                                                              gatgatGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCGGTTTTTTTTGCGAACTATAAACCCAGCCACACATAGAAAAACTCGGGGGCGCCGGTAAGCGCGGATATCTTCCAGGTTTGCCTCCCTGTAGTggc                                                                    >  1:2335730/1‑146 (MQ=255)
                                                                              gatgatGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGTTTGCGCTCTTTTTATGGCGAACTAGACACCCAGCCACACATAGAAAAACTCTGGCGCGCCGGGAAGCGCGGATATCTTCCAGGTTTTCATACCGTTAGTggc                                                                    >  2:978714/1‑146 (MQ=255)
                                                                               atgatgACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCg                                                                   >  2:3779911/1‑148 (MQ=255)
                                                                                 gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTGTTACTCTCTTTTGATGGCGCAAAATAAACCCAGCCACAAAAATAACAACTCAGGCGCGCCGGGAAGCGCGGATATCAAACACTTTTTCATCCCGTTAGTGGCGGt                                                                 >  1:1555969/1‑148 (MQ=255)
                                                                                          tCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTg                                                        >  2:3469816/1‑148 (MQ=255)
                                                                                          tCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTCCGCTCTTTTTATGGCGGACTAGACACCCCGCCACTCATATAACAAAACTGGCGCGCCGGGAAGCGCGTATATATTACAGTTTTTGATCCCGTTAGGGGCGGGAATTTGGTg                                                        >  2:1324933/1‑148 (MQ=255)
                                                                                             cccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTc                                                     <  2:1298062/148‑1 (MQ=255)
                                                                                             cccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTc                                                     <  2:2839283/148‑1 (MQ=255)
                                                                                             cccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCAGGTAATTTGGTGTTc                                                     <  2:927273/148‑1 (MQ=255)
                                                                                             cccGGTGGTGATGGCACATACGGTCGCTGTAGTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTc                                                     <  1:2925818/148‑1 (MQ=255)
                                                                                                 gtggtgATGGCACATACGGTCGCTGTGTTCCTCTCTTTTGATGGCGAACAAGACACCCAGCCACTCATAGAAAAACACGGGCGCGCCGGGAAGCGCGTATATCTTCCCGTTTTGCATCCCGTTAGTGGCGGGACTttgctgttgctgc                                                 >  2:102423/1‑147 (MQ=255)
                                                                                                  tggtgATGGCACATACGGTCGCGGTTTTCCTCTCATTTGATTGCGAACTAGACACCCAGCCACTCATAGAACAACACTGGCGCGCCGGTTAGCGCGGATATCTTCCAGTTTTGCAACCCTTTAGTGCCGGGAATTTGCTGTTACTTaa                                                >  2:2970607/1‑148 (MQ=255)
                                                                                                     tgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCCTAGAACAACACTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAACCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTa                                             >  2:3429881/1‑148 (MQ=255)
                                                                                                     tgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTa                                             <  2:713886/148‑1 (MQ=255)
                                                                                                        tGGCACATACGGTCGCTGTATTGCTCTCTTTTTATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCa                                          >  2:2405698/1‑148 (MQ=255)
                                                                                                          gCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATc                                        <  1:3068264/148‑1 (MQ=255)
                                                                                                               tACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCaa                                   <  1:1695084/148‑1 (MQ=255)
                                                                                                               tACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCaa                                   <  2:3067097/148‑1 (MQ=255)
                                                                                                                aCGGTCGCTGTATTCCTCTATTTTTATGTCGAAATATACAACCACCCACACATATAACAACACAGGCGCGCCCGTAAGCCCGTATATAATACAGTTTTGCATCCCTTTAGTGCCCGTAAATTTCTTTTACTGAATTAACATcacccaa                                  >  2:3732021/1‑142 (MQ=255)
                                                                                                                 cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAg                                 <  1:3779911/148‑1 (MQ=255)
                                                                                                                 cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAg                                 <  2:1623957/148‑1 (MQ=255)
                                                                                                                 cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAg                                 <  1:2795048/148‑1 (MQ=255)
                                                                                                                 cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAg                                 <  2:3625217/148‑1 (MQ=255)
                                                                                                                   gTCGCTGTATTCCTCTCTTTTTATGGCGAACTAGACACCCAGCCACTCATAGAACAAATCTGGCGCGCCGGTAAGCGCGTAAATCTTCCAGTTTTGCATACCTTTAGTGCCGGTAATTTGCTGTTACTGAATTACCATCCGCAAAGCg                               >  2:1722152/1‑148 (MQ=255)
                                                                                                                     cGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGGATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGaa                             <  2:3679429/148‑1 (MQ=255)
                                                                                                                      gCTGTAGTCCTCTCTGTGCATGGCGAACTAGACACCCAGCCACTCAGAGAACACCTCTGGCGCGCCGGTCAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAATCGAAc                            <  2:3570479/148‑1 (MQ=255)
                                                                                                                        tGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTg                          <  1:718125/148‑1 (MQ=255)
                                                                                                                          tATTCCTCTCTTTTGATGGCGAACTAGATACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                        >  2:647054/1‑148 (MQ=255)
                                                                                                                          tATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                        >  1:1119750/1‑148 (MQ=255)
                                                                                                                          tATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                        >  1:1015612/1‑148 (MQ=255)
                                                                                                                           aTTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtg                       <  1:3704935/148‑1 (MQ=255)
                                                                                                                           aTTCATCTCTTTTGATGGCGAACTAGACACCCAGCCACTCTTATAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtg                       <  1:3201850/148‑1 (MQ=255)
                                                                                                                            ttCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtga                      <  1:3429155/148‑1 (MQ=255)
                                                                                                                               ctctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatga                   <  1:426916/148‑1 (MQ=255)
                                                                                                                                tctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGaa                  <  2:2312692/148‑1 (MQ=255)
                                                                                                                                tctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGaa                  <  1:695339/148‑1 (MQ=255)
                                                                                                                                 ctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                 <  2:1892264/148‑1 (MQ=255)
                                                                                                                                   ctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAg               <  2:2505870/148‑1 (MQ=255)
                                                                                                                                    ttttGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAgg              >  2:210474/1‑148 (MQ=255)
                                                                                                                                    ttttGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAgg              >  2:210463/1‑148 (MQ=255)
                                                                                                                                      ttGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGtt            <  2:2367958/148‑1 (MQ=255)
                                                                                                                                        gATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGa          <  2:2904403/148‑1 (MQ=255)
                                                                                                                                        gATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTCATGAACAGGTTGa          <  2:1319376/148‑1 (MQ=255)
                                                                                                                                          tGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAg        >  1:1649354/1‑148 (MQ=255)
                                                                                                                                          tGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAg        <  1:3634025/148‑1 (MQ=255)
                                                                                                                                           ggCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGa       <  1:1324922/148‑1 (MQ=255)
                                                                                                                                             cGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATc     <  1:1266654/148‑1 (MQ=255)
                                                                                                                                             cGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATc     <  1:3515277/148‑1 (MQ=255)
                                                                                                                                               aaCTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa   <  1:112803/148‑1 (MQ=255)
                                                                                                                                               aaCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa   >  2:2956404/1‑148 (MQ=255)
                                                                                                                                                aCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAt  >  2:257168/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAT  >  CP009273/3049593‑3049884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ 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
TCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGC  >  CP009273/3049592‑3049888
                                                                                                                                                    |                                                                                                                                                    
TCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAG                                                                                                                                                     <  SRR19973069.1198133/150‑1 (MQ=60)
TCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAG                                                                                                                                                     <  SRR19973069.1200124/150‑1 (MQ=60)
    CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGTTGGCGAACTAGACAC                                                                                                                                                 >  SRR19973069.1931003/1‑150 (MQ=60)
    CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTGCTCTCTTTTGTTTGGGAACTAGACAC                                                                                                                                                 >  SRR19973069.1117856/1‑150 (MQ=60)
       AATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCA                                                                                                                                              <  SRR19973069.1496942/150‑1 (MQ=60)
           ATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCA                                                                                                                                          <  SRR19973069.2477667/150‑1 (MQ=60)
                    CGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACGCTGGTTTGATTGCGGACTAGACACCCAACCACTCATAGAA                                                                                                                                 >  SRR19973069.1823979/1‑150 (MQ=60)
                     GTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAAC                                                                                                                                <  SRR19973069.3430310/150‑1 (MQ=60)
                     GTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAAC                                                                                                                                <  SRR19973069.3429971/150‑1 (MQ=60)
                     GTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCGTTTTATTGCGAACAAGACACCCAGCCACTCATAGAAC                                                                                                                                >  SRR19973069.235948/1‑150 (MQ=60)
                               GTTACCGCCGGAACAACAGCAGGAATTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCG                                                                                                                      <  SRR19973069.3231573/150‑1 (MQ=60)
                               GTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCG                                                                                                                      <  SRR19973069.1745628/150‑1 (MQ=60)
                               GTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCG                                                                                                                      <  SRR19973069.3231572/150‑1 (MQ=60)
                                 TAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGCACTAGACACCCAGCCCCTCATAGAACAACACTGGCGCG                                                                                                                    >  SRR19973069.2018297/1‑150 (MQ=60)
                                  AACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATTGCGAACTAGACACCCAGCCCCACATAGAACAACACAGGCGCGC                                                                                                                   >  SRR19973069.3048646/1‑150 (MQ=60)
                                  AACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACAAGACACCCACCCACACAAAGAACAACTCTGGCGCGC                                                                                                                   >  SRR19973069.1905383/1‑150 (MQ=60)
                                   ACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCC                                                                                                                  <  SRR19973069.3696438/150‑1 (MQ=60)
                                    CGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCG                                                                                                                 <  SRR19973069.1565125/150‑1 (MQ=60)
                                       CGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATGGCGCACAAGACACCCCCCCACACCAAGAACAACTCTGGCGCGCCGGTA                                                                                                              >  SRR19973069.2367317/1‑150 (MQ=60)
                                              ACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGT                                                                                                       <  SRR19973069.1092511/150‑1 (MQ=60)
                                                AGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTGTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAAATCTGGCGCGCCGGTAAGCGCGGAT                                                                                                     >  SRR19973069.2888413/1‑150 (MQ=60)
                                                  CAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTTTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATAT                                                                                                   <  SRR19973069.1658892/150‑1 (MQ=60)
                                                  CAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCAAATACGGTCGCTTTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCAATCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATAT                                                                                                   <  SRR19973069.891312/150‑1 (MQ=60)
                                                     GAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTT                                                                                                <  SRR19973069.3463300/150‑1 (MQ=60)
                                                      AAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTTGGGAACTAGACACCCCGCCACACATAGAACAACTCGGGCGCGCCGGTAAGCGCGGATATCTTA                                                                                               >  SRR19973069.2012852/1‑150 (MQ=60)
                                                      AAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTTATTGCGAACTAGACACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGgatatatta                                                                                               >  SRR19973069.2384031/1‑141 (MQ=60)
                                                       AATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTTATTGCGAACTAGAAACCCAGCCACAAATAGAACAACAATTGCGCGCCGGTAAGCGCCTATATattac                                                                                              >  SRR19973069.1550087/1‑145 (MQ=60)
                                                       AATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGATTGCGAACTAGACACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGTATATattgc                                                                                              >  SRR19973069.1550132/1‑145 (MQ=60)
                                                               AACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTACTCTCTTTTGTTGGCGGACTAGACACCCACCCCCGCATAGAACAAATCTGGCGCGCCGGGAAGCGCGGATATCTTACAGTTTTTC                                                                                      >  SRR19973069.1533867/1‑150 (MQ=60)
                                                                 CAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCTCGTTTTTTGGGGAACTAGACACCCACCCACACCAAGAACAACACTGGCGCGCCGGGAAGCGCGTATATattaacctttttcat                                                                                    >  SRR19973069.2382152/1‑135 (MQ=60)
                                                                     CCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCT                                                                                <  SRR19973069.487283/150‑1 (MQ=60)
                                                                       GCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCCGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTT                                                                              <  SRR19973069.347798/150‑1 (MQ=60)
                                                                         TACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGGACTAGACACCCAGCCACCCATAGTACAAAACTGGCGCGCCGGTAAGCGCGGATATCTTCCAGTTTTGGATCCCGTTAG                                                                            >  SRR19973069.2101804/1‑150 (MQ=60)
                                                                         TACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAG                                                                            <  SRR19973069.3636060/150‑1 (MQ=60)
                                                                           CCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTG                                                                          <  SRR19973069.3721893/150‑1 (MQ=60)
                                                                             CGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCC                                                                        <  SRR19973069.20423/150‑1 (MQ=60)
                                                                             CGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCC                                                                        <  SRR19973069.20736/150‑1 (MQ=60)
                                                                             CGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCGGTTTTTTTTGCGAACTATAAACCCAGCCACACATAGAAAAACTCGGGGGCGCCGGTAAGCGCGGATATCTTCCAGGTTTGCCTCCCTGTAGTGGC                                                                        >  SRR19973069.2351620/1‑150 (MQ=50)
                                                                              GGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCG                                                                       <  SRR19973069.489236/150‑1 (MQ=60)
                                                                                ATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGT                                                                     <  SRR19973069.1291407/150‑1 (MQ=60)
                                                                                              CCCGGTGGTGATGGCACATACGGTCGCTGTAGTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCT                                                       <  SRR19973069.2945209/150‑1 (MQ=60)
                                                                                                 GGTGGTGATGGCACATACGGTCGCTGTATTCCGCGCTTTTTATTGGGGAAAAGACACCCACCCACACAAATAACAAAACTGGCGCGCCCGTAAGCGCGTATATATTACAGTTTTTCATACCGTTAGTGCCGGGAATTTGCTGTTACGGAA                                                    >  SRR19973069.3055949/1‑150 (MQ=17)
                                                                                                         TGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATC                                            <  SRR19973069.3137088/150‑1 (MQ=60)
                                                                                                           GCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCG                                          <  SRR19973069.3088539/150‑1 (MQ=60)
                                                                                                                TACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAG                                     <  SRR19973069.1706931/150‑1 (MQ=60)
                                                                                                                  CGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCG                                   <  SRR19973069.749850/150‑1 (MQ=60)
                                                                                                                  CGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCG                                   <  SRR19973069.3804453/150‑1 (MQ=60)
                                                                                                                  CGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCG                                   <  SRR19973069.2813625/150‑1 (MQ=60)
                                                                                                                       GCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTG                              >  SRR19973069.3658262/1‑150 (MQ=60)
                                                                                                                         TGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGT                            >  SRR19973069.1127559/1‑150 (MQ=60)
                                                                                                                         TGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGT                            >  SRR19973069.1022619/1‑150 (MQ=60)
                                                                                                                         TGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGT                            <  SRR19973069.723185/150‑1 (MQ=60)
                                                                                                                            ATTCCTCTCTTTTGATGGCGAACTAGCCACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT                         <  SRR19973069.3729314/150‑1 (MQ=60)
                                                                                                                            ATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT                         >  SRR19973069.637622/1‑150 (MQ=60)
                                                                                                                            ATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT                         <  SRR19973069.3729142/150‑1 (MQ=60)
                                                                                                                            ATTCATCTCTTTTGATGGCGAACTAGACACCCAGCCACTCTTATAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT                         <  SRR19973069.3222675/150‑1 (MQ=60)
                                                                                                                             TTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATG                        <  SRR19973069.3451449/150‑1 (MQ=60)
                                                                                                                                CTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAC                     <  SRR19973069.430047/150‑1 (MQ=60)
                                                                                                                                CTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAC                     <  SRR19973069.1954297/150‑1 (MQ=60)
                                                                                                                                 TCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACA                    <  SRR19973069.700223/150‑1 (MQ=60)
                                                                                                                                    CTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGT                 >  SRR19973069.2971591/1‑150 (MQ=60)
                                                                                                                                         GATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAG            >  SRR19973069.1660740/1‑150 (MQ=60)
                                                                                                                                           TGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGAT          <  SRR19973069.3657784/150‑1 (MQ=60)
                                                                                                                                           TGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGAT          <  SRR19973069.2421970/150‑1 (MQ=60)
                                                                                                                                            GGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATC         <  SRR19973069.1334030/150‑1 (MQ=60)
                                                                                                                                            GGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATC         <  SRR19973069.1334041/150‑1 (MQ=60)
                                                                                                                                              CGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA       <  SRR19973069.3537951/150‑1 (MQ=60)
                                                                                                                                              CGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA       <  SRR19973069.1275180/150‑1 (MQ=60)
                                                                                                                                                AACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATG     <  SRR19973069.113906/150‑1 (MQ=60)
                                                                                                                                                AACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATG     <  SRR19973069.114389/150‑1 (MQ=60)
                                                                                                                                                 ACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGA    >  SRR19973069.2736947/1‑150 (MQ=60)
                                                                                                                                                   TAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGC  >  SRR19973069.622287/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                                    
TCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGC  >  CP009273/3049592‑3049888

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