Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,515,788 A→C 5.9% F42L (TTT→TTG dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,7880AC 5.9% 403.4 / 12.3 135F42L (TTT→TTGdamDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (71/56);  new base C (6/2);  total (77/58)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 142 total aligned reads displayed.

AGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTT  >  NC_000913/3515640‑3515935
                                                                                                                                                    |                                                                                                                                                   
aGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                                                                                                                                    <  2:1132121/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                                                                                                                                                 <  2:1181597/150‑1 (MQ=255)
   aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 <  1:2023118/150‑1 (MQ=255)
   aaaCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAaca                                                                                                                                                 <  2:2037325/150‑1 (MQ=255)
     aCCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACAcc                                                                                                                                               >  1:1398524/1‑150 (MQ=255)
       cTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTTCTCAACACCGa                                                                                                                                             <  1:714211/150‑1 (MQ=255)
       cTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGa                                                                                                                                             >  1:625450/1‑150 (MQ=255)
          ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGAccc                                                                                                                                          <  1:1720312/150‑1 (MQ=255)
                   tGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTAc                                                                                                                                 >  2:2053180/1‑150 (MQ=255)
                      ttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaaa                                                                                                                              >  2:888447/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)
                         gggAACACACAGCCCGCGGGCGGCCTTTACGTAATAATCCGTACGCATATTCACATTGTTATAGAGACTGATCCGGTCGCTATTGAAATCGGCATGGATATAACGACAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                                                                                                                           <  2:690824/150‑1 (MQ=255)
                          ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCTTCTTCACATTTTTATAGAGACTTATCAGGTCGCTATTTATATCGGCAAGTATATAACGATAAAAGTCGGTGTTGATAAACACCTACCCGGCACCTACAAAAGGc                                                                                                                          <  1:2044803/150‑1 (MQ=255)
                          ggAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTGTAGAGACTGATCAGGTCGCTATTTATATCGGCAAGGATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                          <  2:2165441/150‑1 (MQ=255)
                          ggAACAAACAGCTCGCGTGCGGCCTGAACGTACTCATCAGTACGCATCTTCACAATTTAATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGTATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGc                                                                                                                          <  2:1850661/150‑1 (MQ=255)
                          cgAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTAATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACTAGAAAAGTCGGTGTAGAGAAACACAGACCCGGCACCTACAAAAGGc                                                                                                                          >  1:179687/2‑150 (MQ=255)
                           gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACAGACCCGGCCCCTACAAAAGGCt                                                                                                                         >  1:2161530/1‑150 (MQ=255)
                           gAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGAACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCt                                                                                                                         <  2:38121/150‑1 (MQ=255)
                                aacaGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc                                                                                                                    >  1:1242668/1‑150 (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:246057/150‑1 (MQ=255)
                                    gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGAc                                                                                                                >  1:1171992/1‑150 (MQ=255)
                                    gCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGAc                                                                                                                >  1:20673/1‑150 (MQ=255)
                                         cgTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTCTAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCg                                                                                                           <  2:1162349/150‑1 (MQ=255)
                                           tGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATAAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCAACAAAAGGCTCAACCAGACATTCGcc                                                                                                         >  1:690973/1‑150 (MQ=255)
                                              ggCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCtt                                                                                                      >  2:1948340/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)
                                                    tacgtacTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGTATATAACGAGAAAAGTCGGTGTTGAGAAACGCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCtgaggcta                                                                                                >  1:318382/1‑143 (MQ=255)
                                                            catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                        >  1:2168576/1‑150 (MQ=255)
                                                            catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                        >  1:2126104/1‑150 (MQ=255)
                                                            catcaGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGtt                                                                                        >  2:1880807/1‑150 (MQ=255)
                                                               cattACGAATCTTAACAATGTTATAGAGAATTATCCGGGCTCTATTTATATCGGCAAGGATATAACGATAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTaa                                                                                     <  1:1900736/147‑1 (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)
                                                                      taTCTTAACAATGTTATAGATACTGTTCAGGTCGCTATTGATATCGGCAAGTATATAACGAGATAAGTCGGTGTTGTTAAACCCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATatcatc                                                                              <  1:1889181/149‑1 (MQ=255)
                                                                       aTCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATAtcatca                                                                             <  1:2001899/150‑1 (MQ=255)
                                                                         cTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAg                                                                           >  2:1717979/1‑150 (MQ=255)
                                                                          ttCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGc                                                                          >  2:803251/1‑150 (MQ=255)
                                                                            cacaATGTTATAGAGACTGATCAGGTCGCTATTTATATCGGCAAGGATATAACGAGAAAAGTCGGAGATGAGAAACACCGACCCGCCACCAACAAAAGGCTCAACCACACAAACGCCCATCGGCAAATGCCGTAAAATATCATCAAGCac                                                                        >  2:857850/1‑149 (MQ=255)
                                                                              caATGTTATATAGACTGATAAGGTCGCTATTGATATGGGCAAGGATATATCGAGAAAAGTCGGTTTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCTCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAggg                                                                      <  2:298478/150‑1 (MQ=255)
                                                                               aaTGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAgggg                                                                     <  2:377858/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)
                                                                                     ataGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACtt                                                               >  1:1207732/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:989177/1‑150 (MQ=255)
                                                                                              gATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGc                                                      >  1:860842/1‑150 (MQ=255)
                                                                                                  aGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCAc                                                  <  2:283133/150‑1 (MQ=255)
                                                                                                   ggTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCGCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACt                                                 >  1:61294/1‑150 (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)
                                                                                                           ttGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCaaaaaa                                         >  2:713099/1‑150 (MQ=255)
                                                                                                            tGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACTCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAg                                        >  1:1009897/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)
                                                                                                             gATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACAAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgc                                       >  1:506585/1‑150 (MQ=255)
                                                                                                              atatCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAgcg                                      >  2:390939/1‑150 (MQ=255)
                                                                                                                atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                    <  2:626459/150‑1 (MQ=255)
                                                                                                                atCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATAATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGcgcg                                    >  2:1688154/1‑150 (MQ=255)
                                                                                                                   ggCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGGTTAATATCATCAAGCAGGGGATACTTGCCCCGTGCCCACTTCAAAAAAGCGCGAtt                                 >  1:846138/1‑150 (MQ=255)
                                                                                                                    gCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAttt                                >  1:973158/1‑150 (MQ=255)
                                                                                                                    gCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGAttt                                <  2:625450/150‑1 (MQ=255)
                                                                                                                         gttattACGAGAAAAGTCGGTGTTGAGAAACTCCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCATTTGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTttctt                           <  2:179687/148‑1 (MQ=255)
                                                                                                                              aaCGAGAAATGTCGGTGTTGAGAAACACCTACCCGGCACCTACAATAGTCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGc                      <  2:1380779/150‑1 (MQ=255)
                                                                                                                               aCGAGAAAAGTCGTTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCt                     <  1:1534690/150‑1 (MQ=255)
                                                                                                                                ctagaAAAGTCGGTGTTTAGAAACACCGACCCTGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTg                    <  2:1397197/148‑1 (MQ=255)
                                                                                                                                cGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTg                    <  1:1047450/150‑1 (MQ=255)
                                                                                                                                cGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGAGCGATTTTTCTTCATGCTg                    >  2:1447400/1‑150 (MQ=255)
                                                                                                                                  agaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGac                  <  2:1014873/150‑1 (MQ=255)
                                                                                                                                  agaAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGac                  >  2:1159616/1‑150 (MQ=255)
                                                                                                                                    ccgggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGacta                >  1:897593/5‑150 (MQ=255)
                                                                                                                                     aaaGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaa               >  1:2205135/1‑150 (MQ=255)
                                                                                                                                      gggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaac              >  2:1697622/3‑150 (MQ=255)
                                                                                                                                       ggTCGGTGTTTCTCAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             >  1:223176/2‑150 (MQ=255)
                                                                                                                                       aggcGGTGTTGATAAACACGGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGAAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGactaact             <  1:1642405/147‑1 (MQ=255)
                                                                                                                                            gtgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTa        <  1:966618/150‑1 (MQ=255)
                                                                                                                                             tgtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTac       >  2:1296302/1‑150 (MQ=255)
                                                                                                                                              gtTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTaca      <  2:1739314/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)
                                                                                                                                                    |                                                                                                                                                   
AGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATTCGCCCTTGGGCAAATGCCGTTTAATATCATCAAGCAGGGGATACTTGCCCCCTGCCCACTTCAAAAAAGCGCGATTTTTCTTCATGCTGACTAACTAATTACACCTT  >  NC_000913/3515640‑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: