Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,894,478 A→T 7.7% intergenic (+46/+74) yoaC → / ← yoaH DUF1889 domain‑containing protein YoaC/UFP0181 family protein YoaH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,894,4780AT 7.7% 366.6 / 15.2 117intergenic (+46/+74)yoaC/yoaHDUF1889 domain‑containing protein YoaC/UFP0181 family protein YoaH
Reads supporting (aligned to +/- strand):  ref base A (52/56);  new base T (4/5);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 123 total aligned reads displayed.

TGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCTCTTCTGCCACCAG  >  NC_000913/1894329‑1894626
                                                                                                                                                     |                                                                                                                                                    
tGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATcaca                                                                                                                                                      <  1:241695/150‑1 (MQ=255)
 ggTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACat                                                                                                                                                     >  1:40077/1‑150 (MQ=255)
         gggCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGTCGGTTTTAATCACATATaaaatg                                                                                                                                             <  2:1962702/150‑3 (MQ=255)
         gggAAAAAAAAATGGAGACATGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAg                                                                                                                                             <  1:868130/150‑1 (MQ=255)
          ggCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGTAATAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGTCGCTTTTTATCATATATAAAAAGc                                                                                                                                            <  2:1256731/150‑1 (MQ=255)
          ggCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGc                                                                                                                                            <  2:1591931/150‑1 (MQ=255)
          ggCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGc                                                                                                                                            >  2:1168553/1‑150 (MQ=255)
           gCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACAtata                                                                                                                                                  <  1:1676445/143‑1 (MQ=255)
           gCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACAtata                                                                                                                                                  >  2:1676445/1‑143 (MQ=255)
            cAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGt                                                                                                                                          >  2:396944/1‑150 (MQ=255)
            cAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGt                                                                                                                                          >  2:1912718/1‑150 (MQ=255)
              aaaaaaaTGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTTATATGTGATTAAAAGCGTcc                                                                                                                                        <  2:1184147/150‑1 (MQ=255)
               aaaaaaTGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTTATATGTGATTAAAAGCGTCCt                                                                                                                                       <  2:32161/150‑1 (MQ=255)
               aaaaaaTGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTTATATGTGATTAAAAGCGTCCt                                                                                                                                       <  2:1149672/150‑1 (MQ=255)
               aaaaaaTGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCt                                                                                                                                       >  2:1651803/1‑150 (MQ=255)
                       gagaCAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTca                                                                                                                               <  1:396944/150‑1 (MQ=255)
                        agaCAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTcac                                                                                                                              <  1:1854945/150‑1 (MQ=255)
                        agaCAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTcac                                                                                                                              <  2:336005/150‑1 (MQ=255)
                        agaCAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTcac                                                                                                                              <  1:987366/150‑1 (MQ=255)
                         taCAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTcaca                                                                                                                             <  2:1248964/149‑1 (MQ=255)
                            aGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATc                                                                                                                          <  2:288494/150‑1 (MQ=255)
                             ggTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCa                                                                                                                         >  2:1253269/1‑150 (MQ=255)
                             ggTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCa                                                                                                                         >  1:126756/1‑150 (MQ=255)
                             ggTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCa                                                                                                                         >  1:679189/1‑150 (MQ=255)
                                  acattCTTTGATTATAAATCCAGAATACTTTTCAATATATATTCTTGAAGAACTAAAATGACTTGTCTTATTTAAATTTATATCATCTTAAATGGTTCAGTACGCTTTTAATCAAATATAAAAAGCGTCCTTTTTCTTACAATCAAATGa                                                                                                                    <  2:893363/147‑1 (MQ=255)
                                  aCGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGa                                                                                                                    <  1:1004121/150‑1 (MQ=255)
                                  aCGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGa                                                                                                                    <  2:704895/150‑1 (MQ=255)
                                       cTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGatcatc                                                                                                               >  2:1345964/1‑150 (MQ=255)
                                        tgtgATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCg                                                                                                              >  1:1135552/1‑150 (MQ=255)
                                        tgtgATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCg                                                                                                              >  1:1128504/1‑150 (MQ=255)
                                         gtgATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACat                                                                                                                                                     >  1:1734846/1‑110 (MQ=255)
                                          tgATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGc                                                                                                            <  1:1596453/150‑1 (MQ=255)
                                          tgATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGc                                                                                                            >  1:923323/1‑150 (MQ=255)
                                           gATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTTGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACTCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGTCt                                                                                                           <  2:747048/150‑1 (MQ=255)
                                             ttAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTc                                                                                                         <  1:1168553/150‑1 (MQ=255)
                                              tAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCt                                                                                                        <  1:1383391/150‑1 (MQ=255)
                                                    tCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACg                                                                                                  >  2:1822827/1‑150 (MQ=255)
                                                        gAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAAt                                                                                              >  1:236111/1‑150 (MQ=255)
                                                        gAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAAt                                                                                              <  1:1521076/150‑1 (MQ=255)
                                                           tACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCg                                                                                                                                           <  2:1537490/102‑1 (MQ=255)
                                                           tACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGa                                                                                           <  1:1912718/150‑1 (MQ=255)
                                                            aCTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGAt                                                                                          <  2:933208/150‑1 (MQ=255)
                                                                 tCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGa                                                                                     >  1:1321201/1‑150 (MQ=255)
                                                                     cATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAg                                                                                 >  1:408314/1‑150 (MQ=255)
                                                                         tatGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAgtgtg                                                                             >  2:919083/1‑150 (MQ=255)
                                                                                    aaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACt                                                                        >  1:1711129/1‑144 (MQ=255)
                                                                                    aaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACt                                                                        <  2:1711129/144‑1 (MQ=255)
                                                                                    aaCTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCtt                                                                  <  2:408314/150‑1 (MQ=255)
                                                                                    aaCTAAAAGCAATTGTGTGTGTTAAATTTATATCATCATAAATGGGTCATGACGGTTTTAATCACATATAAAAATCGTCCTTTTTCCCACAATCAAATGATCATCGGCTTCTATAACGAAATCGTTACGGATTAGTGTGTTACTCGTCtt                                                                  <  2:1005666/150‑1 (MQ=255)
                                                                                      cTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCa                                                                >  2:1404410/1‑150 (MQ=255)
                                                                                       tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCAt                                                               <  1:678840/150‑1 (MQ=255)
                                                                                       tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCAt                                                               >  2:1680974/1‑150 (MQ=255)
                                                                                       tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCAt                                                               >  2:1258429/1‑150 (MQ=255)
                                                                                       tAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCAt                                                               >  2:1414807/1‑150 (MQ=255)
                                                                                       tAAAAGCACTGGTATGAGTTAAATTTATATCAGAATAAATTTGTCAGGACGCTTTTAATAACATATAAAAAGCGTCCTTTTTATCACAATAAAATGATCATCGGATTATATAACGAAATCGATACTGATTATTGTGTTAATCGTCTTCAt                                                               >  2:1813719/1‑150 (MQ=255)
                                                                                         taaGCACTTGTCTTATTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATTATCATCGGGTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATcc                                                             <  2:1508003/149‑1 (MQ=255)
                                                                                         aaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATcc                                                             <  2:40077/150‑1 (MQ=255)
                                                                                         aaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATcc                                                             <  2:1128504/150‑1 (MQ=255)
                                                                                         aaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATcc                                                             >  2:615256/1‑150 (MQ=255)
                                                                                         aaaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATcc                                                             <  1:1822827/150‑1 (MQ=255)
                                                                                          aaGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTAATCGTATTCATcca                                                            >  1:1508003/1‑149 (MQ=255)
                                                                                           aGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTc                                                           <  1:1253269/150‑1 (MQ=255)
                                                                                            gCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCa                                                          <  1:1345964/150‑1 (MQ=255)
                                                                                            gCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCa                                                          <  2:1934353/150‑1 (MQ=255)
                                                                                             cACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaa                                                         >  2:34903/1‑150 (MQ=255)
                                                                                               cTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCaaaa                                                       >  1:1386269/1‑150 (MQ=255)
                                                                                                tGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAAc                                                      <  1:1404410/150‑1 (MQ=255)
                                                                                                   tCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAAcgcg                                                   >  2:283419/1‑150 (MQ=255)
                                                                                                    cTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACgcgc                                                  >  2:1573603/1‑150 (MQ=255)
                                                                                                    cTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACgcgc                                                  >  2:262812/1‑150 (MQ=255)
                                                                                                     tGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGcc                                                 <  2:923323/150‑1 (MQ=255)
                                                                                                     tGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGcc                                                 <  1:895831/150‑1 (MQ=255)
                                                                                                           aaaTTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGAtt                                           >  2:1317388/1‑150 (MQ=255)
                                                                                                            aaTTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCa                                                          <  2:1317551/134‑1 (MQ=255)
                                                                                                            aaTTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATtc                                          <  2:1321201/150‑1 (MQ=255)
                                                                                                            aaTTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATtc                                          <  1:283419/150‑1 (MQ=255)
                                                                                                             aTTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATtcg                                         >  2:41565/1‑150 (MQ=255)
                                                                                                              tttATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATtcgc                                        <  1:1258429/150‑1 (MQ=255)
                                                                                                                   atCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGcc                                   <  1:1651803/150‑1 (MQ=255)
                                                                                                                    tCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCg                                  >  1:608681/1‑150 (MQ=255)
                                                                                                                    tCAGCATAAATGGGTCAGGACGCTTTTAATCACATACAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCg                                  >  2:1170639/1‑150 (MQ=255)
                                                                                                                     cAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCgg                                 <  1:1317388/150‑1 (MQ=255)
                                                                                                                      aGCATAAATGGGTCAGGACGCTTTTAATAACATATAAAAAGCGTCCTTTTTCTCACAATCAAAACATCATCGGCTTCTATAAAGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCCCCACGATTCGCTCCCAGGa                                >  2:1523616/1‑150 (MQ=255)
                                                                                                                            aaTGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATACGATCGCCGGAGTGGtt                          >  1:350793/1‑150 (MQ=255)
                                                                                                                              tGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTgg                        <  2:850624/150‑1 (MQ=255)
                                                                                                                               gggTCAGGACTCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGGGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGc                       <  1:1523616/150‑1 (MQ=255)
                                                                                                                               gggTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGc                       <  1:1414807/150‑1 (MQ=255)
                                                                                                                                 gTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCAc                     >  1:762119/1‑150 (MQ=255)
                                                                                                                                   aaGGACGCTTTTTATATGTGATTAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGc                   >  2:1148812/2‑150 (MQ=255)
                                                                                                                                   aaGGACGCTTTTTATATGTGATTAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGc                   >  2:60899/2‑150 (MQ=255)
                                                                                                                                    aGGACGCTTTTTATATGTGATTAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCa                  >  2:1121923/1‑150 (MQ=255)
                                                                                                                                     ggACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAg                 >  2:1753584/1‑150 (MQ=255)
                                                                                                                                     ggACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAg                 >  1:232494/1‑150 (MQ=255)
                                                                                                                                      gACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGc                <  1:919083/150‑1 (MQ=255)
                                                                                                                                           ttttAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCCGCtcttg           >  2:1400088/1‑149 (MQ=255)
                                                                                                                                            tttAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCtcttct          <  1:752839/150‑1 (MQ=255)
                                                                                                                                              tAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCTCTTCTGc        <  1:1573603/150‑1 (MQ=255)
                                                                                                                                                aTCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCTCTTCTGcca      >  2:1375752/1‑150 (MQ=255)
                                                                                                                                                   acaTATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCTCTTCTGccacca   <  2:679189/150‑1 (MQ=255)
                                                                                                                                                    caTATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGACTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCTCTTCTGccaccac  >  1:1049842/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGTTGGATGGGCAAAAAAAATGGAGACAGGTGAACGTTCTGTGATTAAAAATCCTGAATACTTTTCAACATATATGCAGGAAGAACTAAAAGCACTGGTCTGAGTTAAATTTATATCAGCATAAATGGGTCAGGACGCTTTTAATCACATATAAAAAGCGTCCTTTTTCTCACAATCAAATGATCATCGGCTTCTATAACGAAATCGATACGGATTAGTGTGTTACTCGTCTTCATCCTCAAAACGCGCCACGATTCGCTCGCCGGAGTGGTTGGCACGCAGCTCTTCTGCCACCAG  >  NC_000913/1894329‑1894626

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