Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,1090AC41.3% 115.0 / 14.8 65G219G (GGT→GGG) nanMN‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (31/6);  new base C (0/26);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 3.01e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTG  >  NC_000913/4537963‑4538246
                                                                                                                                                    |                                                                                                                                           
cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑tt                                                                                                                                        <  2:395492/151‑1 (MQ=255)
 ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑ttt                                                                                                                                       >  1:622110/1‑151 (MQ=255)
 ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑ttt                                                                                                                                       >  1:707495/1‑151 (MQ=255)
   gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAt                                                                                                                                     >  1:741462/1‑151 (MQ=255)
    cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTAt                                                                                                                                     <  2:821204/150‑1 (MQ=255)
      aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑ca                                                                                                                                 <  2:738281/151‑1 (MQ=255)
         ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACaa                                                                                                                              <  2:419608/151‑1 (MQ=255)
          cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAAc                                                                                                                             >  2:815120/1‑151 (MQ=255)
          cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAAc                                                                                                                             <  2:248529/151‑1 (MQ=255)
             cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAAccgc                                                                                                                          >  2:590024/1‑151 (MQ=255)
               gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAAccgccg                                                                                                                        <  2:213916/151‑1 (MQ=255)
                      cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGc                                                                                                                 >  2:312978/1‑151 (MQ=255)
                       aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGcc                                                                                                                >  1:648999/1‑151 (MQ=255)
                          tGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGtt                                                                                                             >  2:347891/1‑151 (MQ=255)
                                   gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                      >  2:809043/1‑149 (MQ=255)
                                      gggAGCAAGTTTTTTCCTTTTAAAATTATTCCGGGTAAAATCAATTTAAAAT‑CCGGCATCCTTTCGCAACCCTGT‑TTTGGCTTCGCCTTTAATAAGCCAGTTTTTATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCa                                                                                                     <  1:292657/147‑1 (MQ=255)
                                      gggACCAACTTTTTTCCTTTTTAAATTTTTCCCGGTAAAATAAATTTAAAAT‑CCGGCACCCTTTCCCAATCCTGT‑TTCGGCTCCCCCATTAATACCCCGGTTTTTATCCCT‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                   <  2:87942/149‑1 (MQ=255)
                                      gggAACAAGTTTATTCATTTTAAAATTATCCCCGGTAAAATAAAGTTAAAAT‑CCGGCCTCCGTCCGCAATCGGGT‑TTGGGCTTCCCCTTAAAAAGGCCAGTTTTTATCCCT‑TTTATC‑CAAACCCGCCGCACCAGCCGTTCCGTACCa                                                                                                     <  2:698595/147‑1 (MQ=255)
                                           cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                            >  2:930215/1‑151 (MQ=255)
                                               cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                         >  1:93057/1‑150 (MQ=255)
                                                tttttctttttaAAATTATTCCCGGTGAAATCAAGTAAAAAT‑ACGCCTCCCTTCCGCATTCCTTG‑TTTGCCTTCCCCATAAAAAACCCTTTTTTTTCCCCC‑TTTTTC‑CCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                             <  2:705517/137‑1 (MQ=255)
                                                 tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                      >  2:744003/1‑151 (MQ=255)
                                                   ttCTATTTTAATTTATTCCGGGTAAAATAAATTAAAAAT‑CCGGCACCCGTCCGCAATCTTGG‑TTGGGCTTCCCCTTTAAAAACCCGGTTTTTATCTCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                       <  2:420466/148‑1 (MQ=255)
                                                    tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                    >  2:857568/1‑150 (MQ=255)
                                                       aTTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                >  1:883838/1‑151 (MQ=255)
                                                            aaaTTATTCCCGGTAAAATCAAGTTAAAAT‑ACGGCATCCGTTCGCAATCCTGT‑TTGGGCTCCGCCATTAATAAGCCAGTTTTTATCCCT‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                             <  1:519456/149‑1 (MQ=255)
                                                              attattACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                         >  1:818492/1‑151 (MQ=255)
                                                               ttattaCCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                        >  2:123881/1‑151 (MQ=255)
                                                                   taCCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                    >  2:273005/1‑151 (MQ=255)
                                                                     ccGGTAAAAT‑AAATTCAAATGCCGCCTTCCGTCCGCAA‑CCTGGTTTTGGCTCCGCTATAAATAGGCCGGGTTTTACCTCT‑TTTTTT‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                  <  2:126982/151‑1 (MQ=255)
                                                                             aTCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                           >  2:844435/1‑150 (MQ=255)
                                                                                aaGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                        >  2:840408/1‑150 (MQ=255)
                                                                                aaGTTAAAAT‑ACGCCCCCGGTTCCCAATCCTGG‑TTTGTCTCCCCCATTAATAAGCCGGTTTTTATCCCT‑TTTTTC‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                        <  1:960999/150‑1 (MQ=255)
                                                                                 aGTTAAAAT‑CCGCCATCCGTTCGCAATCCTGT‑TTTGGCTCCGCCATAAAAAAGCCAGTTTTAATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                      <  1:918813/151‑1 (MQ=255)
                                                                                      aaaT‑CCGCCATCCGTTCCCAACCCTGT‑TTTGCCTTCGCCATTAATAAGCCAGTTTTTATCCCT‑TTTATC‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                     <  1:42102/147‑1 (MQ=255)
                                                                                      aaaT‑ACGGCATCCTTTCGCAACCCTGT‑TTTGGTTCCGCCATTAAAAAGCCAGTTTTTATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                  <  1:77070/150‑1 (MQ=255)
                                                                                           aCGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                              >  2:55466/1‑150 (MQ=255)
                                                                                            cGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                             >  1:315816/1‑150 (MQ=255)
                                                                                               cATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                         >  2:558610/1‑151 (MQ=255)
                                                                                                  ccGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                      >  1:940108/1‑151 (MQ=255)
                                                                                                       cgaaaTCCTGG‑TTGGCCTTCGCCATAAAAAAGCCAGTTTTTACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                    <  2:426460/145‑1 (MQ=255)
                                                                                                       caaaacccggg‑TTTGGCTCCGCCTTTAAAAAGCCGGTTTTTATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                 <  1:913265/142‑1 (MQ=255)
                                                                                                       cGCAATCCTGT‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                 <  1:796527/151‑1 (MQ=255)
                                                                                                        ccgaTCTTGG‑TTTGCCTTCCCCATTAAAAAGCCAGTTTTTATCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                  <  1:4357/146‑1 (MQ=255)
                                                                                                          atcccTGT‑TTTGGCTTCCCAATAAAAAAGCCAGGTTTTACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                               <  1:620641/147‑1 (MQ=255)
                                                                                                          aaTCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                              >  2:475932/1‑151 (MQ=255)
                                                                                                            tCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                             >  1:473291/1‑150 (MQ=255)
                                                                                                               tGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                  >  2:4356/1‑142 (MQ=255)
                                                                                                                 gttttGGCTCCGCAATTAAAAAGCCAGTTTTTATCCCT‑TTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                           <  1:660772/146‑1 (MQ=255)
                                                                                                                 g‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                       >  1:379968/1‑151 (MQ=255)
                                                                                                                  ttttGGCTTCCCCATAAATAAGCCGGGTTTTATCCCC‑TTTATT‑CCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                        <  1:193599/149‑1 (MQ=255)
                                                                                                                   tttGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                      >  2:547106/1‑151 (MQ=255)
                                                                                                                   ttgggttcccccATAAAAAACCCGGTTTTTACCCCC‑TTTATC‑CAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                       <  1:575968/140‑1 (MQ=255)
                                                                                                                           ccctataaaaagcccAGGTTTTACCCCT‑TTTATC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt              <  2:93039/138‑1 (MQ=255)
                                                                                                                           ccccATTAAAAACCCGGGTTTTACCCCT‑TTTTTT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt               <  1:767189/148‑1 (MQ=255)
                                                                                                                             ctaTTAATAACCCAGTTTTTATCCCT‑TTTACC‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg             <  2:540461/148‑1 (MQ=255)
                                                                                                                             cctataaaaaGCCTGTTTTTTCCCCC‑TTTTTTCCACACCCCCCGCCCCAGCCGTTCCGTCCCAGGGCGATTCGCCAGCTTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAACCTTATTGAAGAAATAATCTTCTGCttttt               <  2:508543/145‑1 (MQ=255)
                                                                                                                              cATTAAAGAACCGGGTTTTTTCCCCTTTTATC‑CAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                <  2:98675/147‑1 (MQ=255)
                                                                                                                               aTTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc           >  2:83380/1‑150 (MQ=255)
                                                                                                                                  aaaaacccAGTTTTTACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa         <  1:840621/143‑1 (MQ=255)
                                                                                                                                     aaGCCAGTTTTTATCCCC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa     <  1:857782/150‑1 (MQ=255)
                                                                                                                                     aaGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt      >  2:593811/1‑149 (MQ=255)
                                                                                                                                       gCCGGGTTTAACCCCT‑TTTACC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt   <  1:92932/150‑1 (MQ=255)
                                                                                                                                        ccAGGTTTTACCCCC‑TTTATC‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt   <  1:383399/149‑1 (MQ=255)
                                                                                                                                         cAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg  >  2:839449/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAAT‑ACGGCATCCGTTCGCAATCCTGG‑TTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACC‑TTTATT‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTG  >  NC_000913/4537963‑4538246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑