Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,970 A→G 5.7% intergenic (‑128/+111) tcdA ← / ← mltA tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,945,9700AG 5.7% 314.9 / 7.2 105intergenic (‑128/+111)tcdA/mltAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base A (54/45);  new base G (3/3);  total (57/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 109 total aligned reads displayed.

CATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGG  >  NC_000913/2945821‑2946115
                                                                                                                                                     |                                                                                                                                                 
cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACAAAACTCCTCCAATATTTTGGAATTGCATGATTTGTAGGCCGGATa                                                                                                                                                   >  1:1530237/1‑150 (MQ=255)
cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCATAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGtata                                                                                                                                                   >  2:1316553/1‑150 (MQ=255)
cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCATCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATACCTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATa                                                                                                                                                   >  1:1536781/1‑150 (MQ=255)
cATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATACCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATTACTCCTGCAATAGTTTGGATTTTCATGATGATTAGGCCGGATa                                                                                                                                                   >  2:877277/1‑150 (MQ=38)
    aCTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGc                                                                                                                                               <  1:438153/150‑1 (MQ=255)
       aaTTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGtg                                                                                                                                            >  2:80452/1‑149 (MQ=255)
       aaTTACCACAGACATAGCACCCCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAACCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGtt                                                                                                                                            >  2:834149/1‑150 (MQ=255)
         ttACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                          <  2:610931/150‑1 (MQ=255)
          tACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                         >  2:1192243/1‑150 (MQ=255)
           aCCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAcg                                                                                                                                        >  1:335152/1‑150 (MQ=255)
               caGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                    >  1:1682017/1‑150 (MQ=255)
                aGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                                                   <  2:787704/150‑1 (MQ=255)
                  aCATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGTATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                 >  2:1413281/1‑150 (MQ=255)
                   cATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                <  1:538007/150‑1 (MQ=37)
                    aTAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAAGGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                               <  2:25165/150‑1 (MQ=255)
                      aGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                             <  2:918481/150‑1 (MQ=255)
                       gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                            <  1:1510640/150‑2 (MQ=37)
                         aCCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGccta                                                                                                                          <  1:80452/150‑4 (MQ=37)
                         aCCTCAGAATTAAAAAGGGCGAGGATATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                          <  1:36942/150‑1 (MQ=255)
                          ccTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                                                                                         >  1:191011/1‑150 (MQ=255)
                           cTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                                                                                        >  2:597104/1‑150 (MQ=255)
                           cTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                                                                                                                        <  2:223965/150‑1 (MQ=255)
                            tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                       >  2:1093240/1‑150 (MQ=255)
                            tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                       <  1:903296/150‑1 (MQ=255)
                            tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                       <  2:509103/150‑1 (MQ=255)
                            tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                       <  1:661728/150‑1 (MQ=255)
                            tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                       >  1:457732/1‑150 (MQ=255)
                             cataaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca                                                                                                                      <  2:834587/147‑1 (MQ=255)
                               gAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                    >  1:460985/1‑150 (MQ=255)
                               gAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                    >  1:499079/1‑150 (MQ=255)
                                aaTTACAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATGAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                   <  1:1037313/150‑1 (MQ=255)
                                aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                   <  1:1334535/150‑1 (MQ=255)
                                    aaaaaGGGCGAGGTTATCCCTCACCCTACACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAAGATTTTGGATTTGCATGATTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                               >  1:89867/1‑150 (MQ=255)
                                    aaaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCa                                                                                                               >  2:1573992/1‑150 (MQ=255)
                                        aGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTa                                                                                                           >  1:420843/1‑150 (MQ=255)
                                         gggCGAGGTTATCCCACACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCTGCACAACGCGCACTTAg                                                                                                          >  2:319868/1‑150 (MQ=255)
                                           ccGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGAGCTACGCTTATCAGGCCTACATAACTCCTGCAAGATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGAACTTAGTc                                                                                                        >  1:135479/2‑150 (MQ=255)
                                               ggTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaata                                                                                                    >  2:505953/1‑150 (MQ=255)
                                                gTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCagtaa                                                                                                   >  2:283451/1‑150 (MQ=255)
                                                gTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCATACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaataa                                                                                                   <  2:1188993/150‑1 (MQ=255)
                                                   aTCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCa                                                                                                <  1:688770/150‑1 (MQ=255)
                                                      ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                             >  2:1458930/1‑150 (MQ=255)
                                                         cACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgt                                                                                          >  1:114978/1‑150 (MQ=255)
                                                         cACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgt                                                                                          <  1:1093240/150‑1 (MQ=255)
                                                         cACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgt                                                                                          >  1:508077/1‑150 (MQ=255)
                                                          acgcTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACATCGCGCACTTAGTCAATAATCAGGAtgtg                                                                                         <  1:752090/147‑1 (MQ=255)
                                                          aCCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATATTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtctg                                                                                         >  2:698660/1‑150 (MQ=255)
                                                           cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                        >  2:524576/1‑150 (MQ=255)
                                                           cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                        <  2:554642/150‑1 (MQ=255)
                                                           cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                        <  2:316009/150‑1 (MQ=255)
                                                             cTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAgg                                                                                      >  1:882110/1‑150 (MQ=255)
                                                             cTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAgg                                                                                      <  1:2617/150‑1 (MQ=255)
                                                               aaaCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGtt                                                                                    >  2:1473153/1‑150 (MQ=255)
                                                                  cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                                 <  2:418687/150‑1 (MQ=255)
                                                                  cGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcc                                                                                 >  2:911177/1‑150 (MQ=255)
                                                                   gCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAccc                                                                                >  2:1071766/1‑150 (MQ=255)
                                                                    cATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcccc                                                                               >  1:740629/1‑150 (MQ=255)
                                                                      tAAATGCCTCATCCGCTACGCTTATCAGGCCTACAAAACTCCTGCAATATTTTCGATTTGCATGATTTGTAGGCCTGATACGGCGTTCACGCCGCATCCGCGATGAACAACGCGCACTTAGTCATTAATCAGGATGTGAGGTTACGCCTc                                                                             >  2:872267/1‑150 (MQ=25)
                                                                         aTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCCCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAcc                                                                          >  1:1064965/1‑150 (MQ=255)
                                                                          tGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAccc                                                                         <  2:1135101/150‑1 (MQ=255)
                                                                            ccTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTg                                                                       <  1:1586369/150‑1 (MQ=255)
                                                                              tGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTc                                                                     <  2:1452457/150‑1 (MQ=255)
                                                                                       aCGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTa                                                            <  2:882110/150‑1 (MQ=255)
                                                                                        cGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaa                                                           <  1:1473153/150‑1 (MQ=255)
                                                                                          cTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaaa                                                         >  2:4070/1‑150 (MQ=255)
                                                                                            tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                       >  1:1552934/1‑150 (MQ=255)
                                                                                               cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                    >  1:818519/1‑150 (MQ=255)
                                                                                               cAGGCCTACATAACTCCAGCAAAATTTTGGATTTGCATGATTTGTCGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                    <  2:1401148/150‑1 (MQ=255)
                                                                                                tggCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                   <  2:335152/149‑1 (MQ=255)
                                                                                                aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                   <  1:1573992/150‑1 (MQ=255)
                                                                                                 ggCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                                  >  1:160508/1‑150 (MQ=255)
                                                                                                 ggCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCAt                                                  >  2:559230/1‑150 (MQ=255)
                                                                                                     tACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGGACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTATAAACGGCATACAg                                              >  2:257714/1‑150 (MQ=255)
                                                                                                      aCATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAATAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGa                                             >  2:886619/1‑150 (MQ=255)
                                                                                                        aTAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAACCCTCACCCTGGCATATACGAAAAAACTGCATACAGaaa                                           >  1:186184/1‑149 (MQ=37)
                                                                                                           aCTCCTGCAATTTTTTGGATTATCATGATTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGAATGAACAACGCGCACTTAGTAAATAATCAGGATGTGAGGTTACCCCTCACCCTGTAATATCCATAAAAACGGCATACAGAATATc                                        <  1:877277/150‑1 (MQ=255)
                                                                                                              ccTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATcaca                                     <  2:191011/150‑1 (MQ=255)
                                                                                                               cTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCAGCATCCTGCATGAACAAAGCGAACTTAGTCAATAATCAGGATGTGAGGTTACACCTCACACTGTCATATCCGTAAAAACGGCATACATAATATcacaa                                    >  1:280951/1‑149 (MQ=37)
                                                                                                                tGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATc                                   <  2:114978/150‑1 (MQ=255)
                                                                                                                 gCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATcc                                  >  1:1299614/1‑149 (MQ=255)
                                                                                                                 gCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCa                                  >  1:26557/1‑150 (MQ=255)
                                                                                                                 gCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCa                                  >  1:1292746/1‑150 (MQ=255)
                                                                                                                       ttttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCt                            >  1:595480/1‑150 (MQ=255)
                                                                                                                          tGGATTTGCATGATTTGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACTCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACTGCATACAGAATATCACATCAGCCGCTaaa                         <  1:1316553/150‑1 (MQ=255)
                                                                                                                           ggATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg                        <  2:460985/150‑1 (MQ=255)
                                                                                                                            gATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGa                       <  1:834149/150‑1 (MQ=255)
                                                                                                                             aTTTGTATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGAc                      <  1:1413281/150‑1 (MQ=255)
                                                                                                                                   aTGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTAcc                <  2:457732/150‑1 (MQ=255)
                                                                                                                                     gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTg              >  2:1479565/1‑150 (MQ=255)
                                                                                                                                     gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTg              <  1:44866/150‑1 (MQ=255)
                                                                                                                                         tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACACTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGcc          >  2:235516/1‑150 (MQ=255)
                                                                                                                                          gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGccc         >  2:1449434/1‑150 (MQ=255)
                                                                                                                                           catgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACTGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCg        >  1:21895/4‑150 (MQ=38)
                                                                                                                                             tgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCggg      >  1:1442873/2‑150 (MQ=39)
                                                                                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCggg      >  1:1075523/1‑150 (MQ=255)
                                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTCCCCCTCACCCTGTCATATCCGTAAAAACGGCAGACAGAATATCACCTCAGCCGCTAAAGACGTTACCTGCGCCCgggg     >  1:700725/1‑150 (MQ=255)
                                                                                                                                               ccGGATAATGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc    <  2:1530237/150‑1 (MQ=255)
                                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGTATGTGAGGGTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc    <  2:1536781/150‑1 (MQ=255)
                                                                                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc    >  2:787270/1‑150 (MQ=255)
                                                                                                                                                 ggATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCgg  <  2:499079/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CATCACTAATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGG  >  NC_000913/2945821‑2946115

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