Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 729,379 A→C 11.2% intergenic (+39/‑204) ybfA → / → rhsC DUF2517 domain‑containing protein YbfA/rhs element protein RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913729,3790AC11.2% 325.5 / 35.4 125intergenic (+39/‑204)ybfA/rhsCDUF2517 domain‑containing protein YbfA/rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (55/56);  new base C (4/10);  total (59/66)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 129 total aligned reads displayed.

GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATC  >  NC_000913/729231‑729528
                                                                                                                                                    |                                                                                                                                                     
ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAAc                                                                                                                                                      >  2:841468/1‑150 (MQ=255)
     gACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                                                                                                                                                 >  2:150095/1‑150 (MQ=255)
      aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGc                                                                                                                                                <  1:1276701/150‑1 (MQ=255)
      aCCGCGCCCGCGTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGc                                                                                                                                                <  2:648176/150‑1 (MQ=255)
       ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCg                                                                                                                                               <  2:1169554/150‑1 (MQ=255)
        cgcgCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCgg                                                                                                                                              <  2:890973/150‑1 (MQ=255)
         gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGt                                                                                                                                             <  1:813251/150‑1 (MQ=255)
         gcgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGGGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGt                                                                                                                                             <  1:1269349/150‑1 (MQ=255)
          cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                                            <  2:446409/150‑1 (MQ=255)
           gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCTAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGttt                                                                                                                                           >  1:439402/1‑150 (MQ=255)
           gcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGttt                                                                                                                                           <  2:1041240/150‑1 (MQ=255)
            cccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCGGtttt                                                                                                                                          <  2:17462/150‑1 (MQ=255)
             ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGttttt                                                                                                                                         <  2:1224854/150‑1 (MQ=255)
             ccGCTGTTACAGCTACCTGCATCGCGTCTGGGCGAAACCCAGCGATAAACCGGTGGGGATGGATCAGGCCGAAAAGGCAACTGGTGCTTTTTATTGATTGCTGAGTTTTGAATACATTGACAACCGCCAGCGTTTCCGCTGGCGGttttc                                                                                                                                         <  2:1537690/150‑2 (MQ=255)
              cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCTAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                                                        <  2:439402/150‑1 (MQ=255)
               gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTg                                                                                                                                       <  2:81955/150‑1 (MQ=255)
                cTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCAGGCGGTTTTTTGt                                                                                                                                      >  2:1539322/1‑150 (MQ=255)
                 ttttACAGCTACCTGCATCGCGCCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTa                                                                                                                                     <  2:425579/150‑1 (MQ=255)
                   ttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTAta                                                                                                                                   >  2:1500292/1‑149 (MQ=255)
                    tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                                                  >  1:1346364/1‑150 (MQ=255)
                    tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACAGGTGAATTTTATTGATTGCTGAGTTTTGAATACAATGAAAACCGCCAGAGGAAACGCTGGCGCTTTTTTGTATCa                                                                                                                                  >  1:728292/1‑150 (MQ=255)
                     aCAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAg                                                                                                                                 >  1:966518/1‑150 (MQ=255)
                      cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAgg                                                                                                                                <  2:789024/150‑1 (MQ=255)
                       aGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGa                                                                                                                               <  1:431640/150‑1 (MQ=255)
                            ccTGCATCGCGTCTTGTCGAAAACCATCGATAAACCGGTGTTGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg                                                                                                                          <  1:314193/150‑1 (MQ=255)
                                 aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCag                                                                                                                     >  2:422378/1‑150 (MQ=255)
                                 aTCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCAGTTTTTTGTATCAGGAAGATGTTCag                                                                                                                     <  2:1098670/150‑1 (MQ=255)
                                   cgcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                                                   <  1:674885/150‑1 (MQ=255)
                                        cTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAata                                                                                                              >  1:1216326/1‑150 (MQ=255)
                                         tGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAatat                                                                                                             <  2:1536183/150‑1 (MQ=255)
                                         tGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAatat                                                                                                             >  1:89326/1‑150 (MQ=255)
                                          ggTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAtt                                                                                                            >  1:392930/1‑150 (MQ=255)
                                           gtggAAAACCAGCGATAAACCGGTCTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                           <  2:1422752/147‑1 (MQ=255)
                                           gTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATAttt                                                                                                           >  1:818591/1‑150 (MQ=255)
                                             cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                         <  2:1594713/150‑1 (MQ=255)
                                             agAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTCGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                                         >  2:1074503/2‑150 (MQ=255)
                                              gAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGa                                                                                                        >  1:1206342/1‑150 (MQ=255)
                                                  aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                                                    >  2:1230717/1‑150 (MQ=255)
                                                   ccAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                                                   >  1:486034/1‑150 (MQ=255)
                                                   ccAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                                                   >  2:1411379/1‑150 (MQ=255)
                                                      gCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTta                                                                                                >  2:1104663/1‑150 (MQ=255)
                                                        gATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtata                                                                                              <  2:704180/150‑1 (MQ=255)
                                                        gATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTtata                                                                                              >  1:1396008/1‑150 (MQ=255)
                                                         aTAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAc                                                                                             >  2:81950/1‑150 (MQ=255)
                                                            aaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTCGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCaa                                                                                          <  1:1074503/150‑1 (MQ=255)
                                                             aCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                                                         <  1:144198/150‑1 (MQ=255)
                                                                  tgtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                    >  1:525887/1‑150 (MQ=255)
                                                                  tgtgGATGGATCAGGCCGAAAACGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGaaa                                                                                    <  2:1333112/150‑1 (MQ=255)
                                                                   gtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                                                   >  1:275102/1‑150 (MQ=255)
                                                                   gtgGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAAc                                                                                   <  1:514437/150‑1 (MQ=255)
                                                                    tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACa                                                                                  >  2:290628/1‑150 (MQ=255)
                                                                    tggatggaTCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGGTCAGAGGAATATTTAGACATCGTTATAGCAATCGAAACa                                                                                  >  2:1182150/1‑150 (MQ=255)
                                                                     ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                 >  1:371788/1‑150 (MQ=255)
                                                                     ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                 <  1:290628/150‑1 (MQ=255)
                                                                     ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                                                 <  1:1187825/150‑1 (MQ=255)
                                                                        tggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                                              >  2:1084332/1‑150 (MQ=255)
                                                                          gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                                            >  2:1318259/1‑150 (MQ=255)
                                                                             cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCg                                                                                                                                               <  2:481767/80‑1 (MQ=38)
                                                                             cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGTTTCCGCTGGCg                                                                                                                                               >  1:481767/1‑80 (MQ=38)
                                                                               ggCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGa                                                                       >  2:643690/1‑150 (MQ=255)
                                                                               ggCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGa                                                                       <  1:842719/150‑1 (MQ=255)
                                                                                     aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                                                 <  1:150095/150‑1 (MQ=255)
                                                                                     aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                                                 <  1:844929/150‑1 (MQ=255)
                                                                                     aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                                                 <  1:904671/150‑1 (MQ=255)
                                                                                                   aTTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTcc                                                   <  2:236846/150‑1 (MQ=255)
                                                                                                    tttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATTTTTTcgt                                                  >  2:1126432/1‑148 (MQ=255)
                                                                                                    tttttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                                                  >  2:189171/1‑150 (MQ=255)
                                                                                                      tttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTaa                                                <  2:371788/150‑1 (MQ=255)
                                                                                                       ttATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAt                                               >  2:14031/1‑150 (MQ=255)
                                                                                                         attgattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                                             <  1:1084332/150‑1 (MQ=255)
                                                                                                           tgattgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGtt                                           >  1:44361/1‑150 (MQ=255)
                                                                                                              ttgCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCtt                                        <  2:1396008/150‑1 (MQ=255)
                                                                                                                 cTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATc                                     <  2:533469/150‑1 (MQ=255)
                                                                                                                     tttttGAATACATTGAAAACCGTCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGATGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATTTATATAAATATTTTCCTAATTGTTCTTATCTGAc                                 <  2:1472644/149‑1 (MQ=255)
                                                                                                                     gTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGAc                                 >  1:967701/1‑150 (MQ=255)
                                                                                                                       tttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAg                               <  2:1206342/150‑1 (MQ=255)
                                                                                                                        ttGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGa                              <  2:365715/150‑1 (MQ=255)
                                                                                                                              aCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGACGAATATTTAGACATCGTTATACCAATCGAAACATATATTCCTGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAtctc                        <  2:910956/150‑1 (MQ=255)
                                                                                                                                aTTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCAc                      >  1:1588746/1‑150 (MQ=255)
                                                                                                                                aTTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCAc                      >  2:483801/1‑150 (MQ=255)
                                                                                                                                   gAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTa                   >  2:933697/1‑150 (MQ=255)
                                                                                                                                   aaaaaCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTa                   >  1:122864/2‑150 (MQ=255)
                                                                                                                                     aaaCCGCCAGCTGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                 <  1:383775/150‑1 (MQ=255)
                                                                                                                                     aaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                 <  1:422378/150‑1 (MQ=255)
                                                                                                                                     aaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                 >  2:983417/1‑150 (MQ=255)
                                                                                                                                     aaaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAg                 <  2:1216326/150‑1 (MQ=255)
                                                                                                                                      aaCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAgg                >  2:870191/1‑150 (MQ=255)
                                                                                                                                       aCCGCCAGCGTTTCCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGc               >  1:1298329/1‑150 (MQ=255)
                                                                                                                                       aCCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGc               <  1:643690/150‑1 (MQ=255)
                                                                                                                                        ccGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTaa                  >  1:1359019/1‑146 (MQ=255)
                                                                                                                                        ccGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTaa                  <  2:1359019/146‑1 (MQ=255)
                                                                                                                                            cAGCGGAAACGCTGGCTGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt          <  1:1053648/150‑1 (MQ=255)
                                                                                                                                            cAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCt          <  2:44361/150‑1 (MQ=255)
                                                                                                                                             aGCTGAAACTCAGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCtt         <  1:1539322/150‑1 (MQ=255)
                                                                                                                                              gCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTta        <  2:89326/150‑1 (MQ=255)
                                                                                                                                               cGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTtat       >  2:242739/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATaaa    <  2:525887/150‑1 (MQ=255)
                                                                                                                                                  aaaCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATaaa    <  1:1121330/150‑1 (MQ=255)
                                                                                                                                                    aCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATc  <  1:81950/150‑1 (MQ=255)
                                                                                                                                                    aCGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATc  <  2:1588746/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATATCTCACTTAAGGCTTTCTTATAAATC  >  NC_000913/729231‑729528

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