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→GGGnanMN‑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

CCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTG  >  NC_000913/4537963‑4538246
                                                                                                                                                    |                                                                                                                                           
cccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCtt                                                                                                                                        <  2:395492/151‑1 (MQ=255)
 ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                       >  1:622110/1‑151 (MQ=255)
 ccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCttt                                                                                                                                       >  1:707495/1‑151 (MQ=255)
   gCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                     >  1:741462/1‑151 (MQ=255)
    cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAt                                                                                                                                     <  2:821204/150‑1 (MQ=255)
      aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTca                                                                                                                                 <  2:738281/151‑1 (MQ=255)
         ccGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACaa                                                                                                                              <  2:419608/151‑1 (MQ=255)
          cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                             >  2:815120/1‑151 (MQ=255)
          cGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAc                                                                                                                             <  2:248529/151‑1 (MQ=255)
             cAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgc                                                                                                                          >  2:590024/1‑151 (MQ=255)
               gCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAAccgccg                                                                                                                        <  2:213916/151‑1 (MQ=255)
                      cATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGc                                                                                                                 >  2:312978/1‑151 (MQ=255)
                       aTCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGcc                                                                                                                >  1:648999/1‑151 (MQ=255)
                          tGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGtt                                                                                                             >  2:347891/1‑151 (MQ=255)
                                   gaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                      >  2:809043/1‑149 (MQ=255)
                                      gggAGCAAGTTTTTTCCTTTTAAAATTATTCCGGGTAAAATCAATTTAAAATCCGGCATCCTTTCGCAACCCTGTTTTGGCTTCGCCTTTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCa                                                                                                     <  1:292657/147‑1 (MQ=255)
                                      gggACCAACTTTTTTCCTTTTTAAATTTTTCCCGGTAAAATAAATTTAAAATCCGGCACCCTTTCCCAATCCTGTTTCGGCTCCCCCATTAATACCCCGGTTTTTATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                   <  2:87942/149‑1 (MQ=255)
                                      gggAACAAGTTTATTCATTTTAAAATTATCCCCGGTAAAATAAAGTTAAAATCCGGCCTCCGTCCGCAATCGGGTTTGGGCTTCCCCTTAAAAAGGCCAGTTTTTATCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCa                                                                                                     <  2:698595/147‑1 (MQ=255)
                                           cAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTc                                                                                            >  2:930215/1‑151 (MQ=255)
                                               cTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                         >  1:93057/1‑150 (MQ=255)
                                                tttttctttttaAAATTATTCCCGGTGAAATCAAGTAAAAATACGCCTCCCTTCCGCATTCCTTGTTTGCCTTCCCCATAAAAAACCCTTTTTTTTCCCCCTTTTTCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                             <  2:705517/137‑1 (MQ=255)
                                                 tATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                      >  2:744003/1‑151 (MQ=255)
                                                   ttCTATTTTAATTTATTCCGGGTAAAATAAATTAAAAATCCGGCACCCGTCCGCAATCTTGGTTGGGCTTCCCCTTTAAAAACCCGGTTTTTATCTCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg                                                                                       <  2:420466/148‑1 (MQ=255)
                                                    tCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGt                                                                                    >  2:857568/1‑150 (MQ=255)
                                                       aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                >  1:883838/1‑151 (MQ=255)
                                                            aaaTTATTCCCGGTAAAATCAAGTTAAAATACGGCATCCGTTCGCAATCCTGTTTGGGCTCCGCCATTAATAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCa                                                                             <  1:519456/149‑1 (MQ=255)
                                                              attattACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc                                                                         >  1:818492/1‑151 (MQ=255)
                                                               ttattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                        >  2:123881/1‑151 (MQ=255)
                                                                   taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                    >  2:273005/1‑151 (MQ=255)
                                                                     ccGGTAAAATAAATTCAAATGCCGCCTTCCGTCCGCAACCTGGTTTTGGCTCCGCTATAAATAGGCCGGGTTTTACCTCTTTTTTTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg                                                                  <  2:126982/151‑1 (MQ=255)
                                                                             aTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc                                                           >  2:844435/1‑150 (MQ=255)
                                                                                aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                        >  2:840408/1‑150 (MQ=255)
                                                                                aaGTTAAAATACGCCCCCGGTTCCCAATCCTGGTTTGTCTCCCCCATTAATAAGCCGGTTTTTATCCCTTTTTTCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                        <  1:960999/150‑1 (MQ=255)
                                                                                 aGTTAAAATCCGCCATCCGTTCGCAATCCTGTTTTGGCTCCGCCATAAAAAAGCCAGTTTTAATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAg                                                      <  1:918813/151‑1 (MQ=255)
                                                                                      aaaTCCGCCATCCGTTCCCAACCCTGTTTTGCCTTCGCCATTAATAAGCCAGTTTTTATCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa                                                     <  1:42102/147‑1 (MQ=255)
                                                                                      aaaTACGGCATCCTTTCGCAACCCTGTTTTGGTTCCGCCATTAAAAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                  <  1:77070/150‑1 (MQ=255)
                                                                                           aCGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                              >  2:55466/1‑150 (MQ=255)
                                                                                            cGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                             >  1:315816/1‑150 (MQ=255)
                                                                                               cATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                         >  2:558610/1‑151 (MQ=255)
                                                                                                  ccGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                      >  1:940108/1‑151 (MQ=255)
                                                                                                       cgaaaTCCTGGTTGGCCTTCGCCATAAAAAAGCCAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                    <  2:426460/145‑1 (MQ=255)
                                                                                                       caaaacccgggTTTGGCTCCGCCTTTAAAAAGCCGGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                 <  1:913265/142‑1 (MQ=255)
                                                                                                       cGCAATCCTGTTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                 <  1:796527/151‑1 (MQ=255)
                                                                                                        ccgaTCTTGGTTTGCCTTCCCCATTAAAAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                  <  1:4357/146‑1 (MQ=255)
                                                                                                          atcccTGTTTTGGCTTCCCAATAAAAAAGCCAGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                               <  1:620641/147‑1 (MQ=255)
                                                                                                          aaTCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                              >  2:475932/1‑151 (MQ=255)
                                                                                                            tCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                             >  1:473291/1‑150 (MQ=255)
                                                                                                               tGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaag                                  >  2:4356/1‑142 (MQ=255)
                                                                                                                 gttttGGCTCCGCAATTAAAAAGCCAGTTTTTATCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                           <  1:660772/146‑1 (MQ=255)
                                                                                                                 gTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                       >  1:379968/1‑151 (MQ=255)
                                                                                                                  ttttGGCTTCCCCATAAATAAGCCGGGTTTTATCCCCTTTATTCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                        <  1:193599/149‑1 (MQ=255)
                                                                                                                   tttGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                      >  2:547106/1‑151 (MQ=255)
                                                                                                                   ttgggttcccccATAAAAAACCCGGTTTTTACCCCCTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                       <  1:575968/140‑1 (MQ=255)
                                                                                                                           ccctataaaaagcccAGGTTTTACCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttttt              <  2:93039/138‑1 (MQ=255)
                                                                                                                           ccccATTAAAAACCCGGGTTTTACCCCTTTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt               <  1:767189/148‑1 (MQ=255)
                                                                                                                             ctaTTAATAACCCAGTTTTTATCCCTTTTACCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg             <  2:540461/148‑1 (MQ=255)
                                                                                                                             cctataaaaaGCCTGTTTTTTCCCCCTTTTTTCCACACCCCCCGCCCCAGCCGTTCCGTCCCAGGGCGATTCGCCAGCTTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAACCTTATTGAAGAAATAATCTTCTGCttttt               <  2:508543/145‑1 (MQ=255)
                                                                                                                              cATTAAAGAACCGGGTTTTTTCCCCTTTTATCCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtttt                <  2:98675/147‑1 (MQ=255)
                                                                                                                               aTTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc           >  2:83380/1‑150 (MQ=255)
                                                                                                                                  aaaaacccAGTTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa         <  1:840621/143‑1 (MQ=255)
                                                                                                                                     aaGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa     <  1:857782/150‑1 (MQ=255)
                                                                                                                                     aaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt      >  2:593811/1‑149 (MQ=255)
                                                                                                                                       gCCGGGTTTAACCCCTTTTACCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt   <  1:92932/150‑1 (MQ=255)
                                                                                                                                        ccAGGTTTTACCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt   <  1:383399/149‑1 (MQ=255)
                                                                                                                                         cAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg  >  2:839449/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTG  >  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: