Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,336,371 G→T 5.0% Q42K (CAG→AAG)  murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,336,3710GT 5.0% 267.8 / 2.6 79Q42K (CAG→AAG) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (44/31);  new base T (2/2);  total (46/33)
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.77e-01

GCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCA  >  NC_000913/3336224‑3336520
                                                                                                                                                   |                                                                                                                                                     
gCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTTGa                                                                                                                                                     >  1:1492295/1‑150 (MQ=255)
gCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTTGa                                                                                                                                                     >  1:1486300/1‑150 (MQ=255)
 cATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAt                                                                                                                                                    <  1:106813/150‑1 (MQ=255)
 cATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAt                                                                                                                                                    <  2:279710/150‑1 (MQ=255)
  aTGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtc                                                                                                                                                   <  1:212104/150‑1 (MQ=255)
    ggTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtctc                                                                                                                                                 >  1:815023/1‑150 (MQ=255)
    ggTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtctc                                                                                                                                                 <  1:478048/150‑1 (MQ=255)
     gTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATctct                                                                                                                                                >  1:979702/1‑150 (MQ=255)
           aCCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGt                                                                                                                                          >  2:926256/1‑150 (MQ=255)
           aCCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGt                                                                                                                                          <  2:595187/150‑1 (MQ=255)
                aTCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtct                                                                                                                                                  <  2:1614522/137‑1 (MQ=255)
                aTCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGAtct                                                                                                                                                  >  1:1614522/1‑137 (MQ=255)
                 tCGTAAGGTGCGCATAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTcc                                                                                                                                    >  2:889446/1‑150 (MQ=255)
                    tAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGcc                                                                                                                                 <  1:1228422/150‑1 (MQ=255)
                           cgcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTg                                                                                                                          <  1:1279537/150‑1 (MQ=255)
                            gcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGc                                                                                                                         <  1:621803/150‑1 (MQ=255)
                            gcAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGc                                                                                                                         >  2:140725/1‑150 (MQ=255)
                              aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTGTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCgg                                                                                                                       <  2:1429067/150‑1 (MQ=255)
                              aGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCgg                                                                                                                       <  2:918839/150‑1 (MQ=255)
                               gAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGc                                                                                                                      >  1:1497123/1‑150 (MQ=255)
                                    cATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCaaaaa                                                                                                                 >  1:1075360/1‑150 (MQ=255)
                                    cATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCaaaaa                                                                                                                 >  2:1272311/1‑150 (MQ=255)
                                    cATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCaaaaa                                                                                                                 >  1:39711/1‑150 (MQ=255)
                                            tCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGc                                                                                                         <  1:1258212/150‑1 (MQ=255)
                                                  gCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGagca                                                                                                   >  1:1501186/1‑150 (MQ=255)
                                                        atatGCACAGAACCATTACGTGCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTGCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCAttt                                                                                             <  2:1447412/150‑1 (MQ=255)
                                                               caGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGCGCAGCATTTTTCgcgc                                                                                      >  2:1476554/1‑150 (MQ=255)
                                                                 gAACCATTACGTTCTTCTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCg                                                                                    >  2:1579562/1‑150 (MQ=255)
                                                                 gAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCg                                                                                    >  2:98720/1‑150 (MQ=255)
                                                                    ccATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGaa                                                                                 >  2:2714/1‑150 (MQ=255)
                                                                          cGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATtgtg                                                                           >  2:994782/1‑150 (MQ=255)
                                                                          cGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATtgtg                                                                           <  2:1052191/150‑1 (MQ=255)
                                                                          cGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATtgtg                                                                           <  2:514393/150‑1 (MQ=255)
                                                                            ttCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGaa                                                                         >  1:354182/1‑149 (MQ=255)
                                                                            ttCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGAc                                                                         <  1:110352/150‑1 (MQ=255)
                                                                              cTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACtt                                                                       >  1:1479682/1‑150 (MQ=255)
                                                                               tACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTc                                                                      >  1:951715/1‑150 (MQ=255)
                                                                                aCTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCg                                                                     >  1:1316615/1‑150 (MQ=255)
                                                                                aCTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCg                                                                     <  2:232167/150‑1 (MQ=255)
                                                                                  tttCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGcc                                                                   <  2:1497123/150‑1 (MQ=255)
                                                                                  tttCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGcc                                                                   >  1:1363720/1‑150 (MQ=255)
                                                                                   ttCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCGCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGccc                                                                  <  2:1632913/150‑1 (MQ=255)
                                                                                      gCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTgg                                                               >  2:210433/1‑150 (MQ=255)
                                                                                       cACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGa                                                              >  1:489439/1‑150 (MQ=255)
                                                                                         cccAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGc                                                            >  1:1529334/1‑150 (MQ=255)
                                                                                                ggCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCATAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGGCCTGGAGCTTCGTTg                                                     >  2:1047897/1‑150 (MQ=255)
                                                                                                 gCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTgg                                                    >  1:1226309/1‑150 (MQ=255)
                                                                                                    tAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTTGATCTCTACCGGTTCTTCCGCCAGTAGTTCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGccc                                                 <  2:1492295/150‑1 (MQ=255)
                                                                                                    tAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTTGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGccc                                                 <  2:1486300/150‑1 (MQ=255)
                                                                                                     agcagcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGcccc                                                <  1:864329/150‑1 (MQ=255)
                                                                                                      gcagcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCt                                               >  1:298691/1‑150 (MQ=255)
                                                                                                      gcagcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCt                                               >  1:487815/1‑150 (MQ=255)
                                                                                                       aagcTTCAATGAAGTATCCACGTCTTTCAGTTTCGGGACGTTCTGGATATATACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTg                                              >  2:231996/2‑150 (MQ=255)
                                                                                                        agcTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGa                                             <  2:826128/150‑1 (MQ=255)
                                                                                                            tCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAacac                                         >  2:1084970/1‑150 (MQ=255)
                                                                                                            tCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCAGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAacac                                         >  1:715938/1‑150 (MQ=255)
                                                                                                              aTTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGa                                       <  2:815023/150‑1 (MQ=255)
                                                                                                                 gATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAAt                                    >  1:498648/1‑150 (MQ=255)
                                                                                                                      aTCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATc                               >  1:668437/1‑150 (MQ=255)
                                                                                                                       tCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATcc                              >  1:11401/1‑150 (MQ=255)
                                                                                                                        cGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCa                             >  2:616167/1‑150 (MQ=255)
                                                                                                                        cGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCa                             >  2:110797/1‑150 (MQ=255)
                                                                                                                         gACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAt                            >  1:1379714/1‑150 (MQ=255)
                                                                                                                          aCGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAtt                           <  2:979702/150‑1 (MQ=255)
                                                                                                                           cGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCAttt                          <  1:616167/150‑1 (MQ=255)
                                                                                                                             tCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAg                        >  1:1524638/1‑150 (MQ=255)
                                                                                                                              cTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCCCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCAATTTAGt                       >  2:517084/1‑150 (MQ=255)
                                                                                                                               tttCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGtt                      <  2:1347768/150‑1 (MQ=255)
                                                                                                                                ttCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGttt                     <  2:463449/150‑1 (MQ=255)
                                                                                                                                      ttCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTAt                                <  1:1229990/133‑1 (MQ=255)
                                                                                                                                      ttCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTtctc               <  2:1219815/150‑1 (MQ=255)
                                                                                                                                      ttCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTtctc               <  1:768352/150‑1 (MQ=255)
                                                                                                                                      ttCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTtctc               >  1:1219815/1‑150 (MQ=255)
                                                                                                                                      ttCGGGACGTTCTGGATCTCTACAGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTAt                                >  2:1229990/1‑133 (MQ=255)
                                                                                                                                            aCGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTaa         >  1:1003437/1‑150 (MQ=255)
                                                                                                                                             cGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAAc        >  2:842622/1‑150 (MQ=255)
                                                                                                                                              gTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACa       <  1:926256/150‑1 (MQ=255)
                                                                                                                                               ttCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACaa      <  2:1316615/150‑1 (MQ=255)
                                                                                                                                                 cTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAAtt    >  2:1330596/1‑150 (MQ=255)
                                                                                                                                                   ggATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCa  <  1:110797/150‑1 (MQ=255)
                                                                                                                                                   ggATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCa  <  1:210433/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCATGGTTTTAACCAGATCGTAAGGTGCGCAGAATACATTAACGTCGCGGGCATCAATATGCACAGAACCATTACGTTCTACTTTCGCACCCAGCTGGCTTAGCAGCTTCATTGATGTATCGACGTCTTTCAGTTTCGGGACGTTCTGGATCTCTACCGGTTCTTCCGCCAGTAGTGCGGCAAAAAGGATAGGCAGAGCAGCATTTTTAGCGCCGGAAATTGTGACTTCGCCCTGGAGCTTCGTTGGCCCCTGAACACGAAATTTATCCATTTAGTTTGTTCTCAGTTAACAATTCA  >  NC_000913/3336224‑3336520

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: