Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,719,444 C→T Q53Q (CAG→CAA anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,719,4440CT97.4% 265.1 / ‑1.3 78Q53Q (CAG→CAAanmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (42/34);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

CCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGC  >  NC_000913/1719297‑1719589
                                                                                                                                                   |                                                                                                                                                 
ccAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTgg                                                                                                                                                >  1:887987/1‑151 (MQ=255)
 cAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGc                                                                                                                                               >  1:978633/1‑151 (MQ=255)
  aGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGc                                                                                                                                               >  2:131334/1‑150 (MQ=255)
  aGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCa                                                                                                                                              <  1:649496/151‑1 (MQ=255)
    aCGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCaaa                                                                                                                                            <  1:226636/151‑1 (MQ=255)
     cGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAat                                                                                                                                           <  2:820097/151‑1 (MQ=255)
     cGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAat                                                                                                                                           <  1:620069/151‑1 (MQ=255)
      ggTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAat                                                                                                                                           >  2:663794/1‑150 (MQ=255)
       gTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAatat                                                                                                                                         <  2:873057/151‑1 (MQ=255)
       gTTTGACCGTCACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGAAAGGCATTGACCGCATCAGGAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAGCTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAata                                                                                                                                          <  1:249476/150‑1 (MQ=255)
           gACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAg                                                                                                                                     <  1:724585/151‑1 (MQ=255)
            aCCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGt                                                                                                                                    <  1:775517/151‑1 (MQ=255)
             ccGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGt                                                                                                                                    >  1:143958/1‑150 (MQ=255)
               gTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTAc                                                                                                                                  >  1:996532/1‑150 (MQ=255)
                tGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACa                                                                                                                                 <  2:388591/150‑1 (MQ=255)
                 gACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACa                                                                                                                                 <  2:410861/149‑1 (MQ=255)
                       cGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGttt                                                                                                                         >  2:856177/1‑151 (MQ=255)
                         aTCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGttt                                                                                                                         >  1:30060/1‑149 (MQ=255)
                                  atatCCCTTGCCTTCAGATTTTGTACCTTAAGCAAGGCATTTACCGCATCAGCAAAAAGTTGCCTGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATACAGTACAGCCTGTTTCAGAGATAcc                                                                                                               >  1:439104/1‑150 (MQ=255)
                                     tctcTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCggg                                                                                                            >  2:211044/1‑150 (MQ=255)
                                      ctctTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGAt                                                                                                          >  1:78625/1‑151 (MQ=255)
                                       tctTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATc                                                                                                         >  2:890021/1‑151 (MQ=255)
                                        ctTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACTGGGATCg                                                                                                        >  2:210241/1‑151 (MQ=255)
                                             cTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCa                                                                                                    >  2:383788/1‑150 (MQ=255)
                                               gCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACt                                                                                                 <  2:85830/151‑1 (MQ=255)
                                                   aTTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCaaa                                                                                             <  1:999286/151‑1 (MQ=255)
                                                    ttttGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAAc                                                                                            >  2:797988/1‑151 (MQ=255)
                                                      ttGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACtt                                                                                          >  2:985840/1‑151 (MQ=255)
                                                        gTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGc                                                                                        >  1:371906/1‑151 (MQ=255)
                                                         ttCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGcc                                                                                       >  2:379644/1‑151 (MQ=255)
                                                         ttCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGcc                                                                                       >  1:558176/1‑151 (MQ=255)
                                                           ccTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCa                                                                                      <  1:804530/150‑1 (MQ=255)
                                                           ccTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCa                                                                                      <  1:306062/150‑1 (MQ=255)
                                                             ttAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCt                                                                                   >  2:1006270/1‑151 (MQ=255)
                                                              tAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtg                                                                                  >  1:994108/1‑151 (MQ=255)
                                                               aaGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtgt                                                                                 >  1:808455/1‑151 (MQ=255)
                                                                aGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtgtg                                                                                <  2:416466/151‑1 (MQ=255)
                                                                aGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCtgtg                                                                                >  2:524887/1‑151 (MQ=255)
                                                                  cAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCg                                                                              >  2:835306/1‑151 (MQ=255)
                                                                        aTTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCgg                                                                        >  1:719665/1‑151 (MQ=255)
                                                                          tGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgt                                                                       <  1:1138877/150‑1 (MQ=255)
                                                                           gACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtgt                                                                     >  2:433361/1‑151 (MQ=255)
                                                                           gACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGgtgt                                                                     <  1:854845/151‑1 (MQ=255)
                                                                           gACCGCATCAGCAAAAAGATGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAACCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGAGCGACCCGgt                                                                       >  2:1096035/1‑149 (MQ=255)
                                                                            aCCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGtt                                                                    >  2:428640/1‑151 (MQ=255)
                                                                                cATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtcatc                                                                <  2:97400/151‑1 (MQ=255)
                                                                                 aTCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTtcatc                                                                >  1:1028473/1‑150 (MQ=255)
                                                                                    aGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAAtt                                                            <  2:600711/151‑1 (MQ=255)
                                                                                      cAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGt                                                          <  2:74892/151‑1 (MQ=255)
                                                                                       aaaaaGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTc                                                         >  2:15657/1‑151 (MQ=255)
                                                                                        aaaaGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGGCTTGGCAAATAGCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCg                                                        <  1:511068/151‑1 (MQ=255)
                                                                                        aaaaGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCg                                                        <  1:345640/151‑1 (MQ=255)
                                                                                          aaGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGcc                                                      >  1:356769/1‑151 (MQ=255)
                                                                                               gCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaaca                                                  <  2:686465/150‑1 (MQ=255)
                                                                                                cccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacac                                                 >  2:148502/1‑150 (MQ=255)
                                                                                                cccGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCcaacac                                                 >  2:462237/1‑150 (MQ=255)
                                                                                                    aGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACAt                                            >  2:78260/1‑151 (MQ=255)
                                                                                                      ttGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATc                                           >  1:172732/1‑150 (MQ=255)
                                                                                                           tATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACAcc                                     >  1:667346/1‑151 (MQ=255)
                                                                                                              cAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACCCAACATCAACAcc                                     <  2:396125/148‑1 (MQ=255)
                                                                                                                   tGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGCGTTCATCAATTGTCGCCAACACAACATCAACACCATCaa                                >  1:1185395/1‑148 (MQ=255)
                                                                                                                      ccAAACTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTgg                          >  2:1099623/1‑151 (MQ=255)
                                                                                                                        aaaCTGCGAAAGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTg                        >  1:843101/1‑151 (MQ=255)
                                                                                                                                aaaGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACAt                >  2:930644/1‑151 (MQ=255)
                                                                                                                                  aGTGTAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATa               >  1:568662/1‑150 (MQ=255)
                                                                                                                                   gtgtAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAAc             <  1:889966/151‑1 (MQ=255)
                                                                                                                                     gtAAGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACg            <  2:1130697/150‑1 (MQ=255)
                                                                                                                                       aaGCTGCTGNCCTTGGNNANNATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCACCACCATCAAGGCTGGTGCCTGACATAACGCCa         <  1:1132606/151‑1 (MQ=255)
                                                                                                                                       aaGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGcc          <  2:586564/150‑1 (MQ=255)
                                                                                                                                       aaGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCa         <  2:188846/151‑1 (MQ=255)
                                                                                                                                        aGCTGCTGGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCaa        <  2:309305/151‑1 (MQ=255)
                                                                                                                                            gctgGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATaaa    <  1:949158/151‑1 (MQ=255)
                                                                                                                                            gctgGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATaaa    <  1:337782/151‑1 (MQ=255)
                                                                                                                                            gctgGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATaa     >  2:100791/1‑150 (MQ=255)
                                                                                                                                            gctgGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATa      <  1:867930/149‑1 (MQ=255)
                                                                                                                                             ctgGCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATaa     >  1:89588/1‑149 (MQ=255)
                                                                                                                                               ggCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAg   <  1:131328/149‑1 (MQ=255)
                                                                                                                                               ggCCTTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGc  <  1:318193/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCAGACGGTTTGACCGTGACAACCGATCGCAACTATATCTCTTGCCTGCAGATTTTGTTCCTTAAGCAAGGCATTGACCGCATCAGCAAAAAGTTGCCCGAGTTGAGTATCAAGCTGTCCAAACTGCGAAAGTGTAAGCTGCTGGCCCTGGCAAATATCCAGTACAGCCTGTTTCAGAGATACCGGGATCGGCCAACTCAAACTTGCCAGCTGTGCGACCCGGTGTTCATCAATTGTCGCCAACACAACATCAACACCATCAAGGCTGGTGCCTGACATAACGCCAATAAAGC  >  NC_000913/1719297‑1719589

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