Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,672,484 A→T 9.6% L175F (TTA→TTT ydgD → putative serine protease YdgD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,4840AT 9.6% 239.9 / 12.3 83L175F (TTA→TTTydgDputative serine protease YdgD
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base T (3/5);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGG  >  NC_000913/1672337‑1672630
                                                                                                                                                   |                                                                                                                                                  
gCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTaaa                                                                                                                                                  >  1:48816/1‑150 (MQ=255)
  ggTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAg                                                                                                                                                <  1:398719/150‑1 (MQ=255)
      aaaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGc                                                                                                                                            >  1:915262/1‑150 (MQ=255)
      aaaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGc                                                                                                                                            <  2:1066560/150‑1 (MQ=255)
      aaaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGc                                                                                                                                            <  1:1277268/150‑1 (MQ=255)
       aaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcg                                                                                                                                           >  2:1138448/1‑149 (MQ=255)
        aaGCAGATGGGGATGCTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcgg                                                                                                                                          <  1:895799/150‑3 (MQ=255)
         aGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAgg                                                                                                                                         >  2:815352/1‑150 (MQ=255)
          gCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcggta                                                                                                                                        <  2:915262/150‑5 (MQ=255)
          gCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcggta                                                                                                                                        <  1:1138448/150‑5 (MQ=255)
            aGATGGGGATGGTTGGATTGTACCTCCCGCATCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGGCGCGCTTACCGCCGCTTTTAATGCGGcggtaag                                                                                                                                      <  2:741127/150‑7 (MQ=255)
            aGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCTTTTAATGCGGcggtaag                                                                                                                                      <  1:546839/150‑7 (MQ=255)
            aGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCg                                                                                                                                      <  2:712809/150‑1 (MQ=255)
             gATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCTCTCTTACCGCCTCATTAAAAGCGTCAGGTCGt                                                                                                                                     <  2:1592146/150‑1 (MQ=255)
              aTGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTa                                                                                                                                    <  1:424129/150‑1 (MQ=255)
                ggggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTaaa                                                                                                                                  <  2:664755/150‑1 (MQ=255)
                 gggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAg                                                                                                                                 >  1:297009/1‑150 (MQ=255)
                  ggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGt                                                                                                                                <  1:990292/150‑1 (MQ=255)
                  ggATGGGTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGt                                                                                                                                <  2:987047/150‑1 (MQ=255)
                   gATGGTTGGATTGTACCTCCCGCAGCCGCTCCGTGGGACTTCGGATTGATTGTTCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTg                                                                                                                               <  2:673787/150‑1 (MQ=255)
                   gATGGTTGGAGTGTACCTCCCGCAGCCGCGCCGGGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTg                                                                                                                               <  2:888056/150‑1 (MQ=255)
                     tGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTg                                                                                                                               <  1:1026401/148‑1 (MQ=255)
                     tGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTg                                                                                                                               >  2:1026401/1‑148 (MQ=255)
                     tGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGAc                                                                                                                             <  1:404798/150‑1 (MQ=255)
                          ggATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcagg                                                                                                                        <  1:422712/150‑1 (MQ=255)
                                 aCCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTCGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTa                                                                                                                 >  2:1179311/1‑150 (MQ=255)
                                         cAGCCGCGCCGTGGGAGTTCTGATTGATTGTGCTACGTAATCCCCCTTCTGGGATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAg                                                                                                         <  2:721151/150‑1 (MQ=255)
                                         cAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAg                                                                                                         <  2:196216/150‑1 (MQ=255)
                                            ccGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatc                                                                                                      <  2:752174/150‑1 (MQ=255)
                                             cgcgCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatca                                                                                                     <  2:1254037/150‑1 (MQ=255)
                                             cgcgCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatca                                                                                                     <  2:17581/150‑1 (MQ=255)
                                               cgcCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatcatc                                                                                                   >  2:1027023/1‑150 (MQ=255)
                                                  cGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCg                                                                                                >  1:1027635/1‑150 (MQ=255)
                                                     gggACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATa                                                                                             >  1:1625040/1‑150 (MQ=255)
                                                        aCTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGt                                                                                          >  2:124630/1‑150 (MQ=255)
                                                           tCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGt                                                                                       >  2:1586190/1‑150 (MQ=255)
                                                           tCGGATTGATTGTGCTAAGTAATACCCCTTCTGGCATTACGCCCTAGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCAGCATTAAAAGCGGCAGCTCTTAAAGTGCCTCAGGCAGGCAACCCTGAACATCATCTCGATACGTTGt                                                                                       >  2:1468912/1‑150 (MQ=255)
                                                             ggtttgattgTGCTACTTAATCCCCCTTCTGTCATTACGCCGTTGCCGTTATTTGAGGGAGATAAATCCGCGCTTAGCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTAc                                                                                     <  2:614617/147‑1 (MQ=255)
                                                                  gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                                >  1:652010/1‑150 (MQ=255)
                                                                  gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                                <  1:6633/150‑1 (MQ=255)
                                                                  gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                                >  2:1440834/1‑150 (MQ=255)
                                                                      gtgCTACTTAATCCCCCTTCTGGCATTACTCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAATTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCaa                                                                            >  2:1602666/1‑150 (MQ=255)
                                                                        gCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCaaaa                                                                          >  2:93019/1‑150 (MQ=255)
                                                                         cTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAAc                                                                         <  1:1027023/150‑1 (MQ=255)
                                                                          tACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACt                                                                        <  2:571095/150‑1 (MQ=255)
                                                                            cGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgt                                                                      >  2:82572/1‑150 (MQ=255)
                                                                             gTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgtg                                                                     <  2:623497/150‑1 (MQ=255)
                                                                             gTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgtg                                                                     <  1:1249401/150‑1 (MQ=255)
                                                                              tAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGa                                                                    <  1:1341777/150‑1 (MQ=255)
                                                                              tAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACAGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGa                                                                    >  1:1254665/1‑150 (MQ=255)
                                                                                 tCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGt                                                                 >  1:1302930/1‑150 (MQ=255)
                                                                                         cTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActggct                                                         >  2:669477/1‑150 (MQ=255)
                                                                                          tGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActggctg                                                        >  1:434493/1‑150 (MQ=255)
                                                                                            gCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTggg                                                      <  1:669477/150‑1 (MQ=255)
                                                                                              aTTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGgcg                                                    <  1:1575404/150‑1 (MQ=255)
                                                                                                  cGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCaaa                                                >  1:1505454/1‑150 (MQ=255)
                                                                                                   gccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAAc                                               <  1:18078/150‑1 (MQ=255)
                                                                                                    ccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGCCTGGCTGGGCGCAAACg                                              >  2:802981/1‑150 (MQ=255)
                                                                                                    ccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACg                                              >  2:797089/1‑150 (MQ=255)
                                                                                                     cgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGt                                             >  2:769812/1‑150 (MQ=255)
                                                                                                       ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                           >  1:1476606/1‑150 (MQ=255)
                                                                                                        tgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTTAAGTGACTGGCTGGGCGCAAACGTCgg                                          <  1:124630/150‑1 (MQ=255)
                                                                                                          ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTg                                        <  2:652010/150‑1 (MQ=255)
                                                                                                          ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTg                                        >  1:717828/1‑150 (MQ=255)
                                                                                                          ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACGGTGAAGTGACTGGCTGGGCGCAAACGTCGGTg                                        <  2:1361588/150‑1 (MQ=255)
                                                                                                           cgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGa                                       >  1:689032/1‑150 (MQ=255)
                                                                                                               aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTc                                   <  2:60874/150‑1 (MQ=255)
                                                                                                                  tGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTcaca                                >  2:1195457/1‑150 (MQ=255)
                                                                                                                     gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCa                             >  1:245508/1‑150 (MQ=255)
                                                                                                                     gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCa                             >  1:805496/1‑150 (MQ=255)
                                                                                                                     gggAGATAAAGACGCGCTTACCGCCGCATTAAAACCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGCCTGGCTGGGCGCAACCGTCGGTGATGTCACATCa                             >  1:245503/1‑150 (MQ=255)
                                                                                                                         gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGc                         <  1:769812/150‑1 (MQ=255)
                                                                                                                           tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGa                       >  2:172109/1‑150 (MQ=255)
                                                                                                                           tAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGa                       >  2:988508/1‑150 (MQ=255)
                                                                                                                             aaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATa                     <  2:1302930/150‑1 (MQ=255)
                                                                                                                             aaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATa                     >  2:985147/1‑150 (MQ=255)
                                                                                                                                    gcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGcc              <  2:1625040/150‑1 (MQ=255)
                                                                                                                                      gacctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCgg            >  1:1090399/6‑150 (MQ=255)
                                                                                                                                       tACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCg                        >  1:1272852/1‑137 (MQ=255)
                                                                                                                                       tACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCg                        <  2:1272852/137‑1 (MQ=255)
                                                                                                                                        aCCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTTCAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGt          <  1:449706/150‑1 (MQ=255)
                                                                                                                                         ccgccgCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTg         >  1:125984/1‑150 (MQ=255)
                                                                                                                                              gcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGc    <  1:1586190/150‑1 (MQ=255)
                                                                                                                                                aTTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCgg  >  1:1482469/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGG  >  NC_000913/1672337‑1672630

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