Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,505 G→A 15.5% intergenic (+101/‑81) ydcD → / → yncI putative immunity protein for RhsE/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,5050GA15.5% 433.7 / 57.2 173intergenic (+101/‑81)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (72/73);  new base A (12/15);  total (85/88)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 180 total aligned reads displayed.

GAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGC  >  NC_000913/1530406‑1530604
                                                                                                   |                                                                                                   
 aacaacACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGa                                                                                                    >  1:2334835/1‑100 (MQ=255)
  acaacaacaGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAt                                                                                                   <  2:220461/100‑1 (MQ=38)
  acaacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAt                                                                                                   <  1:486976/100‑1 (MQ=38)
    aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAttt                                                                                                 <  2:931159/100‑1 (MQ=255)
     acacCAGAAAGAGAAAAATATAAACAACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                                <  2:399078/100‑1 (MQ=255)
     acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                                >  2:3477279/1‑100 (MQ=255)
      cacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTc                                                                                               >  2:3532369/1‑100 (MQ=255)
       acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCa                                                                                              <  2:600224/100‑1 (MQ=255)
       acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                              <  1:2066913/100‑1 (MQ=38)
        ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                             <  1:4101392/100‑1 (MQ=38)
        ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                             <  1:3860541/100‑1 (MQ=38)
        ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGAGCAAGATCATAAAAATTGTGATGTAAGTCACGATTTTAAt                                                                                             >  2:1749998/1‑100 (MQ=255)
        ccAGAAAGAGAAAAAGATAAAATACATACAGAAAGGCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATGTTCa                                                                                              >  2:839161/1‑99 (MQ=255)
         cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATc                                                                                            >  2:2260246/1‑100 (MQ=255)
         cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                            <  1:1582944/100‑2 (MQ=37)
          aGAAAGATAAAAATATAAAATACATACATAAAGACATGACATCAATGTACGGGTCAGTGCAATATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                             <  2:2069388/98‑1 (MQ=16)
          aGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                            <  2:1060540/99‑2 (MQ=37)
           gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCtt                                                                                          <  1:3292702/100‑1 (MQ=255)
           gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCtt                                                                                          >  1:4191020/1‑100 (MQ=255)
           gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtgat                                                                                          <  2:1658387/100‑4 (MQ=37)
            aaaGAGAAAAATATAAAATACATACAGAAAGCCATGACAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCttt                                                                                         <  2:2822030/100‑1 (MQ=255)
              agagAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGc                                                                                       <  1:577787/100‑1 (MQ=255)
               gagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCt                                                                                      <  1:3363654/100‑1 (MQ=255)
                agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                                     >  1:106422/1‑100 (MQ=255)
                agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAGCTTTGCt                                                                                      >  1:1311121/1‑99 (MQ=255)
                 gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCttt                                                                                    <  1:371271/100‑1 (MQ=255)
                  aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                                   <  2:407686/100‑1 (MQ=255)
                  aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                                   <  1:4008063/100‑1 (MQ=255)
                  aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtgatcttgccc                                                                                   <  1:817845/100‑11 (MQ=21)
                   aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                                  >  2:1744065/1‑100 (MQ=255)
                   aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                                  <  1:778417/100‑1 (MQ=255)
                   aaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTaa                                                                                  <  1:1111409/100‑1 (MQ=255)
                    aaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAc                                                                                 <  1:695718/100‑1 (MQ=255)
                       tataAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGcc                                                                              <  1:2906226/100‑1 (MQ=255)
                           aaaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACa                                                                          <  2:4006694/100‑1 (MQ=255)
                            cgatcttacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                                         <  1:1066330/94‑1 (MQ=255)
                             atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                        <  2:1799755/100‑1 (MQ=255)
                             atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                        <  1:1635929/100‑1 (MQ=255)
                               acatacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTg                                                                      <  2:1406236/100‑1 (MQ=255)
                                catacaGAACGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGCTGa                                                                     >  2:193530/1‑100 (MQ=255)
                                catacaGAAAGACATGAGAGCAAGGGACGGAGCGGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGACTACAGGTGa                                                                     >  2:2763783/1‑100 (MQ=255)
                                 atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGAt                                                                    <  2:3915235/100‑1 (MQ=255)
                                  tacaGCAAGACACGAGAGCAAGGCACGGTGCAGGGCTATATCATGAAAATTATGATGTAAATCACGATTGTCATCTTTGCTTTAACGCCTACAGGTTAt                                                                    <  2:3207016/99‑1 (MQ=255)
                                   acaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                                   <  2:3616071/99‑1 (MQ=255)
                                    caGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaa                                                                 >  1:165439/1‑100 (MQ=255)
                                     aGAAAGACATGCGTGCAAGGGACGGGTCAGGTCAAGATCATGAAAATTGTGATGTAAATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCa                                                                  <  2:3079991/98‑1 (MQ=255)
                                       aaaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAt                                                              <  2:1041975/100‑1 (MQ=255)
                                        aaGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAtt                                                             >  1:1446497/1‑100 (MQ=255)
                                             aTGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAcc                                                        <  1:2440473/100‑1 (MQ=255)
                                              tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                       >  1:3784333/1‑100 (MQ=255)
                                               gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                                      <  1:3805436/100‑1 (MQ=255)
                                               gagaGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCtt                                                      >  2:1239558/1‑100 (MQ=255)
                                                 gagCAAGGGACGGGTCAGGGCAAGATCAGGAAAATTGTGATGTAAAGCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                                    <  1:2681329/100‑1 (MQ=255)
                                                  agCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTgg                                                   <  1:1863545/100‑1 (MQ=255)
                                                   gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                                  >  1:2092006/1‑100 (MQ=255)
                                                    cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                                 >  1:1899115/1‑100 (MQ=255)
                                                     aaGGGACGGGGCAGGGCAAGATCATTAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                                                <  2:3549605/100‑1 (MQ=255)
                                                         gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAGTCACGATTTTCGTCTTTGCTTTAACGACTACAGGTGATCAAAATTTAACCTTTGCAATAAc                                             >  2:2234005/1‑99 (MQ=255)
                                                          aCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCAAGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                                             <  2:906106/98‑1 (MQ=255)
                                                           cGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaa                                          >  1:894540/1‑100 (MQ=255)
                                                            ggggCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaa                                         >  1:266441/1‑100 (MQ=255)
                                                             gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTTGAATAACTaaaa                                        <  1:2763790/100‑1 (MQ=255)
                                                              ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                       >  1:2296957/1‑100 (MQ=255)
                                                              ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                       <  2:2702907/100‑1 (MQ=255)
                                                              ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                       <  1:1506512/100‑1 (MQ=255)
                                                              ggcagagcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                       <  2:3172220/100‑1 (MQ=255)
                                                                    aagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTACATTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATcaa                                 >  2:1932978/9‑100 (MQ=18)
                                                                     caAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                                >  2:1332084/1‑100 (MQ=255)
                                                                     agcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                                >  1:3576321/8‑100 (MQ=25)
                                                                     agcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                                >  2:1549234/8‑100 (MQ=25)
                                                                       caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                              >  1:3513781/6‑100 (MQ=25)
                                                                       caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                              >  1:383376/6‑100 (MQ=25)
                                                                       caaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                              >  1:2171985/6‑100 (MQ=25)
                                                                       aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaaa                              >  1:2978603/1‑100 (MQ=255)
                                                                        gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                             <  2:337971/100‑1 (MQ=255)
                                                                        aaagaTGAAAATCGTGATTTACATCACAATTTTAATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                             >  2:2510045/5‑100 (MQ=25)
                                                                         aagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAgg                            >  1:2805829/4‑100 (MQ=37)
                                                                           cATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                          >  1:1950021/1‑100 (MQ=255)
                                                                            aTGCAAATTGTGATGTAAATCACGCTTTTCATCTTTGCTTTCACGCCTAAAGTAGATCAAGAATTACCCCACGGCATAACCAACAAGCTCCACAAGGAc                          >  2:3892745/1‑99 (MQ=255)
                                                                             tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc                        >  1:3505606/1‑100 (MQ=255)
                                                                                 agcaaagATGAAAATCGTGATTTACATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAgg                    >  1:347617/7‑100 (MQ=37)
                                                                                   ttGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTg                  <  1:1795616/100‑1 (MQ=255)
                                                                                     gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                <  2:1071489/100‑1 (MQ=255)
                                                                                     gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                <  2:4107338/100‑1 (MQ=255)
                                                                                     gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                <  2:1071462/100‑1 (MQ=255)
                                                                                     gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg                >  2:296635/1‑100 (MQ=255)
                                                                                       gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa              <  2:4189402/100‑1 (MQ=255)
                                                                                       gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa              >  1:1870635/1‑100 (MQ=255)
                                                                                        aTGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGATGCCAGGTGAGtct             <  1:2847666/100‑3 (MQ=255)
                                                                                        aTGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTAt             >  2:387195/1‑100 (MQ=255)
                                                                                          gTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTc           >  1:2999719/1‑100 (MQ=255)
                                                                                          gTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTc           <  1:1131738/100‑1 (MQ=255)
                                                                                            aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa         >  1:3662299/1‑100 (MQ=255)
                                                                                              aTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAg       <  2:591973/100‑1 (MQ=255)
                                                                                              aTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAg       <  2:2323591/100‑1 (MQ=255)
                                                                                               tCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGt      >  1:3903200/1‑100 (MQ=255)
                                                                                                cACGATTTTCATCTTTGCGTTAACGCCTACAGGTGATCAAAATTTCATCTTTGGAATACCTAAAAAGATAAAAAAGGACGCCAGGTGAGCATTCAAAGt      >  2:2970046/1‑99 (MQ=255)
                                                                                                 aCGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGttt    <  1:2277409/100‑1 (MQ=255)
                                                                                                 aCGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGttt    >  2:2411900/1‑100 (MQ=255)
                                                                                                   gATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGc  <  1:178229/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTGC  >  NC_000913/1530406‑1530604

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