Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,515,786 A→C 5.5% L43R (CTC→CGC) ‡ dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,7860AC 5.5% 418.4 / 10.8 139L43R (CTC→CGC) ‡damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (72/57);  new base C (6/2);  total (78/61)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 149 total aligned reads displayed.

GATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTT  >  NC_000913/3515637‑3515935
                                                                                                                                                     |                                                                                                                                                     
gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                                                                       <  1:1179521/150‑1 (MQ=255)
gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                                                                       <  1:2162194/150‑1 (MQ=255)
gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                                                                       <  2:1252803/150‑1 (MQ=255)
gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                                                                       <  1:2170570/150‑1 (MQ=255)
gATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTga                                                                                                                                                       >  2:1429646/1‑150 (MQ=255)
 aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                                                                      <  1:1995604/150‑1 (MQ=255)
 aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTgag                                                                                                                                                      <  1:593276/150‑1 (MQ=255)
   aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                    <  1:85600/150‑1 (MQ=255)
   aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                    <  1:804552/150‑1 (MQ=255)
   aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                    <  2:1132121/150‑1 (MQ=255)
      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 >  2:1889181/1‑150 (MQ=255)
      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 <  1:2023118/150‑1 (MQ=255)
      aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 <  2:2037325/150‑1 (MQ=255)
        aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACAcc                                                                                                                                               >  1:1398524/1‑150 (MQ=255)
          cTCGGCTCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGa                                                                                                                                             <  2:1669709/150‑1 (MQ=255)
          cTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGa                                                                                                                                             >  1:625450/1‑150 (MQ=255)
             ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                                                                          <  1:1720312/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)
                         ttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaaa                                                                                                                              >  2:888447/1‑150 (MQ=255)
                          tCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGCAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaaaa                                                                                                                             >  2:128828/1‑150 (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)
                             cgAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTAATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACTAGAAAAGTCGGTGTAGAGAAACACAGACCCGGCACCTACAAAAGGc                                                                                                                          >  1:179687/2‑150 (MQ=255)
                              gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACAGACCCGGCCCCTACAAAAGGCt                                                                                                                         >  1:2161530/1‑150 (MQ=255)
                                  aaacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAc                                                                                                                     <  1:239979/150‑1 (MQ=255)
                                   aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                                                                    >  2:892819/1‑150 (MQ=255)
                                   aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                                                                    >  1:1242668/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)
                                            cgTGCGGCCTGGACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                           <  2:173217/150‑1 (MQ=255)
                                              tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATAAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCAACAAAAGGCTCAACCAGACATTCGcc                                                                                                         >  1:690973/1‑150 (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                                                                                        >  2:1880807/1‑150 (MQ=255)
                                                                   aGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAACGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTATAat                                                                                    >  1:1336910/1‑150 (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)
                                                                            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)
                                                                                  aaTGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgggg                                                                     <  2:377858/150‑1 (MQ=255)
                                                                                      ttATAGAGCCTGATCAGGTCGCTATTGATATCGGCAAGGATAAAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATAc                                                                 <  1:1772816/150‑1 (MQ=255)
                                                                                       tataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACt                                                                >  2:1356403/1‑150 (MQ=255)
                                                                                           gagaCTGATCAGGGCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGcc                                                            <  2:1322374/150‑1 (MQ=255)
                                                                                            agaCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGccc                                                           <  1:264843/150‑1 (MQ=255)
                                                                                             gaCTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAAGATCATCAAGCAGGGGATACTTGcccc                                                          >  2:435704/1‑150 (MQ=255)
                                                                                                 gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc                                                      >  1:860842/1‑150 (MQ=255)
                                                                                                 gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc                                                      >  1:818116/1‑150 (MQ=255)
                                                                                                  aTCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAACGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGcc                                                     >  2:61629/1‑150 (MQ=255)
                                                                                                     aGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc                                                  <  2:283133/150‑1 (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)
                                                                                                             aTTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAATTGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaaa                                          <  2:2210741/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:1136618/150‑1 (MQ=255)
                                                                                                               tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                        <  1:892819/150‑1 (MQ=255)
                                                                                                                gATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                       >  2:1992466/1‑150 (MQ=255)
                                                                                                                gATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACAAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                       >  1:506585/1‑150 (MQ=255)
                                                                                                                 atatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgcg                                      >  2:390939/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)
                                                                                                                       gCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAttt                                >  1:973158/1‑150 (MQ=255)
                                                                                                                           ggATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttct                            >  1:1267164/1‑150 (MQ=255)
                                                                                                                           ggATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttct                            <  2:1171992/150‑1 (MQ=255)
                                                                                                                            gttattACGAGAAAAGTCGGTGTTGAGAAACTCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCATTTGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttctt                           <  2:179687/148‑1 (MQ=255)
                                                                                                                                taACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATg                       >  1:269782/1‑150 (MQ=255)
                                                                                                                                  aCGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTACAAAAAAGCGCGATTTTTCTTCATGCt                     >  1:1368042/1‑150 (MQ=255)
                                                                                                                                   ctagaAAAGTCGGTGTTTAGAAACACCGACCCTGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTg                    <  2:1397197/148‑1 (MQ=255)
                                                                                                                                   cGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTg                    <  1:1047450/150‑1 (MQ=255)
                                                                                                                                     agaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGac                  >  2:1159616/1‑150 (MQ=255)
                                                                                                                                     agaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGac                  <  2:1014873/150‑1 (MQ=255)
                                                                                                                                       ccgggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGacta                >  1:897593/5‑150 (MQ=255)
                                                                                                                                          ggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             >  1:223176/2‑150 (MQ=255)
                                                                                                                                          ggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             >  1:1469468/2‑150 (MQ=255)
                                                                                                                                          aggcGGTGTTGATAAACACGGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGAAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             <  1:1642405/147‑1 (MQ=255)
                                                                                                                                              ggTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAAtt         <  2:634859/150‑1 (MQ=255)
                                                                                                                                               gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa        >  1:654080/1‑150 (MQ=255)
                                                                                                                                               gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa        <  1:966618/150‑1 (MQ=255)
                                                                                                                                                  ttGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTacac     <  2:101035/150‑1 (MQ=255)
                                                                                                                                                    gagaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCt   <  2:1155226/150‑1 (MQ=255)
                                                                                                                                                     agaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCtt  <  2:1625587/150‑1 (MQ=255)
                                                                                                                                                     agaAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCACTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCtt  >  2:2138969/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTT  >  NC_000913/3515637‑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: