Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 379,246 (G)10→9 intergenic (‑120/+47) yaiN ← / ← yaiO putative alpha helix chain/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913379,2370G.84.3% 208.6 / 38.5 71intergenic (‑111/+56)yaiN/yaiOputative alpha helix chain/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (7/4);  new base . (43/16);  total (50/21)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

ATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCA  >  NC_000913/379147‑379327
                                                                                          |                                                                                           
aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                           <  1:2280752/93‑1 (MQ=255)
aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                           <  1:2280753/93‑1 (MQ=255)
 ttCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                          <  2:3100820/93‑1 (MQ=255)
 ttCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                          >  2:182664/1‑93 (MQ=255)
  tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                         >  2:1142574/1‑93 (MQ=255)
  tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                         <  1:4587040/93‑1 (MQ=255)
  tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                         >  1:3153410/1‑93 (MQ=255)
  tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                         <  2:4703083/93‑1 (MQ=255)
  tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                         >  2:2468558/1‑93 (MQ=255)
   cTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggg                                                                                        <  1:4343404/93‑1 (MQ=255)
    tataTGTCAATGCATACCCCCTTATAGTATATAACTCTGGAGAATAAAAGTTATCGGGTGAAGGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                       <  1:733439/93‑1 (MQ=255)
    tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                       >  2:2132405/1‑93 (MQ=255)
          tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggctaa                                                                                <  2:4670920/93‑5 (MQ=255)
          tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggctaa                                                                                <  2:1491234/93‑5 (MQ=255)
             aTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATc                                                                             >  1:3192347/1‑93 (MQ=255)
             aTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCg                                                                            <  1:3540161/93‑1 (MQ=255)
               gCATACCCCCCTATATTAAAAAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGc                                                                          <  1:3522301/93‑1 (MQ=255)
               gCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGc                                                                          <  1:1383618/93‑1 (MQ=255)
                 aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAg                                                                        >  1:3578892/1‑93 (MQ=255)
                 aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAg                                                                        >  1:4269238/1‑93 (MQ=255)
                 aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAgg                                                                       >  1:1414167/1‑93 (MQ=255)
                  tACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAgg                                                                       >  1:3479349/1‑93 (MQ=255)
                   aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAgg                                                                       >  2:710752/1‑93 (MQ=255)
                   aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAggg                                                                      >  1:1995313/1‑93 (MQ=255)
                      ccccTGTAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGTTTAATGGGAGGGGag                                                                    >  2:4710704/1‑92 (MQ=255)
                      ccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAg                                                                   >  2:977195/1‑93 (MQ=255)
                      ccccTATAGTATATAAAACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAg                                                                   <  1:3737032/93‑1 (MQ=255)
                        ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgc                                                                 >  2:2970718/1‑93 (MQ=255)
                         cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGgc                                                                 >  1:4024051/1‑93 (MQ=255)
                         cTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAAAAAAGAGAAAAGGGGGGGGGGCTAATCGGCAGGGAAGgc                                                                 <  2:4527514/93‑1 (MQ=255)
                          tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgggcc                                                               >  2:1816156/1‑90 (MQ=255)
                          tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCACGGAAggtcc                                                               >  1:3846777/1‑90 (MQ=255)
                           ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccgc                                                              >  2:3883463/1‑93 (MQ=255)
                            taGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGtgcgcc                                                             >  2:2799629/1‑93 (MQ=255)
                            taGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGccc                                                            >  2:90352/1‑93 (MQ=255)
                             aGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGccc                                                            >  1:2504837/1‑93 (MQ=255)
                              atatataACACGGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGcccc                                                           <  2:2996924/92‑1 (MQ=255)
                               tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGGccccccc                                                          >  2:3333679/1‑92 (MQ=255)
                                 tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCggc                                                        >  2:1017080/1‑92 (MQ=255)
                                 tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCCCCCCGGa                                                        >  2:1017081/1‑93 (MQ=255)
                                  aaaaCACTGGAAAAAAAAATTTAACCGGGGAATGTTGGCCGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGa                                                        <  1:2462179/91‑1 (MQ=255)
                                   taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATa                                                      >  2:1854028/1‑93 (MQ=255)
                                   taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATa                                                      >  2:1615362/1‑93 (MQ=255)
                                   taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATa                                                      >  1:2301673/1‑93 (MQ=255)
                                     acacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc                                                    <  1:2671139/93‑1 (MQ=255)
                                      cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                   >  1:1632425/1‑93 (MQ=255)
                                      cacTGAAGAATAAAATTAATCCGGTGAATGTGGCCGGAAAAAAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                   <  1:3536933/93‑1 (MQ=255)
                                       acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGGCCCGGAGAgcgg                                                  >  1:635698/1‑93 (MQ=255)
                                       acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                  <  1:650327/93‑1 (MQ=255)
                                        cTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCGAATCGGCAGGGAAAGAAACCCCCGATAgcggg                                                 >  2:380677/1‑93 (MQ=255)
                                         tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACCAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGGCGCCCCGGGCAgcggg                                                 >  2:229341/1‑93 (MQ=255)
                                             gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCTCCCCTGATAGCGGGCGGCa                                            >  2:2469601/1‑93 (MQ=255)
                                             gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCa                                            >  2:1649249/1‑93 (MQ=255)
                                             gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCCGGAAACGGGGGGCa                                            >  2:4323117/1‑93 (MQ=255)
                                                 aaaaTTTATCCGGTGAATGTGGTCGAAAAACAAAGAGGAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAg                                        <  2:2779901/93‑1 (MQ=255)
                                                    aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGGCCCCCCGGATAGCGGGCGGCCGAAggga                                     >  1:3164837/1‑91 (MQ=255)
                                                    aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGaa                                     >  1:546935/1‑93 (MQ=255)
                                                    aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAGAGCGGGCGGCAGAAggga                                     >  2:1973561/1‑91 (MQ=255)
                                                     tttATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCCCCCCGGATAGCGGGGGGCAAAAGGaa                                     >  2:3626653/1‑93 (MQ=255)
                                                     tttATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                    >  2:3777256/1‑93 (MQ=255)
                                                      ttATCCGGTGAAGGTGGTCGAAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc                                   <  2:4541526/93‑1 (MQ=255)
                                                       tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                  >  2:511668/1‑93 (MQ=255)
                                                       tATCCGGTGAATGCGGTCGGAAAACAAAGAGGAAAGGGGGGGGACTAATCGGCTGGGAAGGGCGGCCCGGTAAACGCGCGGCAGAAGGAATCa                                  >  1:4076803/1‑93 (MQ=255)
                                                        agccGGTAAATGTGGTCGGAAAACAGAGAGGAAAGGGGGGAGGCTAATCGGCAGGGAAGGCCGCGCCGGATAGCGGGCGGCAGAAGGAATCAg                                 <  2:4520400/91‑1 (MQ=255)
                                                          ccGGTGAATGTGGTCGGAAAACATAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa                               <  1:2179612/93‑1 (MQ=255)
                                                          ccGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATTAGaa                               >  2:3145344/1‑93 (MQ=255)
                                                             gTGAATGTGGTTGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                           >  2:2369099/1‑93 (MQ=255)
                                                             gTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGCCCTCCCAGGATAGCGGGGGGCAGAAGGAATCAGAAttt                            >  2:3039412/1‑93 (MQ=255)
                                                                aaTGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCAGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                        >  2:4392656/1‑93 (MQ=255)
                                                                 aTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                        >  1:671600/1‑93 (MQ=255)
                                                                 aTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                        >  1:671601/1‑93 (MQ=255)
                                                                  tgtgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAgg                       >  1:444420/1‑93 (MQ=255)
                                                                    tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGGTAGCGGGCGGCCGAAGGAATCAGAATTTCCAGGTc                     >  2:1087522/1‑93 (MQ=255)
                                                                    gggTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc                     <  1:3851325/92‑1 (MQ=255)
                                                                        cgaaaaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                 <  1:88267/90‑1 (MQ=255)
                                                                        cGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa                  <  1:767816/93‑1 (MQ=255)
                                                                        cGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                 <  1:2369583/93‑1 (MQ=255)
                                                                            aaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc             <  1:757603/93‑1 (MQ=255)
                                                                            aaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc             <  2:440844/93‑1 (MQ=255)
                                                                             aagaaaGAGGAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt            <  2:169253/90‑1 (MQ=255)
                                                                             aaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCt            >  1:702150/1‑93 (MQ=255)
                                                                              aCAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc          >  1:4125357/1‑93 (MQ=255)
                                                                              aCAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc          >  1:4663802/1‑93 (MQ=255)
                                                                               cAAAGAGGAAAGGGGGGGGTGCCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGTCGGCAGAAGGAATCAGAATTTCCAGGGCAGACGGGCt            <  2:2344767/93‑1 (MQ=255)
                                                                                aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc          >  1:4103641/1‑93 (MQ=255)
                                                                                 aaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGcca        >  1:1064307/1‑91 (MQ=255)
                                                                                  agagGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt      <  1:406745/93‑1 (MQ=255)
                                                                                  agagGAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg       <  1:876157/93‑1 (MQ=255)
                                                                                    agGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt     <  1:2480437/93‑1 (MQ=255)
                                                                                     ggAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg    >  2:474173/1‑93 (MQ=255)
                                                                                      gAAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg    <  1:3429897/93‑1 (MQ=255)
                                                                                        aaGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa  <  1:1593768/93‑1 (MQ=255)
                                                                                          |                                                                                           
ATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCA  >  NC_000913/379147‑379327

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