Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,403,948 A→T 5.3% D108V (GAC→GTC)  PP_3014 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,403,9480AT 5.3% 266.9 / 3.0 87D108V (GAC→GTC) PP_3014hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base T (2/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCGGATAGCT  >  NC_002947/3403806‑3404096
                                                                                                                                              |                                                                                                                                                     
ccACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTa                                                                                                                                                >  1:1493820/1‑150 (MQ=255)
ccACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTa                                                                                                                                                <  1:711902/150‑1 (MQ=255)
 cacaGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTAc                                                                                                                                              <  1:2223434/150‑1 (MQ=255)
    aGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGcc                                                                                                                                           >  2:1671242/1‑150 (MQ=255)
     ggCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCa                                                                                                                                          <  2:1733817/150‑1 (MQ=255)
      gCCACGCAGACGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAg                                                                                                                                         >  2:1165903/1‑150 (MQ=255)
      gCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAg                                                                                                                                         <  2:1962896/150‑1 (MQ=255)
      gCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAg                                                                                                                                         <  1:1839850/150‑1 (MQ=255)
       ccACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAgg                                                                                                                                        >  1:1009194/1‑150 (MQ=255)
           gCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTATGCCAGGCTGt                                                                                                                                    >  1:366502/1‑150 (MQ=255)
             aGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCa                                                                                                                                  >  1:97629/1‑150 (MQ=255)
             aGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGGGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAATTTACGCCACGCTGTCa                                                                                                                                  >  2:1225022/1‑150 (MQ=255)
              gAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAg                                                                                                                                 >  1:196863/1‑150 (MQ=255)
                    gATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATg                                                                                                                           >  1:18070/1‑150 (MQ=255)
                     aTAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGc                                                                                                                          >  2:886078/1‑150 (MQ=255)
                      tAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGcc                                                                                                                         <  2:97629/150‑1 (MQ=255)
                       aGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCt                                                                                                                        <  2:1215021/150‑1 (MQ=255)
                        gCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTg                                                                                                                       <  1:1735473/150‑1 (MQ=255)
                          ccACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGa                                                                                                                     <  1:1441569/150‑1 (MQ=255)
                            aCGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGa                                                                                                                                                       <  1:2201910/115‑1 (MQ=255)
                            aCGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGa                                                                                                                                                       >  2:2201910/1‑115 (MQ=255)
                              ggCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGa                                                                                                                 <  1:984304/150‑1 (MQ=255)
                              ggCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGa                                                                                                                 <  1:984327/150‑1 (MQ=255)
                                  aCCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCa                                                                                                             <  2:1928586/150‑1 (MQ=255)
                                  aCCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCa                                                                                                             <  1:1165903/150‑1 (MQ=255)
                                   ccGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAg                                                                                                            <  1:1866799/150‑1 (MQ=255)
                                   ccGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAg                                                                                                            <  1:2080247/150‑1 (MQ=255)
                                    cGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCCTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGt                                                                                                           <  2:1261496/150‑1 (MQ=255)
                                        tCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGAcc                                                                                                       >  1:215292/1‑150 (MQ=255)
                                             gcgcCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGc                                                                                                  <  2:18070/150‑1 (MQ=255)
                                                 cAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGt                                                                                              >  2:1032967/1‑150 (MQ=255)
                                                      gcCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTgc                                                                                         >  1:734185/1‑150 (MQ=255)
                                                      gcCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTgc                                                                                         <  2:1009194/150‑1 (MQ=255)
                                                              aaGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGa                                                                                 >  2:151466/1‑150 (MQ=255)
                                                                       tgtgCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCa                                                                        <  1:1032967/150‑1 (MQ=255)
                                                                           cTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTtgt                                                                                           <  1:94470/127‑1 (MQ=255)
                                                                           cTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTtgt                                                                                           >  2:94470/1‑127 (MQ=255)
                                                                           cTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTa                                                                    >  2:1221696/1‑150 (MQ=255)
                                                                             tGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTAcc                                                                  <  1:151466/150‑1 (MQ=255)
                                                                                tgtgGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCa                                                                                                                                          >  2:1410796/1‑75 (MQ=255)
                                                                                tgtgGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCa                                                                                                                                          <  1:1410796/75‑1 (MQ=255)
                                                                                 gtgGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCaa                                                              <  1:609910/150‑1 (MQ=255)
                                                                                    gTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTATAAGGTCAGCTTTCTGTCTGGCGACAAGTTACGGTATGCTGTCAGGCTACTCCTGGATTGAATTAGTTACCCGAGCGGGTTGTGCCGCCCCGATCGACTGCTGCGACTCCAACTc                                                           >  1:1355711/1‑150 (MQ=255)
                                                                                       gACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATc                                                        >  2:1874085/1‑150 (MQ=255)
                                                                                        aCGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCTACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATcc                                                       >  1:2295222/1‑150 (MQ=255)
                                                                                            tCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCt                                                   >  2:2283424/1‑150 (MQ=255)
                                                                                             cACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCtt                                                  <  2:1493820/150‑1 (MQ=255)
                                                                                                gAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCTCCTGCAGTTACCCCAACTCATCCCGCTTGtt                                               >  2:1924849/1‑150 (MQ=255)
                                                                                                    gCAGCCGCGCGAGCTTTTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGACTCGGTGACCCGAGCGGGTTGTGCCGACGAGGCCGCCTGCAGTTACCCCTTCGCAGCCCGTTTGTTCaaa                                           <  2:2449883/150‑1 (MQ=255)
                                                                                                     cAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAAt                                          <  2:215292/150‑1 (MQ=255)
                                                                                                         cgcgcgACCTTCTACAAGCGCAGCTTTCTCGCTGGCGCCCAGTTAAGCCAGGCTGTCAGGCAATGCCTGGACTCAATCAGTGCCCAGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTc                                      <  1:624008/150‑1 (MQ=255)
                                                                                                            gcgACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCa                                   <  2:196863/150‑1 (MQ=255)
                                                                                                            gcgACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCa                                   >  2:2202933/1‑150 (MQ=255)
                                                                                                              gACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACt                                 >  1:1630130/1‑150 (MQ=255)
                                                                                                                  ttCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCg                                     >  1:38064/1‑142 (MQ=255)
                                                                                                                  ttCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCg                                     <  2:38064/142‑1 (MQ=255)
                                                                                                                  ttCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCaa                             >  1:818239/1‑149 (MQ=255)
                                                                                                                    tccaggcattgcctgacagcCTGGCGTAACTTGTCGCCAGGCTGTCGGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCt                           <  1:1134867‑M1/130‑1 (MQ=255)
                                                                                                                       cAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTgcg                        >  1:780468/1‑150 (MQ=255)
                                                                                                                        aaGGGCAGCTTTCTGTCTGGCTACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTgcgc                       <  2:2295222/150‑1 (MQ=255)
                                                                                                                          gggCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCgg                                                                                                <  1:533243/75‑1 (MQ=255)
                                                                                                                          gggCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCgg                                                                                                >  2:533243/1‑75 (MQ=255)
                                                                                                                          cattgcctgacagcCTGGCGTAACTTGTCGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGtt                                                                                             >  1:1000554‑M1/15‑78 (MQ=255)
                                                                                                                          cattgcctgacagcCTGGCGTAACTTGTCGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGtt                                                                                             >  1:1000616‑M1/15‑78 (MQ=255)
                                                                                                                          cattgcctgacagcCTGGCGTAACTTGTCGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGtt                                                                                             <  2:1000616‑M1/64‑1 (MQ=255)
                                                                                                                          cattgcctgacagcCTGGCGTAACTTGTCGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGtt                                                                                             <  2:1000554‑M1/64‑1 (MQ=255)
                                                                                                                            gCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTcc                   <  2:1630130/150‑1 (MQ=255)
                                                                                                                             cAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCa                  >  1:941680/1‑150 (MQ=255)
                                                                                                                              aGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCAt                 >  1:270474/1‑150 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCTCCGATCGCCTGCa                                                                        >  1:1333980/1‑93 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCa                                                                        <  2:1333980/93‑1 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTc               >  1:26747/1‑150 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTc               <  1:2283424/150‑1 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTc               <  1:2365239/150‑1 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTc               >  1:824043/1‑150 (MQ=255)
                                                                                                                                cTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGCGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTc               >  1:1173465/1‑150 (MQ=255)
                                                                                                                                      gtctgGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTAc                                                                   >  2:1124047/1‑92 (MQ=255)
                                                                                                                                      gtctgGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTAc                                                                   <  1:1124047/92‑1 (MQ=255)
                                                                                                                                       tctgGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCgg        >  2:340232/1‑150 (MQ=255)
                                                                                                                                         taGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCGGAt      <  1:442584‑M1/148‑1 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTAccc                                                                 <  2:1486221‑M1/91‑1 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTAccc                                                                 >  1:1486221‑M1/1‑91 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAAc                                                             <  2:1950598‑M1/95‑1 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAAc                                                             >  1:1950598‑M1/1‑95 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTc                                      >  2:2094185‑M1/1‑118 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTc                                      <  1:2094185‑M1/118‑1 (MQ=255)
                                                                                                                                         TGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCGGAt      <  1:2310467‑M1/150‑1 (MQ=255)
                                                                                                                                           GCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCGGATAg    <  1:1267069‑M1/150‑1 (MQ=255)
                                                                                                                                             GACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCGGATAGCt  >  2:26747‑M1/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CCACAGGCCACGCAGAAGTCGATAGCCCACGGCTACCGGCTCACGGCGCCAAGCGCCGTGGCAAGCGACGATGTGCTTGATGTGGTTGACGCTCACGAGGGCAGCCGCGCGACCTTCTACAAGGGCAGCTTTCTGTCTGGCGACAAGTTACGCCAGGCTGTCAGGCAATGCCTGGATTGAATCAGTGACCCGAGCGGGTTGTGCCGCCCCGATCGCCTGCAGTTACCCCAACTCATCCCGCTTGTTCAAATGTTCGCACTTCACCTGCGCCTCCATTCGGCTCGGATAGCT  >  NC_002947/3403806‑3404096

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