Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,126 T→A 13.1% intergenic (‑371/+26) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,1260TA13.1% 218.9 / 18.7 90intergenic (‑371/+26)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (38/40);  new base A (4/8);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

TCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAA  >  NC_002947/5078983‑5079270
                                                                                                                                                |                                                                                                                                                
tCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTg                                                                                                                                              <  5:178561/148‑1 (MQ=255)
    tGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCgg                                                                                                                                         <  2:75581/149‑1 (MQ=255)
      cgacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCggg                                                                                                                                        <  2:42831/148‑1 (MQ=255)
      cgacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCggg                                                                                                                                        <  6:128682/148‑1 (MQ=255)
        acgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtt                                                                                                                                     <  3:142639/149‑1 (MQ=255)
        acgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtt                                                                                                                                     >  4:133825/1‑149 (MQ=255)
          gaTTGTTGTGGATGGGAAAGCAGATTTGATTTTTGAGGCCGCTCAGCTGCTCCAAACTACTCCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCtttt                                                                                                                                   <  5:317990/149‑1 (MQ=255)
           aTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                  <  2:22354/149‑1 (MQ=255)
           aTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                  <  6:159396/149‑1 (MQ=255)
            ttgttgTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                  <  8:122873/148‑1 (MQ=255)
               ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGAt                                                                                                                               <  7:118410/148‑1 (MQ=255)
               ttgtGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGAt                                                                                                                               <  7:107960/148‑1 (MQ=255)
                 gtgGATGGGAACGCACATTTGATTTTTGAGGCCGCTCAGAGGCTCCAAACTACACCGCAATCCTTACAGAAGAGCGTAGCTCGACAGGCGCAAGAATTCATTGGCGCACAAACAAACAGCCCGCTCACTAGAGCGGGCTTTTTgtttg                                                                                                                             <  3:281318/148‑5 (MQ=255)
                 gtgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                            >  6:175678/1‑149 (MQ=255)
                 gtgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCTTTTTgtttg                                                                                                                             <  7:182657/148‑5 (MQ=255)
                  tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgggg                                                                                                                           <  3:246130/149‑1 (MQ=255)
                  tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                            <  1:224275/148‑1 (MQ=255)
                   ggATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgggg                                                                                                                           >  7:233200/1‑148 (MQ=255)
                       gggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAAta                                                                                                                                               >  8:221619/1‑123 (MQ=255)
                       gggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAAta                                                                                                                                               <  7:221619/124‑2 (MQ=255)
                        ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGc                                                                                                                     <  7:220813/149‑1 (MQ=255)
                        ggAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGc                                                                                                                     <  3:281321/149‑1 (MQ=255)
                         gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGc                                                                                                                                       >  2:351648/1‑130 (MQ=255)
                         gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGc                                                                                                                                       <  1:351648/130‑1 (MQ=255)
                         gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCg                                                                                                                    >  1:131443/1‑149 (MQ=255)
                         gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCg                                                                                                                    >  1:136449/1‑149 (MQ=255)
                          aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgg                                                                                                                             <  8:109096/139‑1 (MQ=255)
                          aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATgg                                                                                                                             >  7:109096/1‑139 (MQ=255)
                             gCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTg                                                                                                                                 <  6:300911/132‑1 (MQ=255)
                             gCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTg                                                                                                                                 >  5:300911/1‑132 (MQ=255)
                               acaTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCTGGCTTTTTGATGGGGTATGGCGGACt                                                                                                                >  8:142633/1‑147 (MQ=255)
                                 aTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagc                                                                                                            <  1:179925/149‑1 (MQ=255)
                                 aTTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcag                                                                                                             >  2:179925/1‑148 (MQ=255)
                                  tttGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                           <  2:213280/149‑1 (MQ=255)
                                     gATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGa                                                                                                        >  6:242601/1‑149 (MQ=255)
                                      aTTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGa                                                                                                        <  3:52538/148‑1 (MQ=255)
                                       tttttGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                               <  1:263027/108‑1 (MQ=255)
                                       tttttGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAtt                                                                                                                                               >  2:263027/1‑108 (MQ=255)
                                       tttttGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAtt                                                                                                      >  8:2503/1‑149 (MQ=255)
                                                 cGCTCAGCGGCTCCAAACTGCACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGc                                                                                             >  4:321343/1‑148 (MQ=255)
                                                   cTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGAc                                                                                          <  1:323211/149‑1 (MQ=255)
                                                      aGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCg                                                                                                    <  7:292870/136‑1 (MQ=255)
                                                      aGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCg                                                                                                    >  8:292870/1‑136 (MQ=255)
                                                           cTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGt                                                                                  <  5:175678/149‑1 (MQ=255)
                                                           cTCAAAACTACACCTCAATAATTACAGAAGAGCGTATATCGACAGGCTCAATAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGt                                                                                  <  8:233200/149‑1 (MQ=39)
                                                            tCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGtt                                                                                 >  7:281026/1‑149 (MQ=255)
                                                              cAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGAATGGCGGACTCAGCAGGATTCGAACCTGCGAACGCTCGGTTCg                                                                               >  7:297705/1‑149 (MQ=255)
                                                               aaaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGt                                                                              >  4:42444/1‑149 (MQ=255)
                                                                 aCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAg                                                                            >  7:85540/1‑149 (MQ=255)
                                                                   tACACCGCAATCATTACAGAAGAGCGTAGATCGTCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGcc                                                                          >  2:342936/1‑149 (MQ=255)
                                                                   tACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                           <  5:110319/116‑1 (MQ=255)
                                                                   tACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcagca                                                                                                           >  6:110319/1‑116 (MQ=255)
                                                                      acCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAg                                                                       <  3:321343/149‑1 (MQ=255)
                                                                          catacATTACAGAAGAGCGTAGATCGACAGGCTCAAGAAATCATTTGCGCACAAACAAAAAGCCAGGTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAct                                                                   <  8:297705/145‑1 (MQ=37)
                                                                             tCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActct                                                                 >  4:303395/1‑148 (MQ=255)
                                                                             tCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                >  3:151084/1‑149 (MQ=255)
                                                                               aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                              <  5:254819/149‑1 (MQ=255)
                                                                               aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                              <  5:242601/149‑1 (MQ=255)
                                                                                ttACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACGAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                              >  1:14360/1‑148 (MQ=255)
                                                                                    agaagaGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                         <  8:85540/149‑1 (MQ=255)
                                                                                    agaagaGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                         >  5:338520/1‑149 (MQ=255)
                                                                                      aagaGCGTAGATCGTCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                       <  1:342936/149‑1 (MQ=255)
                                                                                      aagaGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                         >  5:381824/1‑147 (MQ=255)
                                                                                        gagCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                                     <  7:2503/149‑1 (MQ=255)
                                                                                           cGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                  >  6:246574/1‑149 (MQ=255)
                                                                                           cGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                  <  8:281026/149‑1 (MQ=255)
                                                                                           cGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                  <  3:42444/149‑1 (MQ=255)
                                                                                                aTCGACAGGCTCAATAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAtt                                                                                                      <  1:241524/92‑1 (MQ=255)
                                                                                                aTCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGAtt                                                                                                      >  2:241524/1‑92 (MQ=255)
                                                                                                  cGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCATTTTGATGGGGTAAGGCGGACTCAGCAGGATTAGAACCTGCGACCCCTCGGTTAGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGtacc                                           >  5:278861/1‑146 (MQ=25)
                                                                                                   gACAGGCTCAAGAATTCATTGGCGCACAAACGAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGt                                          <  2:14360/149‑1 (MQ=255)
                                                                                                    tcAGGATCAAGAATTCATTGGCTCACAAACTAAAAGCCGGCTCTATTGTGTGGGCTTTTTGTTGGGTTATGGCGGACTCAGCAGTATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTg                                         <  6:278861/148‑1 (MQ=18)
                                                                                                    aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTg                                         <  6:76794/149‑1 (MQ=255)
                                                                                                          tCAAGAATTCATTGGCGCACAAACAAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATgg                                     >  1:194439/1‑148 (MQ=255)
                                                                                                            aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGAGATTGTaag                                 >  6:371027/1‑147 (MQ=255)
                                                                                                            aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACAGTGACCACGTCAGATATGAGTCCGTGaacgaagt                                 >  7:129570/1‑142 (MQ=25)
                                                                                                              gAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGttt                               <  3:303395/149‑1 (MQ=255)
                                                                                                              gAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGAGAttggat                                  >  2:109648/1‑141 (MQ=255)
                                                                                                                    aTTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGAc                          >  4:311896/1‑148 (MQ=255)
                                                                                                                    aTTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGAc                          >  3:243828/1‑148 (MQ=255)
                                                                                                                    aTTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGAc                          >  3:168440/1‑148 (MQ=255)
                                                                                                                             ccatcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATg                                                 <  7:69152/112‑1 (MQ=255)
                                                                                                                             ccatcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGCTTCGTAGCCGAGTACTCAATCCAGCTCAGCAATg                                                 >  8:69152/5‑116 (MQ=18)
                                                                                                                             caaacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCa                >  5:133557/1‑149 (MQ=255)
                                                                                                                              catcaaaAAGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCGATGAGTCCGTGATGGTAGtttt                              >  6:144511/4‑134 (MQ=255)
                                                                                                                                      aaGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTa        <  5:246574/149‑1 (MQ=255)
                                                                                                                                        gCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACt      >  7:217002/1‑149 (MQ=255)
                                                                                                                                        gCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTATATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGAGCAGGTTTCCACTGGCTTACt      >  1:195823/1‑149 (MQ=21)
                                                                                                                                        gCCCGCTATATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACt      <  7:142633/149‑1 (MQ=255)
                                                                                                                                            gCTCTATTGAGCGGGCTTTTTGATGGTGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGaaa  <  4:138518/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAA  >  NC_002947/5078983‑5079270

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