Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,636,877 T→C V41A (GTC→GCC)  ynfO → uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,8770TC98.5% 239.9 / ‑3.9 67V41A (GTC→GCC) ynfOuncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (36/30);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

GAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCCCCTC  >  NC_000913/1636739‑1637024
                                                                                                                                          |                                                                                                                                                   
gAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCt                                                                                                                                          <  1:456430/150‑1 (MQ=14)
  aTACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTgct                                                                                                                                       <  1:1207178/151‑2 (MQ=14)
    aCGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTgccgc                                                                                                                                     <  1:188077/151‑1 (MQ=17)
     cGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTgccgcc                                                                                                                                    >  2:1176847/1‑151 (MQ=18)
        ggCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGga                                                                                                                                 <  1:218445/151‑1 (MQ=18)
             ccTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACtg                                                                                                                            <  2:51197/151‑1 (MQ=37)
               tCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACtgtg                                                                                                                          >  2:1093863/1‑151 (MQ=38)
                cATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGc                                                                                                                         >  1:1163130/1‑151 (MQ=39)
                   gTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCa                                                                                                                      <  2:637255/151‑1 (MQ=255)
                     cAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTc                                                                                                                       >  2:135822/1‑148 (MQ=255)
                      aaCGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAg                                                                                                                     >  2:845646/1‑149 (MQ=255)
                        cGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGaaa                                                                                                                  >  1:431167/1‑150 (MQ=255)
                         gaagaaCAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGCTTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGaaaa                                                                                                                 >  2:1124986/1‑150 (MQ=255)
                           agaacagaacCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAag                                                                                                              <  2:610669/151‑1 (MQ=255)
                              acagaacCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGt                                                                                                           >  1:877316/1‑151 (MQ=255)
                                agaacCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGttt                                                                                                         <  2:1025769/151‑1 (MQ=255)
                                       gCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTg                                                                                                   >  1:915621/1‑150 (MQ=255)
                                       gCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGa                                                                                                  <  2:477041/151‑1 (MQ=255)
                                         aGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGa                                                                                                  >  1:494473/1‑149 (MQ=255)
                                         aGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGa                                                                                                  >  2:769684/1‑149 (MQ=255)
                                              aacaacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaac                                                                                           >  1:1221161/1‑151 (MQ=255)
                                                 aacCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaaca                                                                                          >  1:16153/1‑149 (MQ=255)
                                                  acCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGcaaacaa                                                                                         >  1:315021/1‑149 (MQ=255)
                                                     cGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaa                                                                                     >  1:888712/1‑150 (MQ=255)
                                                      gCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaaa                                                                                   <  2:268469/151‑1 (MQ=255)
                                                       cAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaaag                                                                                  <  1:173760/151‑1 (MQ=255)
                                                        aaCATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACaaagaaa                                                                                    >  1:364936/1‑148 (MQ=255)
                                                             ccGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGa                                                                             <  1:1075707/150‑1 (MQ=255)
                                                                cTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCACCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACAtt                                                                         >  2:780525/1‑151 (MQ=255)
                                                                 tttCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACAtt                                                                         >  2:221065/1‑150 (MQ=255)
                                                                    ccTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTg                                                                     <  2:653019/151‑1 (MQ=255)
                                                                     cTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTg                                                                     >  1:753588/1‑150 (MQ=255)
                                                                      tAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGa                                                                    >  1:932064/1‑150 (MQ=255)
                                                                       aaCCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGaa                                                                   <  1:566712/150‑1 (MQ=255)
                                                                       aaCCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAAt                                                                  <  2:1091223/151‑1 (MQ=255)
                                                                        aCCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGaa                                                                   <  2:412439/149‑1 (MQ=255)
                                                                         ccAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTg                                                                <  1:882/151‑1 (MQ=255)
                                                                         ccAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTg                                                                <  2:877311/151‑1 (MQ=255)
                                                                         ccAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTg                                                                >  1:917154/1‑151 (MQ=255)
                                                                             aTGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTc                                                             >  2:796049/1‑150 (MQ=255)
                                                                             aTGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCa                                                            <  2:324053/151‑1 (MQ=255)
                                                                                 ttGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGaa                                                         >  1:552084/1‑150 (MQ=255)
                                                                                 ttGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGaa                                                         >  2:1009018/1‑150 (MQ=255)
                                                                                  tGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGaaa                                                        >  1:658675/1‑150 (MQ=255)
                                                                                    aaCAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACa                                                     <  1:1036068/151‑1 (MQ=255)
                                                                                      cAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAgg                                                   >  1:587402/1‑151 (MQ=255)
                                                                                       aaaTTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGt                                                  >  1:1217895/1‑151 (MQ=255)
                                                                                           tAACATCGCTCTAGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGATTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAg                                               <  2:221667/150‑1 (MQ=255)
                                                                                            aaCATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGTTCAGAAAACAGGTCAGAc                                             >  1:548126/1‑151 (MQ=255)
                                                                                              cATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGAcc                                            >  1:304405/1‑150 (MQ=255)
                                                                                                  gCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGAtt                                       <  1:569098/151‑1 (MQ=255)
                                                                                                  gCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGAtt                                       >  1:376762/1‑151 (MQ=255)
                                                                                                       tGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAg                                   <  1:1148377/150‑1 (MQ=255)
                                                                                                        gAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGc                                  >  1:306398/1‑150 (MQ=255)
                                                                                                               aaGGGTCCGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGa                          <  2:990923/151‑1 (MQ=255)
                                                                                                                     ccGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATc                     <  1:473136/150‑1 (MQ=255)
                                                                                                                     ccGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCg                    <  1:135815/151‑1 (MQ=255)
                                                                                                                      cGGGAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGt                   <  1:1020832/151‑1 (MQ=255)
                                                                                                                       gggAATTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTg                  >  1:314503/1‑151 (MQ=255)
                                                                                                                           aTTTCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGa              <  2:381151/151‑1 (MQ=255)
                                                                                                                            tttCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGAc             >  1:13165/1‑151 (MQ=255)
                                                                                                                             ttCTCAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACt            >  1:769646/1‑151 (MQ=255)
                                                                                                                                 cAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTAc          <  2:970934/149‑1 (MQ=255)
                                                                                                                                 cAGCCTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACg         <  1:1045288/150‑1 (MQ=255)
                                                                                                                                     cTGGGCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGccccc    >  1:517585/1‑151 (MQ=255)
                                                                                                                                       gggCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCCCCt   >  1:1032023/1‑150 (MQ=255)
                                                                                                                                       gggCCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCCCCTc  >  1:110579/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAACATCGCTCTTGAGCAAAAAGGGTCCGGGAATTTCTCAGCCTGGGTCATTGAAGCCTGCCGCCGGAGACTGTGCTCAGAAAAAAGAGTTTCTTCTGAAGCAAACAAAGAAAAGAGTGACATTACTGAATTGCTCAGAAAACAGGTCAGACCAGATTGAAGCAATTTAGATAATCGTGCAGACTACGCCCCCTC  >  NC_000913/1636739‑1637024

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: