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

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378653‑4378952
                                                                                                                                                     |                                                                                                                                                       
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTac                                                                                                                                                        >  2:728897/1‑151 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTa                                                                                                                                                         <  2:447609/150‑1 (MQ=255)
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGCTTACCCTGATAAATTac                                                                                                                                                        >  2:737338/1‑151 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTacc‑g                                                                                                                                                     >  1:147654/1‑151 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTacc‑gcc                                                                                                                                                   <  1:1426271/150‑1 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTacc‑gcc                                                                                                                                                   <  2:947458/150‑1 (MQ=255)
           aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTacc‑gccaccgcc                                                                                                                                             >  2:897820/1‑150 (MQ=255)
           aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTacc‑gccaccgcc                                                                                                                                             >  1:421678/1‑150 (MQ=255)
             aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTCCC‑GCCACCGCCAt                                                                                                                                           <  2:103764/150‑1 (MQ=255)
              gCGTCGCCCCCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATAc                                                                                                                                         <  1:447353/151‑1 (MQ=255)
                    ncaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGNNCGNNTAGCCCCAGGTAAAGTAATAAGGTTNNCCCTGATAAATTACC‑GCCACCGCCATACCCggg                                                                                                                                    >  1:105850/2‑150 (MQ=255)
                    ccaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCggg                                                                                                                                    >  2:961101/1‑150 (MQ=255)
                      cccaacccccccAATAATTTTTTTGCTGACCGAACCTAACACAAACACCGTTTGCTGTGGGACGGGCTGCTTTTTGGCGATGTCCGCATACCCCCAGAAAAAGAAAAAAGGTTTCCCTTAATAATTTCCC‑GCCACCGCCATACCCGGGATc                                                                                                                                 <  1:1192169/141‑1 (MQ=255)
                      cccaAGCCCGCCAGTAAATGTTTTGCTGACCGACCCAAACTAAAACACCTTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGCCCCCATACCCCCAGGTAAAGAAATAAGGTTTCCCTTAATAAATTCCC‑GCCACCGCCATACCCGGGATc                                                                                                                                 <  1:298538/150‑1 (MQ=255)
                       ccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCCTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCt                                                                                                                                >  2:236826/1‑151 (MQ=255)
                           gCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTttg                                                                                                                            >  1:430470/1‑151 (MQ=255)
                             ccGCCAGTAAATTTTTTGCTGCCCGAACCTAACTCAAACACCGTTTGCTGTGTGACGGGTTGTTTTTGGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGGTTTCCCCTAATAAATTCCC‑GCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          <  2:1349667/150‑1 (MQ=255)
                               gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                        >  2:1028686/1‑151 (MQ=255)
                               gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                        >  1:342510/1‑151 (MQ=255)
                                  aGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                     >  1:913905/1‑151 (MQ=255)
                                    tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:594763/1‑151 (MQ=255)
                                     aatttttttGCGGACCGACCCAAACTCAAACACCTTTTGTTGTGTGAGGGGTTGTTTTTGGGCGATGCCCGCATAGCCCAAGATAAAGTAATAAGTTTTCCCTTGAAAAATCCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   <  1:595198/145‑1 (MQ=255)
                                     aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:603824/1‑150 (MQ=255)
                                     aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:1007000/1‑150 (MQ=255)
                                     aaaTGTTTTGCTGACCAAACCTAACTAAAACAACGTTTGCTGTGTGGCGGGCTTTTTTTGGGCGATCCCCCCATAGCCCCAGGTAAAGTAATAAGTTTTCCCCTAAAAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  <  1:331413/151‑1 (MQ=255)
                                     aaaTGTTTGGCTGACCAAACCAAACTAAAAAAGGTTTTGCTGTGGGGCGGGTTGTTTTTGGGGGATGTCCGAATACCCCGGGAAAAAGAAAAAGGGTTTCCCTTGATAAATGACC‑GCCGCCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                  <  2:939298/151‑1 (MQ=255)
                                         tttttGCTGACCGACCCTAACTAAAACAACTTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGCCCGCATAGCCCCAGGAAAAGTAATAAGTTTTCCCCTAATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                               <  2:1140360/149‑1 (MQ=255)
                                              gCTGCCCGACCCTACCCCAAACAGCTTTTGCTGTGTGAGGGGTTGTTTTTGGGCAATGCCCGCATAGCCCCAGTAAAATTAATAA‑TTTTTCCCGTAAAAATTACCGGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         <  1:1416459/151‑1 (MQ=255)
                                              gCTGACCGACCCTAACTCAAACAACTTTTGCTGTGTGAGGGGCTGTTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGTTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                          <  2:587227/150‑1 (MQ=255)
                                                 gACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                      >  2:868091/1‑151 (MQ=255)
                                                     gAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                                  >  1:1319340/1‑151 (MQ=255)
                                                       aCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                                 >  2:658873/1‑150 (MQ=255)
                                                           aacaaaaaCACCTTTTGCTGGGTGACGGGCTGTTTTTTGGCGATGTCCGCATACCCCCAGTAAAATTAAAAAGTTTTCCCTTAATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                            <  2:490911/146‑1 (MQ=255)
                                                           aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             >  1:200723/1‑150 (MQ=255)
                                                            acaaaaaCAACTTTTGCTGTGTGAGGGGCTTTTTTTTGGCGATGTCCGCATACCCCAAGGTAAAGAAAAAGGGTTTCCCCTGAAAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                           <  2:41971/147‑1 (MQ=255)
                                                            aCTAAAACAACTTTTGTTGTGTACCGGGTTTTTTTTTGGCGTTGCCCCCATACCCCCAGGAAAAGTAATAAGTTTTCCCCTGAAAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                           <  1:1353355/151‑1 (MQ=255)
                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                          >  1:700516/1‑150 (MQ=255)
                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                          >  1:882515/1‑150 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:843279/1‑151 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:237938/1‑151 (MQ=255)
                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  2:331346/1‑151 (MQ=255)
                                                               aaaaCACCGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGGTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        <  2:1167699/150‑1 (MQ=255)
                                                                aaccaccttttGCTGTGTGAGGGGCTTTTTTTTGGCTATGTCCCCATCCCCCAAGGAAAATAAATAAGGTTTCCCCTGAAAAATCCCC‑CCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                       <  2:943690/143‑1 (MQ=255)
                                                                 aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                      >  1:20649/1‑151 (MQ=255)
                                                                  acagcGTTTGCTGGGTGACGGGCTGCTTTGGGGCGATGGCGGAATCGCCCCGGGAAAAATAATAAGTTTTCCCCTAATAAATGCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                     <  2:1138151/151‑1 (MQ=255)
                                                                       gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                 >  1:727929/1‑150 (MQ=255)
                                                                          tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                             >  2:492126/1‑151 (MQ=255)
                                                                          tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                             >  2:181206/1‑151 (MQ=255)
                                                                           gCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                            >  2:963712/1‑151 (MQ=255)
                                                                                gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         >  1:351145/1‑149 (MQ=255)
                                                                                    ggggtttttttttGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGTTTTCCCCTAAAAAATTCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   <  2:1129024/143‑1 (MQ=255)
                                                                                    ggggttttttttgGGCGATGTCCGCATAGCCCAAGGAAAAGTAAAAGGTTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   <  2:1006024/143‑1 (MQ=255)
                                                                                     gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   >  2:1006913/1‑150 (MQ=255)
                                                                                       gcttttttttGGCGTGGTCCGAATAGCCCCGGGAAAAGAAATAAGGTTTCCCCTAATAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                                <  2:1400060/146‑1 (MQ=255)
                                                                                         ttttttttGGCGATGCCCGCATACCCCCAGAAAAATTAATAAGTTTTCCCCTGATAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               <  1:773313/147‑1 (MQ=255)
                                                                                         tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               >  2:1345765/1‑150 (MQ=255)
                                                                                          gttttttGGCGATGCCCGCATACCCCGAGGTAAAGTAATATGTTTTCCCCTAATAATTTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               <  1:222620/147‑1 (MQ=255)
                                                                                             tttgGGCGTTGCCCGCATAGCCCCGGGAAAAGTAATAAGTTTTCCCCTAATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                           <  2:52838/150‑1 (MQ=255)
                                                                                               ttGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGGTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                         <  2:1376375/150‑1 (MQ=255)
                                                                                               tgggCGATGTCCGCATAGCCCCAGGTAAAGTAAAAAGGTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                         <  2:811784/148‑1 (MQ=255)
                                                                                                 ggCGATGTCCGCATAGCCCAAGGTAAAGTAATAAGGTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAg                                                      <  2:493826/151‑1 (MQ=255)
                                                                                                  gCGATGCCCGCATAGCCCCAGGTAAAGTAATAAGGTTTCCCTTAAAAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                     <  2:546994/151‑1 (MQ=255)
                                                                                                  gCGATGCCCGCATACCCCCAGGTAAAGAAATAAGTTTTCCCCTAATAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                     <  2:1215714/151‑1 (MQ=255)
                                                                                                    gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                    >  2:41310/1‑150 (MQ=255)
                                                                                                    aaggTCCGCATACCCCGGGAAAAAATAATAGGTTTTCCCATTAAAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                   <  2:788448/148‑1 (MQ=255)
                                                                                                     aTGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                    >  1:344417/1‑149 (MQ=255)
                                                                                                        tCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCACCAAATGTGGAGCAAGAgg                                               >  1:911255/1‑151 (MQ=255)
                                                                                                        cccacATACCCCCAGGAAAAGTAATAGGTTTTCCCTTAATAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                               <  1:1055220/147‑1 (MQ=255)
                                                                                                         ccGCATAGCCCCAGGAAAAGTAATAGGTTTTCCCCTGATAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                              <  2:534510/151‑1 (MQ=255)
                                                                                                           gCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                            >  2:1451575/1‑151 (MQ=255)
                                                                                                             aTAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                           >  1:1397297/1‑150 (MQ=255)
                                                                                                               aGCCCGAGGAAAATAAAAAAGTTTTCCCCTAAAAAATCCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                        <  1:401913/151‑1 (MQ=255)
                                                                                                               aGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                        >  2:400238/1‑151 (MQ=255)
                                                                                                                gCCCCAGGAAAATAAATAGGTTTTCCCCTGATAATTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                       <  2:1309857/151‑1 (MQ=255)
                                                                                                                 ccccAGGAAAAGAAAAAAGTTTTCCCCTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                      <  2:1039390/151‑1 (MQ=255)
                                                                                                                  cccAGGAAAAGAAATAGGGTTTCCCTTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                     <  2:175217/151‑1 (MQ=255)
                                                                                                                   caggttAAAGTAATAGGTTTTCCCTTGAAAATTCCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    <  2:407546/146‑1 (MQ=255)
                                                                                                                   caagatAAAGAAATAAGTTTTCCCCTAATAAATAACC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    <  1:742548/146‑1 (MQ=255)
                                                                                                                      ggaaaaaTAATAAGTTTTACCCTGATAAATTACC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                                 <  1:803426/148‑1 (MQ=255)
                                                                                                                      ggTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                   >  1:1319197/1‑149 (MQ=255)
                                                                                                                      gaaaaaGTAAAAAGGTTTCCCTGGATATTTTCCC‑GCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                                 <  2:319886/148‑1 (MQ=255)
                                                                                                                       gTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  2:553353/1‑151 (MQ=255)
                                                                                                                       gTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  2:321778/1‑151 (MQ=255)
                                                                                                                            gTAATAGGTTTTCCCCTAATAAATTCCC‑GCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                           <  1:475903/151‑1 (MQ=255)
                                                                                                                            gTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  2:361272/1‑149 (MQ=255)
                                                                                                                            gTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:190828/1‑149 (MQ=255)
                                                                                                                             taataaGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:1289951/1‑148 (MQ=255)
                                                                                                                              aataaGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                           >  2:353206/1‑149 (MQ=255)
                                                                                                                               ataaGGTTTCCCCTGATAAATTCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                        <  2:438556/151‑1 (MQ=255)
                                                                                                                                  aGGTTTANCTTGATAAATTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGNNANNCACAATATCGTTGATNNNNNNNGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                     <  2:105846/151‑1 (MQ=255)
                                                                                                                                       tccccTAAAAAATTCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                 <  1:1418795/148‑1 (MQ=255)
                                                                                                                                        aCCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg               >  1:780918/1‑151 (MQ=255)
                                                                                                                                             gATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt          >  2:119403/1‑151 (MQ=255)
                                                                                                                                             gATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt          >  2:334236/1‑151 (MQ=255)
                                                                                                                                              aTAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGTACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATCAGGCGCAGAGCGTCGTTTTGAACAt          >  1:876075/1‑150 (MQ=255)
                                                                                                                                                aaaTTCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg       <  1:737345/151‑1 (MQ=255)
                                                                                                                                                aaaTTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg       >  2:339245/1‑151 (MQ=255)
                                                                                                                                                 aaTTCCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg      <  2:353707/151‑1 (MQ=255)
                                                                                                                                                  aTTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt     >  1:70872/1‑151 (MQ=255)
                                                                                                                                                     acc‑gccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  2:172281/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378653‑4378952

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