Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,062,355 T→G 10.4% intergenic (+20/+36) asnV → / ← ldtA tRNA‑Asn/L,D‑transpeptidase ErfK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,062,3550TG10.4% 221.1 / 17.9 77intergenic (+20/+36)asnV/ldtAtRNA‑Asn/L,D‑transpeptidase ErfK
Reads supporting (aligned to +/- strand):  ref base T (37/32);  new base G (3/5);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

ACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAA  >  NC_000913/2062208‑2062504
                                                                                                                                                   |                                                                                                                                                     
aCCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCAGTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCt                                                                                                                                                     <  1:506009/150‑1 (MQ=255)
 cccTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTa                                                                                                                                                    >  1:1609413/1‑150 (MQ=255)
     ttGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGTCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCATAGTAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAg                                                                                                                                                <  2:1519695/150‑1 (MQ=255)
            cACCATTCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGTCttttt                                                                                                                                         <  1:1567919/150‑1 (MQ=17)
            cACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAACAAGCCTGCTAGAAAGCAGGCttttt                                                                                                                                         <  1:501199/150‑1 (MQ=17)
             aCCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTg                                                                                                                                        <  1:142837/150‑1 (MQ=255)
              ccATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCttttttg                                                                                                                                       <  1:598193/150‑3 (MQ=255)
              ccATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCttttttg                                                                                                                                       <  1:985810/150‑3 (MQ=255)
                  gCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTc                                                                                                                                   <  2:722432/150‑1 (MQ=255)
                  gCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTc                                                                                                                                   <  2:763666/150‑1 (MQ=255)
                    aCTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTa                                                                                                                                 <  1:14312/150‑1 (MQ=255)
                       gCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAAtt                                                                                                                              >  1:1676090/1‑150 (MQ=255)
                        ccATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTa                                                                                                                             <  1:1290441/150‑1 (MQ=255)
                          aTTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTAcc                                                                                                                           >  2:73229/1‑150 (MQ=255)
                               tatGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt                                                                                                                                        <  2:1964425/132‑2 (MQ=255)
                               tatGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTAGAAAGCAGGCtttttt                                                                                                                                        >  1:1964425/1‑131 (MQ=255)
                                  gcgcCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGctct                                                                                                                   <  1:1408140/150‑1 (MQ=255)
                                    gcCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCAGTATGTCACTGGGTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAACGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTa                                                                                                                 >  2:1376224/1‑150 (MQ=255)
                                     ccccGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTaa                                                                                                                >  1:1036279/1‑150 (MQ=255)
                                      cccGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTaaa                                                                                                               >  2:389155/1‑150 (MQ=255)
                                      cccGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTaaa                                                                                                               <  1:110703/150‑1 (MQ=255)
                                            cacaCGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCt                                                                                                         >  2:1019233/1‑150 (MQ=255)
                                                  aTTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTg                                                                                                   <  2:528748/150‑1 (MQ=255)
                                                    tCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGaa                                                                                                 <  2:1676090/150‑1 (MQ=255)
                                                     ccTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAAc                                                                                                <  2:792978/150‑1 (MQ=255)
                                                      ctctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAAcc                                                                                               >  2:176183/1‑150 (MQ=255)
                                                      ctctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAAcc                                                                                               <  1:457430/150‑1 (MQ=255)
                                                        ctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAg                                                                                             >  2:1643135/1‑150 (MQ=255)
                                                          gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCttt                                                                                                                                    <  2:741479/109‑1 (MQ=255)
                                                          gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCttt                                                                                                                                    >  1:741479/1‑109 (MQ=255)
                                                          gTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGaa                                                                                           >  2:317396/1‑150 (MQ=255)
                                                            aGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACt                                                                                         <  1:1778351/150‑1 (MQ=255)
                                                                cAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAAtt                                                                                     >  2:282693/1‑150 (MQ=255)
                                                                  gTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTg                                                                                   >  1:566818/1‑150 (MQ=255)
                                                                  gTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTATAATTACAAACTCTCTTAAAACATCTGTCTTGAACCAGACCTAATTTg                                                                                   >  2:1849526/1‑150 (MQ=255)
                                                                   tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGc                                                                                  <  1:389155/150‑1 (MQ=255)
                                                                               cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTccc                                                                      <  1:1182658/150‑1 (MQ=255)
                                                                               cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTccc                                                                      <  1:73229/150‑1 (MQ=255)
                                                                                    tGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCg                                                                 <  1:1423076/150‑1 (MQ=255)
                                                                                        aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGt                                                             >  2:475019/1‑150 (MQ=255)
                                                                                           ccGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAACCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGc                                                          <  1:622679/150‑1 (MQ=255)
                                                                                              tATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAAc                                                       <  2:566818/150‑1 (MQ=255)
                                                                                                tGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgc                                                     >  2:657314/1‑150 (MQ=255)
                                                                                                 gTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgca                                                    >  2:289776/1‑150 (MQ=255)
                                                                                                 gTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACgca                                                    <  1:1643135/150‑1 (MQ=255)
                                                                                                     cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAt                                                >  1:1548102/1‑150 (MQ=255)
                                                                                                       ggTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGaa                                                                                                 >  2:1063117/1‑99 (MQ=255)
                                                                                                       ggTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGaa                                                                                                 <  1:1063117/99‑1 (MQ=255)
                                                                                                       ggTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAttt                                              >  2:1432237/1‑150 (MQ=255)
                                                                                                       ggTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAttt                                              >  1:1884063/1‑150 (MQ=255)
                                                                                                       ggTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAttt                                              >  1:509648/1‑150 (MQ=255)
                                                                                                       ggTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCAttt                                              >  2:1999171/1‑150 (MQ=255)
                                                                                                         ttCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTgc                                            <  1:1619622/150‑1 (MQ=255)
                                                                                                          tctaGTCCAGTCAGAGGAGCCAAATTCAAAAAAGGCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATATGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTgcg                                           <  2:1022713/147‑1 (MQ=255)
                                                                                                          tCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTgcg                                           >  2:904570/1‑150 (MQ=255)
                                                                                                            gAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGcgcg                                         >  2:1880341/1‑150 (MQ=255)
                                                                                                            gAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGcgcg                                         >  2:756999/1‑150 (MQ=255)
                                                                                                             aGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGa                                        >  1:1817126/1‑150 (MQ=255)
                                                                                                              gTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGAt                                       >  1:1684983/1‑150 (MQ=255)
                                                                                                                  agacAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTAc                                   >  1:1419698/4‑150 (MQ=255)
                                                                                                                  aGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTAc                                   <  2:607615/150‑1 (MQ=255)
                                                                                                                   gTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGACCCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACa                                  >  1:1360337/1‑150 (MQ=255)
                                                                                                                      agagGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCa                               <  1:1782742/150‑1 (MQ=255)
                                                                                                                      agagGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCa                               >  2:818278/1‑150 (MQ=255)
                                                                                                                      agagGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCa                               <  1:1789821/150‑1 (MQ=255)
                                                                                                                         ggagCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAAct                            <  1:1678368/150‑1 (MQ=255)
                                                                                                                              cAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTg                       >  2:1106751/1‑150 (MQ=255)
                                                                                                                              cAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACACGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTg                       >  2:324188/1‑150 (MQ=255)
                                                                                                                               aaaTTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTATCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAAATTGCACAGGCATTCCCGATCGACGTAGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGc                      >  2:1656413/1‑150 (MQ=255)
                                                                                                                                   tCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGt                  >  1:630063/1‑150 (MQ=255)
                                                                                                                                    cAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGtt                 >  2:604067/1‑150 (MQ=255)
                                                                                                                                     caaaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTg                >  1:307272/2‑150 (MQ=255)
                                                                                                                                     aaaaaaGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTg                <  1:1019233/150‑1 (MQ=255)
                                                                                                                                       aaaaGCCTGCTAGAAAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGAt              >  2:1349896/1‑150 (MQ=255)
                                                                                                                                             cTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGc        <  1:904570/150‑1 (MQ=255)
                                                                                                                                              tGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGcc       >  1:1606371/1‑150 (MQ=255)
                                                                                                                                                cTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCg     <  1:604067/150‑1 (MQ=255)
                                                                                                                                                   tCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCaa  >  2:1960402/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCCTTTGACATCACCATGCACTGCCATTAATATGCGCCCCGTTCACACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCAAAAAAGCCTGCTTTCTAGCAGGCTTTTTGCTTTCTAATTACCAACGCTCTTAAAACATCTGTCTTGAACCAGAACTAATTTGCACAGGCATTCCCGATCGACGTTGCAACGCAGCATTTGCGCGATTTACATCAACTTCTTGCCCGTTGATAAACGCCCGCAA  >  NC_000913/2062208‑2062504

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