Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 702,598 A→G S52P (TCC→CCC)  nagA ← N‑acetylglucosamine‑6‑phosphate deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913702,5980AG98.9% 337.4 / ‑3.9 90S52P (TCC→CCC) nagAN‑acetylglucosamine‑6‑phosphate deacetylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (34/54);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTAATGC  >  NC_000913/702455‑702745
                                                                                                                                               |                                                                                                                                                   
aaGCGTCGGCAGATAGTTAGTCCAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTGCAGCGGTGTCGTTAAACTGTACGCCGCCGTAGCCGTTTAACTGCACATCGATAAAGCCGGGGGGGAGAAt                                                                                                                                               <  2:1147591/150‑1 (MQ=255)
 aGCTTCGGCAGATAGTTAGTACATCCTGATTTCTCATTTGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATg                                                                                                                                              <  2:434369/150‑1 (MQ=255)
 aGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATgg                                                                                                                                             >  2:907730/1‑151 (MQ=255)
     tCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGcccc                                                                                                                                         >  2:734795/1‑151 (MQ=255)
      cGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCg                                                                                                                                        >  1:172946/1‑151 (MQ=255)
        ggagatagatAGGACAGCCTGATGACACATGGGCTTTCTGCATGATTTCTAGCGTTTCCACGCTGCCCGCTTCAACGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTCACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGt                                                                                                                                       >  2:766726/3‑150 (MQ=255)
        gCAGATAGTTAGTGCAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGtt                                                                                                                                      <  1:374077/151‑1 (MQ=255)
         cAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTc                                                                                                                                     <  2:1171515/151‑1 (MQ=255)
          aGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTc                                                                                                                                     >  1:731445/1‑150 (MQ=255)
          aGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTc                                                                                                                                     >  2:155940/1‑150 (MQ=255)
          aGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCa                                                                                                                                    <  1:1166631/151‑1 (MQ=255)
           gATAGTTAGTACAGCCTTATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAg                                                                                                                                   <  2:1065287/151‑1 (MQ=255)
             tagttagtACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAg                                                                                                                                   >  2:830602/1‑149 (MQ=255)
               gttagtACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGaa                                                                                                                               <  1:482240/151‑1 (MQ=255)
                 tagtACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACg                                                                                                                             >  1:912708/1‑151 (MQ=255)
                       aGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTc                                                                                                                       >  1:189323/1‑151 (MQ=255)
                        gCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCg                                                                                                                      >  1:1215227/1‑151 (MQ=255)
                           tGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGAtc                                                                                                                   >  2:67857/1‑151 (MQ=255)
                            gATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGAtc                                                                                                                   <  1:500577/150‑1 (MQ=255)
                                tctcATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTg                                                                                                               <  1:1120557/150‑1 (MQ=255)
                                   cATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTggcg                                                                                                            <  2:1207909/150‑1 (MQ=255)
                                      tGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTggcggc                                                                                                          >  1:1177041/1‑149 (MQ=255)
                                       ggCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGt                                                                                                       <  2:1111368/151‑1 (MQ=255)
                                        gCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGtt                                                                                                      >  2:883366/1‑151 (MQ=255)
                                           ttCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTcc                                                                                                    <  2:274013/150‑1 (MQ=255)
                                           ttCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTcc                                                                                                    <  2:638986/150‑1 (MQ=255)
                                            tCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGc                                                                                                  >  2:941503/1‑151 (MQ=255)
                                             cTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCt                                                                                                 >  2:290349/1‑151 (MQ=255)
                                              tGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCt                                                                                                 >  1:366451/1‑150 (MQ=255)
                                               gCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCt                                                                                                 <  1:461056/149‑1 (MQ=255)
                                                cATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTAcc                                                                                              <  1:802534/151‑1 (MQ=255)
                                                  tGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTGGATCTCTGCCGGCAGTTCCGCTCCCg                                                                                             >  1:643449/1‑150 (MQ=255)
                                                      ttCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACa                                                                                         <  1:378957/150‑1 (MQ=255)
                                                        ccAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGa                                                                                       <  1:471999/150‑1 (MQ=255)
                                                         cAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAACCCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACg                                                                                     <  1:260416/151‑1 (MQ=255)
                                                         cAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACg                                                                                     >  1:945826/1‑151 (MQ=255)
                                                         cAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTCTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGa                                                                                       >  2:210012/1‑149 (MQ=255)
                                                            cGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCtt                                                                                  >  1:535000/1‑151 (MQ=255)
                                                            cGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCtt                                                                                  >  1:854862/1‑151 (MQ=255)
                                                             gTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCttt                                                                                 <  1:468659/151‑1 (MQ=255)
                                                              tttCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCtttt                                                                                >  1:111753/1‑151 (MQ=255)
                                                               ttCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCtttt                                                                                <  2:286533/150‑1 (MQ=255)
                                                                tCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTa                                                                               <  1:584384/150‑1 (MQ=255)
                                                                 ccACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAAt                                                                             >  1:827416/1‑151 (MQ=255)
                                                                  cACGCTGACCGCGTCAGCGGTGTCGTTAAACTTTACGCCGCCGCAGCCGTTTAACTGCATATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGATCGTTGTTCGATCTCTGGCGGTAGTTCCGCTCCCGGACAGACGCTTTTat                                                                              <  2:994412/149‑2 (MQ=255)
                                                                    cGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCa                                                                           <  1:583301/150‑1 (MQ=255)
                                                                    cGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAg                                                                          <  2:1115510/151‑1 (MQ=255)
                                                                     gCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAgg                                                                         <  2:359188/151‑1 (MQ=255)
                                                                     gCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAgg                                                                         >  2:973709/1‑151 (MQ=255)
                                                                       tGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGcc                                                                       >  2:656455/1‑151 (MQ=255)
                                                                        gACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCa                                                                      <  1:1160896/151‑1 (MQ=255)
                                                                            gCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATc                                                                    >  1:632993/1‑149 (MQ=255)
                                                                             cTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCGGGCCATCAg                                                                  <  1:388093/150‑1 (MQ=255)
                                                                              ttCAGCGGTGTCGTTAACCTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGTTGAGAATGGCCGCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGCCAGACGCTTTTAATCAgg                                                                         <  2:1200492/142‑1 (MQ=255)
                                                                              ttCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCg                                                                >  1:1065130/1‑151 (MQ=255)
                                                                               tCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGc                                                                 <  1:459876/149‑1 (MQ=255)
                                                                               tCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGc                                                                 <  1:907730/149‑1 (MQ=255)
                                                                               tCAGCGGTGTCGTTAAACTGGACGCCGCCGCAGCCGTTTGACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGATCGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGc                                                                 <  2:632972/149‑1 (MQ=255)
                                                                                cAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGAt                                                              <  2:292416/151‑1 (MQ=255)
                                                                                 aGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATa                                                             <  2:565453/151‑1 (MQ=255)
                                                                                   cGGTGTCGTTAAACTGTACGCCGCCGCTGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTGCAGTGAGCGTTGTTCGATCTCTGGCGGCATTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATaa                                                            <  2:1207286/150‑1 (MQ=255)
                                                                                     gtgtCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATaacaa                                                         <  2:832800/151‑1 (MQ=255)
                                                                                      tgtCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAACCGCCCCGCTCAGGGAACGTAGTTCGATCTCTGGCGGCAGTTCCGCCACCGGACAAACGCTTTTAATCAGGCCATCAGCGATaacaac                                                        >  2:678561/1‑151 (MQ=255)
                                                                                            tAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCGGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTg                                                  >  1:1090329/1‑151 (MQ=255)
                                                                                             aaaCTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTgg                                                 <  2:661421/151‑1 (MQ=255)
                                                                                             aaaCTGTACGCCGACGCAGCCGTTGAACTGCACATCGATAAAACCGGGAGGGAGAATGGCCCCGTTCAGTGAACGTTGTGCAATCTCTAGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCGGCGATAACAACCGAGTg                                                  <  2:541727/150‑1 (MQ=255)
                                                                                               aCTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGtc                                               <  2:391082/151‑1 (MQ=255)
                                                                                                 tGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGtcat                                             >  1:676993/1‑151 (MQ=255)
                                                                                                  gTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGtca                                              <  2:1236064/149‑1 (MQ=255)
                                                                                                     cgccgccgcAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCaa                                          <  1:347800/150‑1 (MQ=255)
                                                                                                      gccgccgcAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGCCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGccatca                                           >  2:366965/1‑148 (MQ=255)
                                                                                                         gccgcAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAAt                                     <  2:512073/151‑1 (MQ=255)
                                                                                                           cgcAGCCGTTTAACTGCACATCGATAAAGCCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTc                                   <  2:492110/151‑1 (MQ=255)
                                                                                                           cgcAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTc                                   <  2:795987/151‑1 (MQ=255)
                                                                                                            gcAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCg                                  >  2:832621/1‑151 (MQ=255)
                                                                                                             cAGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGt                                 <  1:108550/151‑1 (MQ=255)
                                                                                                              aGCCGTTTAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTg                                >  1:1056365/1‑151 (MQ=255)
                                                                                                                     tAACTGCACATCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGt                          <  1:441035/150‑1 (MQ=255)
                                                                                                                            acaTCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATcc                  <  1:440168/151‑1 (MQ=255)
                                                                                                                             caTCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCg                 <  1:658566/151‑1 (MQ=255)
                                                                                                                              aTCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCgg                <  2:291780/151‑1 (MQ=255)
                                                                                                                              aTCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCgg                <  1:155936/151‑1 (MQ=255)
                                                                                                                              aTCGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCgg                <  2:110105/151‑1 (MQ=255)
                                                                                                                                cGATAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGc               <  1:674839/150‑1 (MQ=255)
                                                                                                                                  aTAAAACCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCt            <  1:317309/151‑1 (MQ=255)
                                                                                                                                    aaaaCCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTgg          >  2:134855/1‑151 (MQ=255)
                                                                                                                                     aaaCCGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTggg         <  1:835040/151‑1 (MQ=255)
                                                                                                                                         cGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTaa     <  1:734798/151‑1 (MQ=255)
                                                                                                                                         cGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTaa     <  1:14582/151‑1 (MQ=255)
                                                                                                                                         cGGGGGGGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTa      >  2:646638/1‑150 (MQ=255)
                                                                                                                                          ggggggGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTAAt    <  2:1028835/151‑1 (MQ=255)
                                                                                                                                          ggggggGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTAAt    >  2:1223711/1‑151 (MQ=255)
                                                                                                                                            ggggGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTAATGc  <  2:534977/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AAGCGTCGGCAGATAGTTAGTACAGCCTGATTTCTCATTGGCTTTCTGCATGATTTCCAGCGTTTCCACGCTGACCGCTTCAGCGGTGTCGTTAAACTGTACGCCGCCGCAGCCGTTTAACTGCACATCGATAAAACCGGGGGAGAGAATGGCCCCGTTCAGTGAACGTTGTTCGATCTCTGGCGGCAGTTCCGCTACCGGACAGACGCTTTTAATCAGGCCATCAGCGATAACAACCGCGTGGTCATCAAGAAATTCGTGGCCGGTAAAGATCCGGCCCTGGGTTAATGC  >  NC_000913/702455‑702745

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