Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 939,379 G→A 6.0% intergenic (+42/‑49) rarA → / → serS recombination factor/serine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913939,3790GA 6.0% 264.1 / 4.3 84intergenic (+42/‑49)rarA/serSrecombination factor/serine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base A (2/3);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGA  >  NC_000913/939235‑939528
                                                                                                                                                |                                                                                                                                                     
cTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCACAGGc                                                                                                                                                  >  2:29034/1‑150 (MQ=255)
     attTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTgg                                                                                                                                             >  1:467456/1‑150 (MQ=255)
        tCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCCTGCGACca                                                                                                                                          >  1:249074/1‑150 (MQ=255)
         cccGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCcac                                                                                                                                         >  2:1064881/1‑150 (MQ=255)
          ccGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCcaca                                                                                                                                        >  1:71514/1‑150 (MQ=255)
          ccGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCcaca                                                                                                                                        <  2:1271864/150‑1 (MQ=255)
           cGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACAt                                                                                                                                       <  1:1293461/150‑1 (MQ=255)
            gACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCCTGCGACcacagg                                                                                                                                      <  2:249074/150‑3 (MQ=255)
            gACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGCCACAGCATGCGACcacagg                                                                                                                                      <  2:153106/150‑3 (MQ=255)
              cAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACAtctc                                                                                                                                    >  2:954222/1‑150 (MQ=255)
                aaCAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTccc                                                                                                                                  >  2:478866/1‑150 (MQ=255)
                  cAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCAt                                                                                                                                <  1:1051001/150‑1 (MQ=255)
                    ggggCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCAttt                                                                                                                              >  1:1089143/1‑150 (MQ=255)
                        ccTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAAtt                                                                                                                          >  1:376777/1‑150 (MQ=255)
                         cTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTc                                                                                                                         >  2:553141/1‑150 (MQ=255)
                         cTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTc                                                                                                                         <  2:1089143/150‑1 (MQ=255)
                            gAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCAtt                                                                                                                               >  2:1041014/1‑141 (MQ=255)
                            gAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCAtt                                                                                                                               <  1:1041014/141‑1 (MQ=255)
                                 cAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGca                                                                                                                 >  1:379239/1‑150 (MQ=255)
                                   aGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGcaca                                                                                                               >  2:1299064/1‑150 (MQ=255)
                                   aGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGcaca                                                                                                               <  2:941420/150‑1 (MQ=255)
                                    gATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAg                                                                                                              <  2:828920/150‑1 (MQ=255)
                                      ttGGCGAGAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGa                                                                                                            >  1:1076272/1‑150 (MQ=255)
                                        ggCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATa                                                                                                          <  2:358408/150‑1 (MQ=255)
                                          cGAGAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAg                                                                                                        <  2:1076272/150‑1 (MQ=255)
                                           gAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGc                                                                                                       <  2:69460/150‑1 (MQ=255)
                                             aaaGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCAt                                                                                                     <  2:71514/150‑1 (MQ=255)
                                               aGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGc                                                                                                   >  1:601711/1‑150 (MQ=255)
                                               aGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGCCACATCTcc                                                                                                                                   <  2:196837/117‑1 (MQ=255)
                                               aGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGCCACATCTcc                                                                                                                                   >  1:196837/1‑117 (MQ=255)
                                                gCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCt                                                                                                  >  2:387671/1‑150 (MQ=255)
                                                     ccTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCAGAGGATAAGCATGCTCGATc                                                                                             <  2:1163129/150‑1 (MQ=255)
                                                     ccTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATc                                                                                             >  1:109715/1‑150 (MQ=255)
                                                      ctggctggCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATcc                                                                                            <  2:1067747/150‑1 (MQ=255)
                                                             gctgAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATctg                                                                                     >  2:972905/1‑150 (MQ=255)
                                                                 aaCAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCtgctgc                                                                                 >  1:528738/1‑150 (MQ=255)
                                                                  aCAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGTCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCg                                                                                <  2:509113/150‑1 (MQ=255)
                                                                    aGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTa                                                                              <  2:379239/150‑1 (MQ=255)
                                                                     ggATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTaa                                                                             >  1:699357/1‑150 (MQ=255)
                                                                       aTCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTTCGGTAATGTTGTTACTGTATCCCTGTGGTCTCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATg                                                                           <  2:1437893/150‑1 (MQ=255)
                                                                       aTCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATg                                                                           <  2:109715/150‑1 (MQ=255)
                                                                        tCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGa                                                                          >  1:1118584/1‑150 (MQ=255)
                                                                          aaaaTAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGc                                                                        >  1:741209/1‑150 (MQ=255)
                                                                             aTAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAg                                                                     <  2:1032575/150‑1 (MQ=255)
                                                                                gCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACg                                                                  <  1:302936/150‑1 (MQ=255)
                                                                                 ccccATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGc                                                                 <  2:924501/150‑1 (MQ=255)
                                                                                   ccATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAg                                                               >  2:489254/1‑150 (MQ=255)
                                                                                        aaaCGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCt                                                          >  2:1295898/1‑150 (MQ=255)
                                                                                          aCGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGa                                                        <  1:886326/150‑1 (MQ=255)
                                                                                           cGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGaa                                                       >  1:1255479/1‑150 (MQ=255)
                                                                                             cTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCTCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGaaaa                                                     <  2:793780/150‑1 (MQ=255)
                                                                                                 cGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTg                                                 >  1:182531/1‑150 (MQ=255)
                                                                                                   ttAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGc                                               <  1:523608/150‑1 (MQ=255)
                                                                                                    tAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCACAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCa                                              <  1:29034/150‑1 (MQ=255)
                                                                                                       tGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGc                                           >  1:1385876/1‑150 (MQ=255)
                                                                                                           aTCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCggg                                       <  2:817635/150‑1 (MQ=255)
                                                                                                             cGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGc                                     >  2:1391715/1‑150 (MQ=255)
                                                                                                                ttctgTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAACCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGTGCttt                                  <  2:525802/146‑1 (MQ=255)
                                                                                                                tGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCttt                                  <  1:114145/150‑1 (MQ=255)
                                                                                                                 gCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTa                                 <  2:182531/150‑1 (MQ=255)
                                                                                                                   gttAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAg                               <  1:923394/148‑1 (MQ=255)
                                                                                                                        tgttgtTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGa                          >  2:195532/1‑150 (MQ=255)
                                                                                                                        tgttgtTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGa                          <  2:817727/150‑1 (MQ=255)
                                                                                                                        tgttgtTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGa                          >  2:726513/1‑150 (MQ=255)
                                                                                                                        tgttgtTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGa                          <  1:489254/150‑1 (MQ=255)
                                                                                                                          ttgttACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATg                        <  1:553141/150‑1 (MQ=255)
                                                                                                                           tgttACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCTATCCGAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGt                       >  2:155917/1‑150 (MQ=255)
                                                                                                                           tgttACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGt                       <  1:1299064/150‑1 (MQ=255)
                                                                                                                             ttACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAg                     <  1:1295898/150‑1 (MQ=255)
                                                                                                                              tACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGa                    >  2:1438551/1‑150 (MQ=255)
                                                                                                                               aCTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGAt                   >  2:790338/1‑150 (MQ=255)
                                                                                                                                  gTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAg                >  1:812771/1‑150 (MQ=255)
                                                                                                                                    aTCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCt              >  1:1210008/1‑150 (MQ=255)
                                                                                                                                     tCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTg             >  2:566968/1‑150 (MQ=255)
                                                                                                                                      cccTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTgg            >  2:1537887/1‑150 (MQ=255)
                                                                                                                                        cTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGgc          <  1:478866/150‑1 (MQ=255)
                                                                                                                                        cTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGgc          <  2:827479/150‑1 (MQ=255)
                                                                                                                                        cTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGgc          <  2:832958/150‑1 (MQ=255)
                                                                                                                                        cTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGgc          <  1:1064881/150‑1 (MQ=255)
                                                                                                                                           tgGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGct       >  2:276805/1‑150 (MQ=255)
                                                                                                                                             gTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGctct     <  1:987337/150‑1 (MQ=255)
                                                                                                                                               cGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTg   >  2:1049842/1‑150 (MQ=255)
                                                                                                                                               cGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTg   >  1:862705/1‑150 (MQ=255)
                                                                                                                                               cGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTg   >  1:1087544/1‑150 (MQ=255)
                                                                                                                                                gCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCttca  <  1:972905/150‑3 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTATTATTTCCCGACAAACAGGGGCCTTGAAGGCAAGATTGGCGAAAAGCTCGCCTGGCTGGCTGAACAGGATCAAAATAGCCCCATAAAACGCTACCGTTAATGTTATCGTTGCGGTAATGTTGTTACTGTATCCCTGTGGTCGCAGGCTGTGGCCACATCTCCCATTTAATTCGATAAGCACAGGATAAGCATGCTCGATCCCAATCTGCTGCGTAATGAGCCAGACGCAGTCGCTGAAAAACTGGCACGCCGGGGCTTTAAGCTGGATGTAGATAAGCTGGGCGCTCTTGA  >  NC_000913/939235‑939528

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: