Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 32,912 A→C 5.2% E699A (GAA→GCA) ‡ carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091332,9120AC 5.2% 466.7 / 8.3 151E699A (GAA→GCA) ‡carBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (72/69);  new base C (4/5);  total (76/75)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

CTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGC  >  NC_000913/32764‑33061
                                                                                                                                                    |                                                                                                                                                     
cTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGaa                                                                                                                                                      >  1:795310/1‑150 (MQ=255)
cTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGAGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGaa                                                                                                                                                      >  2:803041/1‑150 (MQ=255)
 tGGCGCGCGCGCTGGAAGCTGCTGGCGTACCTGTTATAGGCACCATCCCGGATGATATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGaaa                                                                                                                                                     <  1:156948/150‑1 (MQ=255)
      ggcgcgcGGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCg                                                                                                                                                <  1:143195/149‑1 (MQ=255)
       gcgcgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCgg                                                                                                                                               <  2:1040703/150‑1 (MQ=255)
       gcgcgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCgg                                                                                                                                               <  1:541355/150‑1 (MQ=255)
       gcgcgcTGGAAGCTGCGGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCgg                                                                                                                                               <  1:708843/150‑1 (MQ=255)
          cgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTag                                                                                                                                            >  2:869375/1‑150 (MQ=255)
           gcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaac                                                                                                                                           <  2:1108271/150‑3 (MQ=255)
            gtGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                          <  2:543693/149‑1 (MQ=255)
            cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                          >  1:1022400/1‑150 (MQ=255)
            cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                          >  2:529355/1‑150 (MQ=255)
            cTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTagag                                                                                                                                          >  1:896321/1‑150 (MQ=255)
              ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGaa                                                                                                                                        <  2:1683276/150‑1 (MQ=255)
              ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGaa                                                                                                                                        >  1:1278199/1‑150 (MQ=255)
              ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggt                                                                                                                                        <  1:1489593/150‑6 (MQ=255)
              ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggt                                                                                                                                        <  1:819562/150‑6 (MQ=255)
                aaGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGGTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggtgg                                                                                                                                      <  1:485407/150‑8 (MQ=255)
                   cggctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCCGAAGACCGGGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGTAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGa                                                                                                                                   <  2:64106/148‑1 (MQ=255)
                     gctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGCATCAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGaaa                                                                                                                                 >  1:417994/1‑150 (MQ=255)
                      ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                                >  1:447228/1‑150 (MQ=255)
                      ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                                <  1:1746223/150‑1 (MQ=255)
                       tgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAaga                                                                                                                               <  2:1279824/150‑1 (MQ=255)
                        ggCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCg                                                                                                                                    <  2:122434/144‑1 (MQ=255)
                          cGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGAGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGCGCGTCGGAACCTGAAACAACCGGCGCACGCCCCCGTGACCGCTATTGAAATGGCGGTAGAGACGGCGATAGAGAt                                                                                                                            >  1:298106/1‑150 (MQ=255)
                            tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                          <  1:1306051/150‑1 (MQ=255)
                            tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                          <  1:1529183/150‑1 (MQ=255)
                               cGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCt                                                                                                                       <  2:903822/150‑1 (MQ=255)
                                ggTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTTTTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTa                                                                                                                      >  1:1623689/1‑150 (MQ=255)
                                    aTCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGGACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACCACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCg                                                                                                                  >  1:1158837/1‑150 (MQ=255)
                                     tCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGAGAAAGAGATTGGCTACCCGc                                                                                                                 >  2:673322/1‑150 (MQ=255)
                                       ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                               >  2:1205512/1‑150 (MQ=255)
                                       ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                               >  2:1573787/1‑150 (MQ=255)
                                        gCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtgg                                                                                                              <  1:198020/150‑1 (MQ=255)
                                           ccAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggtgg                                                                                                           >  1:1077578/1‑150 (MQ=255)
                                             aGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTa                                                                                                         <  1:1413232/150‑1 (MQ=255)
                                                   gATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccg                                                                                                   <  1:1580043/150‑1 (MQ=255)
                                                      gCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGATTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCt                                                                                                >  2:33908/1‑150 (MQ=255)
                                                       cTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCtt                                                                                               >  1:1064049/1‑150 (MQ=255)
                                                       cTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCtt                                                                                               >  2:1192382/1‑150 (MQ=255)
                                                        tATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGACACACCCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTa                                                                                              >  2:513034/1‑150 (MQ=255)
                                                         aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                             >  1:612710/1‑150 (MQ=255)
                                                         aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                             <  2:1757888/150‑1 (MQ=255)
                                                           cGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGt                                                                                           >  2:1646173/1‑150 (MQ=255)
                                                                gTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcg                                                                                      >  2:1641738/1‑150 (MQ=255)
                                                                gTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcg                                                                                      >  2:1535362/1‑150 (MQ=255)
                                                                gTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCAGGTGGTACGTCCGTCTTACGTTCTcg                                                                                      >  1:985880/1‑150 (MQ=255)
                                                                 tGCAGAAGACCGTGAACGCTTCCTGCATGCGGTTGCGCGTCTGAAACTGAAACAACCGGCGAACGCCCCCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGATACCCGCTGGGGTTACGTCCGTCTTACGTTCTcgg                                                                                     >  2:1127326/1‑150 (MQ=255)
                                                                  gCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggc                                                                                    >  1:1732004/1‑150 (MQ=255)
                                                                   cAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcg                                                                                   <  2:339244/150‑1 (MQ=255)
                                                                    agaagaCCGAGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcgg                                                                                  <  1:1809317/150‑1 (MQ=255)
                                                                       cgaCCGGGAACGCTTCCAGCATGCGGTGGATCGTCTGAACCTGAAACACCCGGCGACCGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGTTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                               <  2:1562168/149‑1 (MQ=255)
                                                                       agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                               <  2:610097/150‑1 (MQ=255)
                                                                       agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                               <  2:1300713/150‑1 (MQ=255)
                                                                          ccGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGc                                                                            >  2:1295999/1‑150 (MQ=255)
                                                                              gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                        <  2:896321/150‑1 (MQ=255)
                                                                              gAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATg                                                                        <  2:1706309/150‑1 (MQ=255)
                                                                               aaCGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATgg                                                                       <  2:234245/150‑1 (MQ=255)
                                                                                      ccAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGATATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCTGTCGGGCGATGGAAAACGt                                                                >  2:1388567/1‑150 (MQ=255)
                                                                                       cAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTc                                                               <  2:746099/150‑1 (MQ=255)
                                                                                        tgCATGCGGTTGAGCGGCTGAAACTGAAACAACGGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTTGAGACGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCt                                                              <  2:301565/149‑1 (MQ=255)
                                                                                         gCATGCGGTTGAGAGTCTGAAACTGAAACAACCGGCGAACGCCACCGTAACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTAGCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTa                                                             <  1:617191/150‑1 (MQ=255)
                                                                                          cATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTAt                                                            >  2:413668/1‑150 (MQ=255)
                                                                                             gCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGAc                                                         >  2:504012/1‑150 (MQ=255)
                                                                                               ggTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGa                                                       >  2:203112/1‑150 (MQ=255)
                                                                                                  tGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGc                                                    <  1:1409713/150‑1 (MQ=255)
                                                                                                          tgaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGc                                            >  1:1304197/1‑150 (MQ=255)
                                                                                                           gaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCg                                           >  2:1049416/1‑150 (MQ=255)
                                                                                                             atctgaaacACCCGGCGAAGGCCACCGTTACCGCTAGTGTAATGGCGGTAGTGAAGGCGTAAGTGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAATTCGTCTATGACGATGCTGACCTGCGTc                                         <  1:1247066/148‑1 (MQ=255)
                                                                                                             aactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTc                                         >  1:725272/1‑150 (MQ=255)
                                                                                                              actgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCg                                        >  1:1181829/1‑150 (MQ=255)
                                                                                                               ctgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGc                                       >  2:702308/1‑150 (MQ=255)
                                                                                                                tgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCt                                      <  1:513625/150‑1 (MQ=255)
                                                                                                                tgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCt                                      <  1:1099857/150‑1 (MQ=255)
                                                                                                                 gaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTa                                     >  2:124298/1‑150 (MQ=255)
                                                                                                                  aaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTAc                                    >  2:1473423/1‑150 (MQ=255)
                                                                                                                  aaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACACGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATAGTCTATGTCGAAGCTGACCTGCGTCGCTAc                                    >  2:746637/1‑150 (MQ=255)
                                                                                                                     caacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTc                                 >  2:1461282/1‑150 (MQ=255)
                                                                                                                      aacCGGCGTACGCCACCGTTACCGCTATTGAATTGGCGGTAGAGACTGCGAAAGAGATTGGCTCCCCGCTGGTGGTACGTCCGTCTTACTGTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                                <  2:277171/150‑1 (MQ=255)
                                                                                                                      aacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                                <  2:464975/150‑1 (MQ=255)
                                                                                                                       acCGGCGAACGCCACTGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCa                               >  2:1728177/1‑150 (MQ=255)
                                                                                                                        ccGGCGAACGCCACCGTTACAGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAAGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGAGTCGCTACTTCCAg                              >  2:413063/1‑150 (MQ=255)
                                                                                                                         cGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGa                             >  1:243599/1‑150 (MQ=255)
                                                                                                                            cGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcgg                          >  1:521232/1‑150 (MQ=255)
                                                                                                                             gAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcggc                         <  2:1722305/150‑1 (MQ=255)
                                                                                                                                     tctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg                 >  2:146585/5‑150 (MQ=255)
                                                                                                                                     ccGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGATGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg                 <  1:1834294/150‑1 (MQ=255)
                                                                                                                                     ccGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg                 <  1:610686/150‑1 (MQ=255)
                                                                                                                                      ctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgt                >  2:976802/4‑150 (MQ=255)
                                                                                                                                       gTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtg               <  1:1049994/150‑1 (MQ=255)
                                                                                                                                       gTTACCGCTAATGAACTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtg               <  1:862146/150‑1 (MQ=255)
                                                                                                                                           ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa           <  1:560097/150‑1 (MQ=255)
                                                                                                                                           ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa           <  1:1620288/150‑1 (MQ=255)
                                                                                                                                           ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa           >  1:1028307/1‑150 (MQ=255)
                                                                                                                                             gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc         >  1:1126578/1‑150 (MQ=255)
                                                                                                                                             gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGGGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc         <  1:291790/150‑1 (MQ=255)
                                                                                                                                             gCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGATGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCTATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAAc         >  1:1132414/1‑150 (MQ=255)
                                                                                                                                                aTTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGCCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGAt      >  1:1638932/1‑150 (MQ=255)
                                                                                                                                                    aaaTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCGCGGCTGTCGGGCGATGGAAATCGTCTATGACGCAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGCTgcgc  >  2:488800/1‑150 (MQ=255)
                                                                                                                                                    aaaTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGGCTGGCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcgc  <  1:1205512/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGCGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGC  >  NC_000913/32764‑33061

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