Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,515,787 G→T 7.7% L43I (CTC→ATC) ‡ dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,7870GT 7.7% 374.5 / 15.6 139L43I (CTC→ATC) ‡damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (72/54);  new base T (6/6);  total (78/61)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 144 total aligned reads displayed.

ATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTT  >  NC_000913/3515638‑3515935
                                                                                                                                                     |                                                                                                                                                    
aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                                                                      <  1:1995604/150‑1 (MQ=255)
  aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                    <  1:804552/150‑1 (MQ=255)
  aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                    <  1:85600/150‑1 (MQ=255)
   gTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGACTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTTTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaaa                                                                                                                                                   >  1:1380779/1‑150 (MQ=255)
     aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 <  1:2023118/150‑1 (MQ=255)
     aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 >  2:1889181/1‑150 (MQ=255)
       aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACAcc                                                                                                                                               >  1:1398524/1‑150 (MQ=255)
         cTCGGCTCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGa                                                                                                                                             <  2:1669709/150‑1 (MQ=255)
         cTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGa                                                                                                                                             <  1:714211/150‑1 (MQ=255)
                aattttGTTTCGGTAACAAACTGCTCGGGGGCGGCCTGTACGTACTCATCAGTACGCATATTCACAATTTTATAGAGACTGATCAGGTCGCTATTTATATCGCCAAGTATATAACGATAAACGTCGGTGTTGAGAAACACCGACCCGGCa                                                                                                                                      <  2:1815094/147‑1 (MQ=255)
                     tGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTAc                                                                                                                                 >  2:2053180/1‑150 (MQ=255)
                      gTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCTGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACa                                                                                                                                >  1:2042560/1‑150 (MQ=255)
                      gTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACa                                                                                                                                >  2:1136618/1‑150 (MQ=255)
                        ttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaaa                                                                                                                              >  2:888447/1‑150 (MQ=255)
                         tCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGCAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaaaa                                                                                                                             >  2:128828/1‑150 (MQ=255)
                          cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATATAGACTGATCAGGTCGCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                                                                            <  1:655737/150‑1 (MQ=255)
                          cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAg                                                                                                                            >  2:2001899/1‑150 (MQ=255)
                          cGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCAGGCACCTACAAAAg                                                                                                                            >  1:1890124/1‑150 (MQ=255)
                            ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTGTAGAGACTGATCAGGTCGCTATTTATATCGGCAAGGATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                          <  2:2165441/150‑1 (MQ=255)
                            ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                          <  1:2107217/150‑1 (MQ=255)
                             gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACAGACCCGGCCCCTACAAAAGGCt                                                                                                                         >  1:2161530/1‑150 (MQ=255)
                             gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGAACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCt                                                                                                                         <  2:38121/150‑1 (MQ=255)
                                caaacaGCTCGCGTGCGGCCTGTACGTACACATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAACTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGCATCaa                                                                                                                      >  1:547754/1‑150 (MQ=255)
                                 aaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAc                                                                                                                     <  1:239979/150‑1 (MQ=255)
                                  aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                                                                    >  1:1271141/1‑150 (MQ=255)
                                  aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                                                                    >  2:892819/1‑150 (MQ=255)
                                    ccgCTCGGGGGCGGCCTTTACGTACTCATCTGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCAATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAg                                                                                                                  <  1:527031/148‑1 (MQ=255)
                                    caGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAg                                                                                                                  >  1:1206215/1‑150 (MQ=255)
                                      gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGAc                                                                                                                >  1:20673/1‑150 (MQ=255)
                                      gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGAc                                                                                                                >  1:1171992/1‑150 (MQ=255)
                                        tCGCGTGCGGCCTGTACGTAATCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAt                                                                                                              <  2:333413/150‑1 (MQ=255)
                                           cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTCTAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                           <  2:1162349/150‑1 (MQ=255)
                                           cgTGCGGCCTGGACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                           <  2:173217/150‑1 (MQ=255)
                                                ggCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCtt                                                                                                      >  2:1047450/1‑150 (MQ=255)
                                                   cTGTAGTTACTCACCTGTAAGCATGTTCACAATGTTATACAGACTCTTCAGGGCGCTATTTATATGGGCACGGTTATAAGGTGAAATTTCGGTGTTGTCAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTggg                                                                                                   <  2:726058/150‑1 (MQ=255)
                                                    tGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGc                                                                                                  >  2:1175737/1‑150 (MQ=255)
                                                        cgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAAt                                                                                              >  1:982368/1‑150 (MQ=255)
                                                              catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                        >  1:2168576/1‑150 (MQ=255)
                                                              catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                        >  2:1880807/1‑150 (MQ=255)
                                                                  aGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAACGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTATAat                                                                                    >  1:1336910/1‑150 (MQ=255)
                                                                  aGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAat                                                                                    <  2:166017/150‑1 (MQ=255)
                                                                   gTACGCATCTTCACAATGTTATCGAGACTGATCAGGGCGCTATTGATAACGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata                                                                                   <  1:601096/150‑1 (MQ=255)
                                                                   gTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAata                                                                                   <  1:321331/150‑1 (MQ=255)
                                                                    tACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTTGGCAAATGCCGTTTAaaat                                                                                  >  1:1940007/1‑150 (MQ=255)
                                                                        taTCTTAACAATGTTATAGATACTGTTCAGGTCGCTATTGATATCGGCAAGTATATAACGAGATAAGTCGGTGTTGTTAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcatc                                                                              <  1:1889181/149‑1 (MQ=255)
                                                                          tCTTCACAATGTTATAGAGACAGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCaa                                                                            >  1:173106/1‑150 (MQ=255)
                                                                           cTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAg                                                                           >  2:1717979/1‑150 (MQ=255)
                                                                            ttCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                                                                          >  2:803251/1‑150 (MQ=255)
                                                                             tCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCACGCa                                                                         >  2:385106/1‑150 (MQ=255)
                                                                              cacaATGTTATAGAGACTGATCAGGTCGCTATTTATATCGGCAAGGATATAACGAGAAAAGTCGGAGATGAGAAACACCGACCCGCCACCAACAAAAGGCTCAACCACACAAACGCCCATCGGCAAATGCCGTAAAATATCATCAAGCac                                                                        >  2:857850/1‑149 (MQ=255)
                                                                                caATGTTATATAGACTGATAAGGTCGCTATTGATATGGGCAAGGATATATCGAGAAAAGTCGGTTTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCTCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                                                                      <  2:298478/150‑1 (MQ=255)
                                                                                   tGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGAt                                                                   >  2:433152/1‑150 (MQ=255)
                                                                                    gTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATa                                                                  >  2:639341/1‑150 (MQ=255)
                                                                                     ttATAGAGCCTGATCAGGTCGCTATTGATATCGGCAAGGATAAAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATAc                                                                 <  1:1772816/150‑1 (MQ=255)
                                                                                      tataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACt                                                                >  2:1356403/1‑150 (MQ=255)
                                                                                       ataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt                                                               >  1:1207732/1‑150 (MQ=255)
                                                                                           agaCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc                                                           <  1:264843/150‑1 (MQ=255)
                                                                                            gaCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAAGATCATCAAGCAGGGGATACTTGcccc                                                          >  2:435704/1‑150 (MQ=255)
                                                                                                gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc                                                      >  1:989177/1‑150 (MQ=255)
                                                                                                  tCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGccc                                                    >  2:2014326/1‑150 (MQ=255)
                                                                                                    aGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc                                                  <  2:283133/150‑1 (MQ=255)
                                                                                                    aGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc                                                  >  1:941979/1‑150 (MQ=255)
                                                                                                     ggTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCGCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACt                                                 >  1:61294/1‑150 (MQ=255)
                                                                                                     ggTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACt                                                 <  1:2169609/150‑1 (MQ=255)
                                                                                                        cGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCa                                              >  1:101035/1‑150 (MQ=255)
                                                                                                         gCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaa                                             >  2:651236/1‑150 (MQ=255)
                                                                                                           tATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaa                                           <  1:338982/150‑1 (MQ=255)
                                                                                                            aTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaaa                                          <  1:69655/150‑1 (MQ=255)
                                                                                                             ttGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaaaa                                         >  2:713099/1‑150 (MQ=255)
                                                                                                              tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACTCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                        >  1:1009897/1‑150 (MQ=255)
                                                                                                              tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCTTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                        >  2:1581530/1‑150 (MQ=255)
                                                                                                              tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                        <  1:892819/150‑1 (MQ=255)
                                                                                                              tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                        <  1:1136618/150‑1 (MQ=255)
                                                                                                               gATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                       >  2:1992466/1‑150 (MQ=255)
                                                                                                                 tatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgcgc                                     >  2:684528/1‑150 (MQ=255)
                                                                                                                  atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                    <  2:626459/150‑1 (MQ=255)
                                                                                                                  atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATAATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                    >  2:1688154/1‑150 (MQ=255)
                                                                                                                     ggCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAtt                                 >  2:102267/1‑150 (MQ=255)
                                                                                                                     ggCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGGTTAATATCATCAAGCAGGGGATACTTGCCCCGTGCCCACTTCAAAAAAGCGCGAtt                                 >  1:846138/1‑150 (MQ=255)
                                                                                                                      gCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAttt                                <  2:625450/150‑1 (MQ=255)
                                                                                                                          ggATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttct                            <  2:1171992/150‑1 (MQ=255)
                                                                                                                          ggATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttct                            >  1:1267164/1‑150 (MQ=255)
                                                                                                                           gttattACGAGAAAAGTCGGTGTTGAGAAACTCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCATTTGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttctt                           <  2:179687/148‑1 (MQ=255)
                                                                                                                               taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATg                       >  1:269782/1‑150 (MQ=255)
                                                                                                                                aaCGAGAAATGTCGGTGTTGAGAAACACCTACCCGGCACCTACAATAGTCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGc                      <  2:1380779/150‑1 (MQ=255)
                                                                                                                                 aCGAGAAAAGTCGTTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCt                     <  1:1534690/150‑1 (MQ=255)
                                                                                                                                  cGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTg                    <  1:1047450/150‑1 (MQ=255)
                                                                                                                                  cGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGAGCGATTTTTCTTCATGCTg                    >  2:1447400/1‑150 (MQ=255)
                                                                                                                                    agaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGac                  <  1:888447/150‑1 (MQ=255)
                                                                                                                                       cgggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaa               >  1:1861196/4‑150 (MQ=255)
                                                                                                                                         ggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             >  2:2175063/2‑150 (MQ=255)
                                                                                                                                         ggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             >  1:223176/2‑150 (MQ=255)
                                                                                                                                         aggcGGTGTTGATAAACACGGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGAAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             <  1:1642405/147‑1 (MQ=255)
                                                                                                                                            cGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAAt          <  1:433152/150‑1 (MQ=255)
                                                                                                                                              gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa        <  1:966618/150‑1 (MQ=255)
                                                                                                                                              gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa        >  1:654080/1‑150 (MQ=255)
                                                                                                                                               tgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTac       >  2:1296302/1‑150 (MQ=255)
                                                                                                                                                gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTaca      <  2:2078842/150‑1 (MQ=255)
                                                                                                                                                 ttGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTacac     <  2:101035/150‑1 (MQ=255)
                                                                                                                                                   gagaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCt   <  2:1155226/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTT  >  NC_000913/3515638‑3515935

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