Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC40.8% 198.7 / 25.6 98V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (57/1);  new base C (0/40);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 8.34e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378653‑4378952
                                                                                                                                                             |                                                                                                                                                      
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ac                                                                                                                                                       <  4:135002/151‑1 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ac                                                                                                                                                       >  3:301508/1‑151 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑a                                                                                                                                                        <  2:291757/150‑1 (MQ=255)
 cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑acc                                                                                                                                                      >  2:211882/1‑151 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑accg                                                                                                                                                     <  1:39499/151‑1 (MQ=255)
       aGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑accgccacc                                                                                                                                                >  3:751763/1‑151 (MQ=255)
          aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑accgccaccgcc                                                                                                                                             <  2:638049/151‑1 (MQ=255)
               cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATAcc                                                                                                                                        >  2:265698/1‑151 (MQ=255)
               cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATAcc                                                                                                                                        >  3:1233708/1‑151 (MQ=255)
                 tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCg                                                                                                                                      >  3:1402502/1‑151 (MQ=255)
                     caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGAt                                                                                                                                  >  2:4304/1‑151 (MQ=255)
                     caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGAt                                                                                                                                  >  2:459075/1‑151 (MQ=255)
                      accaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATc                                                                                                                                 >  1:33478/1‑151 (MQ=255)
                         aggcccGCCAGAAAATTTTTTGCTGCCCGAACCAAACTCAAACATCTTTTGCTGTGTGGCGGGCTTTTTTTTGGCGT‑TGTCCGCATACCCCCGGGA‑AAATAAAAAAGGTTT‑CCCG‑TG‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATc                                                                                                                                 <  1:237740/143‑1 (MQ=255)
                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCttt                                                                                                                              >  3:926569/1‑151 (MQ=255)
                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCttt                                                                                                                              >  2:1163178/1‑151 (MQ=255)
                            cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                           >  2:418657/1‑151 (MQ=255)
                                ccAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                       >  2:1207859/1‑151 (MQ=255)
                                     aatttttttGCGGCCCGAACCAAACTAAAACAGGTTTTGTTGTGTG‑GGGGGTTTTTTTTGGGGATTGCCCGAATACCCCAAGAA‑AAAGAAAAAGGTTTT‑CCCC‑TA‑AAAA‑T‑TT‑CCCCCCACCGCCATACCCGGGATCtttt                                                                                                                             <  3:285849/135‑1 (MQ=255)
                                       aTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  >  4:255507/1‑149 (MQ=255)
                                           tttGCTGACAAACCCTACCTAAAACAACTTTTGTTGTGTGGCGGGTTGTTTTTGGGCGA‑GGCCCGCATACCCCCAGGA‑AAAG‑AAAAAGGTTTCCCCC‑TT‑AAAA‑A‑TC‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                             <  3:325191/150‑1 (MQ=255)
                                            ttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                           >  2:754923/1‑151 (MQ=255)
                                               cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGTCCTTTTGTTGCTTTTTAAGCggggaa                                                                                                         >  3:766518/1‑148 (MQ=255)
                                               cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         >  3:253080/1‑150 (MQ=255)
                                               cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         >  2:655662/1‑150 (MQ=255)
                                               cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         >  4:113726/1‑150 (MQ=255)
                                                 gcccacaCCTAACCAAACCAGCGTTTGTGGTGTGCCGGCTTGTTTTTGGGCGA‑TGCCCGAATGGCCCCGGAA‑AAAGTAATAGGTTTT‑TCCT‑TT‑AAAA‑C‑TA‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                       <  3:935615/144‑1 (MQ=255)
                                                 gacagacccaaaccaaaaCAGCTTTTGCTGTGTGAGGGGCTTTTTTTTGGGGT‑TGACC‑CAACACCCCATGTAAAAATAAAAAGTTTT‑CCCT‑AT‑ATAA‑T‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                       <  2:831784/135‑1 (MQ=255)
                                                    cGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                   >  1:376036/1‑151 (MQ=255)
                                                         cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                              >  2:60261/1‑151 (MQ=255)
                                                           aaCACAAACACCGTTTTCTGTGGGAGGGTCTGTTTTTTGGGGA‑TGCCCGAATAGCCCGAGGA‑AAAGTATTAAGTTTT‑CCCT‑TA‑ATAA‑T‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                              <  1:126971/149‑1 (MQ=255)
                                                           aaCACAAACAACTTTTGTGTTGTGGGGGGCTTTTTTTTGGCGT‑TGCCCGCATGGCCCGAGAA‑AAAGTATTAAGGTTT‑CCCC‑TA‑AAAA‑A‑TT‑ACCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                                                                                <  3:531307/147‑1 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:886667/1‑151 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  1:425019/1‑151 (MQ=255)
                                                                aaacatgttttGTTGTGGGAGGGTTTGTTTTTGGGCGT‑TGCCCCCAGACCCCAAGGA‑AAATAAATAGGTTTC‑CCCC‑TA‑AAAA‑T‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                          <  2:358090/140‑1 (MQ=255)
                                                                aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                               >  3:722667/1‑143 (MQ=255)
                                                                  acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                     >  2:188868/1‑151 (MQ=255)
                                                                    aacGTTTGCTGTGTGGGGGTTTTTTTTTTGGCGT‑TCCCCCCATACCCCCGGTA‑AAAGTAATAGGTTTT‑CCCC‑TA‑ATAA‑A‑TT‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                   <  2:692794/151‑1 (MQ=255)
                                                                    aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                   >  4:259426/1‑151 (MQ=255)
                                                                     acttttGCGGTGGGGGGGGCTTTTTTTGGGCGA‑GCCCCGCACACCCCCAGGA‑AAAGTAAAAAGTTTT‑CCCT‑TA‑AAAA‑A‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                  <  1:275170/148‑1 (MQ=255)
                                                                          tgtttggtgAGGGGCTTTTTTTGGGCGG‑TTCCCGCATCCCCCCAGAT‑AAATAAATAGGTTTT‑CCCT‑TA‑AAAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                              <  2:956678/144‑1 (MQ=255)
                                                                           gCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                            >  3:589839/1‑151 (MQ=255)
                                                                            cTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                           >  2:592297/1‑151 (MQ=255)
                                                                             tgtgtgACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                          >  2:501037/1‑151 (MQ=255)
                                                                              gtgggAGGGGCTGTTTTTTGGCAA‑TGCCCGCATCGCCCAAGGA‑AAAGAAATAAGTTTT‑CCCT‑TA‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                            <  2:1287643/148‑1 (MQ=255)
                                                                                gtggggggttttttttgGGCGA‑TGCCCGCATCCCCCCAGAA‑AAAATAATAAGTTTT‑CCCC‑TA‑AAAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         <  3:1189433/137‑1 (MQ=255)
                                                                                  gaggggTTGTTTTTTGGCGA‑TGCCCGCATACCCCCAGAA‑AAAGTAA‑AAGGTTT‑TCCCTTA‑AAAA‑A‑TT‑CCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                        <  4:128086/145‑1 (MQ=255)
                                                                                  gaggggCTGTTTTTGGGGGT‑TGCCCGCATACCCCCAGGA‑AAAGTAATAAGTTTT‑CCCC‑TA‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                       <  1:72893/146‑1 (MQ=255)
                                                                                     ggttttttttttGGCTA‑TCCCCCCATCCCCCCAGAA‑AAAGAAATAGGTTTT‑CCCT‑TA‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   <  4:387407/143‑1 (MQ=255)
                                                                                     gggCTGTTTTTTGGCGA‑TGCCCGCATACCCCGAGAA‑AAAGTAATAGGTTTT‑CCCT‑TA‑ATAA‑T‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         <  3:754957/144‑1 (MQ=255)
                                                                                      ggCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   >  2:589399/1‑149 (MQ=255)
                                                                                        ctgctTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                                >  4:365404/1‑150 (MQ=255)
                                                                                            ttttgGGCGA‑TTTCCCAAAACCCCAAGAT‑AAAGTAATAGGTTTT‑CCCC‑TG‑ATAA‑A‑TC‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                              <  2:195021/148‑1 (MQ=255)
                                                                                             tttgGGCGA‑TGCCCGCATAGCCCCAGAT‑AAAGAAATAGGTTTT‑CCCC‑TA‑AAAA‑A‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                           <  1:248672/150‑1 (MQ=255)
                                                                                                gggCAA‑TGTCCCCATAGCCCCAGG‑‑AAAATTATAATTTTT‑TCCC‑TT‑AAAACC‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                       <  3:317248/150‑1 (MQ=255)
                                                                                                          ccGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                              >  2:645002/1‑151 (MQ=255)
                                                                                                          ccGCAGACCCCAAGAA‑AAAGTAATAGGTTTC‑CCCC‑TG‑TAAA‑A‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                 <  3:99341/148‑1 (MQ=255)
                                                                                                           cGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                             >  2:980101/1‑151 (MQ=255)
                                                                                                           cGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                             >  3:1289175/1‑151 (MQ=255)
                                                                                                            gCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                            >  3:1162237/1‑151 (MQ=255)
                                                                                                            gCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                            >  3:459416/1‑151 (MQ=255)
                                                                                                             cgaaaacccccagaaaaaTAATTAAGTTT‑ACCC‑TGAAAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                             <  4:354381/137‑1 (MQ=255)
                                                                                                               tAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                          >  2:1199666/1‑150 (MQ=255)
                                                                                                               ataccccAAGA‑AAAATAAAAAGGTTT‑TCCC‑TT‑AAAA‑ATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                            <  2:1293272/146‑1 (MQ=255)
                                                                                                                 gCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                       >  3:533409/1‑151 (MQ=255)
                                                                                                                   ccaaGGT‑AAAGAAATATGTTTT‑CCCA‑TA‑ATAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                         <  3:943367/144‑1 (MQ=255)
                                                                                                                    caaGGA‑AAATTAATAAGTTTT‑ACCT‑TA‑ATAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    <  3:317965/149‑1 (MQ=255)
                                                                                                                      aGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                   >  3:174536/1‑150 (MQ=255)
                                                                                                                       ggT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                   >  2:70101/1‑149 (MQ=255)
                                                                                                                        ggaaaaGAAATAAGGTTT‑CCCC‑TG‑AAAA‑A‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                                  <  2:926535/147‑1 (MQ=255)
                                                                                                                        gaaaaaGAAATAGGTTTT‑CCCC‑AA‑AAAA‑T‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                                 <  4:312856/148‑1 (MQ=255)
                                                                                                                         t‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  4:127795/1‑150 (MQ=255)
                                                                                                                         t‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  3:797926/1‑150 (MQ=255)
                                                                                                                           aaataaatagGTTTT‑CCCA‑TA‑ATAA‑T‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                <  4:299014/144‑1 (MQ=255)
                                                                                                                           aaaGAAAAAAGTTTT‑CCCT‑TA‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                               <  1:212697/150‑1 (MQ=255)
                                                                                                                            aaGTAATAAGTTTT‑CCCC‑TA‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             <  2:379780/151‑1 (MQ=255)
                                                                                                                            aaGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:87583/1‑151 (MQ=255)
                                                                                                                            aaGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  2:912591/1‑151 (MQ=255)
                                                                                                                             aGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  4:406715/1‑150 (MQ=255)
                                                                                                                                aataaGGTTT‑CCCC‑TA‑AAAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                         <  1:142424/151‑1 (MQ=255)
                                                                                                                                aaaaaGTTTT‑CCCC‑TA‑ATAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                            <  2:341359/145‑1 (MQ=255)
                                                                                                                                aaaaaGTTTT‑CCCC‑AG‑ATAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                            <  2:535647/145‑1 (MQ=255)
                                                                                                                                 aaaaGGTTT‑TCCC‑CT‑AAAA‑A‑TTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                         <  3:530645/149‑1 (MQ=255)
                                                                                                                                    aGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                       >  2:820410/1‑149 (MQ=255)
                                                                                                                                           cccc‑TA‑AAAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                 <  1:379228/148‑1 (MQ=255)
                                                                                                                                           aCCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                >  2:839337/1‑150 (MQ=255)
                                                                                                                                           aCCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg               >  2:1025478/1‑151 (MQ=255)
                                                                                                                                            ccc‑TA‑ATAA‑T‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                <  2:1233649/149‑1 (MQ=255)
                                                                                                                                            ccc‑TA‑ATAA‑A‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                 <  1:147637/148‑1 (MQ=255)
                                                                                                                                            ccc‑TA‑AAAA‑A‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg               <  2:301480/150‑1 (MQ=255)
                                                                                                                                              c‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa             >  2:604156/1‑150 (MQ=255)
                                                                                                                                                   aTAA‑A‑TT‑CCCGCACCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc            <  3:265724/148‑1 (MQ=255)
                                                                                                                                                    tAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt          >  3:1356961/1‑149 (MQ=255)
                                                                                                                                                     aa‑a‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg       <  1:136492/151‑1 (MQ=255)
                                                                                                                                                     aa‑a‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAGATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg       >  3:818282/1‑151 (MQ=255)
                                                                                                                                                        a‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg      >  2:185641/1‑150 (MQ=255)
                                                                                                                                                          tt‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTc    >  1:269778/1‑151 (MQ=255)
                                                                                                                                                           t‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTCATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt   >  2:658946/1‑151 (MQ=255)
                                                                                                                                                           t‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGGAATTAATAAGGCGCAGAGCGCCGTTTTGAACATAGGGTc    >  2:573567/1‑150 (MQ=255)
                                                                                                                                                             accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  2:1217135/1‑151 (MQ=255)
                                                                                                                                                             |                                                                                                                                                      
CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTT‑ACCC‑TG‑ATAA‑A‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378653‑4378952

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: ‑