Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,011,749 A→G A369A (GCA→GCG metE → tetrahydropteroyltriglutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,7490AG99.6% 817.2 / NA 257A369A (GCA→GCGmetEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (131/125);  minor base C (1/0);  total (132/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838
                                                                                          |                                                                                         
tGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACggcggc                                                                                         <  2:4489458/93‑1 (MQ=255)
  cAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctc                                                                                       >  1:789020/1‑93 (MQ=255)
   aGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctct                                                                                      <  1:6112120/93‑1 (MQ=255)
   aGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctct                                                                                      >  1:3061727/1‑93 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTg                                                                                     >  1:3324129/1‑93 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTg                                                                                     >  1:3980604/1‑93 (MQ=255)
         gAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCaga                                                                                <  2:1763350/93‑1 (MQ=255)
          aaaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCagag                                                                               >  1:4884726/1‑93 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              <  1:5717068/93‑1 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              >  1:4957727/1‑93 (MQ=255)
            aaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTg                                                                             <  1:2250297/93‑1 (MQ=255)
             aGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTgg                                                                            <  2:1286670/93‑1 (MQ=255)
             aGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTgg                                                                            >  1:5606562/1‑93 (MQ=255)
             aGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTAGTGACACGGCGGCTCTGGCAGAGTgg                                                                            >  2:5093969/1‑93 (MQ=255)
              gCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGa                                                                           >  1:1576235/1‑93 (MQ=255)
              gCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGa                                                                           <  1:4658244/93‑1 (MQ=255)
              gCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGa                                                                           <  2:3985918/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         >  1:5886137/1‑93 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         >  2:3923597/1‑93 (MQ=255)
                 ggTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGATCAGTGGTGACACGGCGGCTCTGGCAGAGAGGAgcg                                                                        >  2:1842146/1‑93 (MQ=255)
                  gTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgcgc                                                                       >  2:5158768/1‑93 (MQ=255)
                    ttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGccc                                                                     <  2:3663244/93‑1 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    >  1:2657260/1‑93 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    <  2:2993586/93‑1 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    >  2:4813598/1‑93 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    <  2:2831838/93‑1 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCGCGAtt                                                                >  1:3808276/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  1:315315/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  2:4136853/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  2:1931008/1‑93 (MQ=255)
                          tCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTc                                                               <  2:768094/93‑1 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              <  2:2300935/93‑1 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGCTGCGCTGAACAGTTGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             <  2:6098536/93‑1 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             <  1:1374571/93‑1 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             >  2:5569222/1‑93 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             >  1:2429112/1‑93 (MQ=255)
                             cccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAgg                                                            >  1:2468972/1‑93 (MQ=255)
                              ccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGc                                                           <  1:1236682/93‑1 (MQ=255)
                              ccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGc                                                           <  2:2277662/93‑1 (MQ=255)
                               cTACAAAAATGCCATGAACTGTCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          <  2:4766841/93‑1 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          >  2:4763393/1‑93 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          >  1:5492146/1‑93 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAc                                                         <  1:5684174/93‑1 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAc                                                         >  1:1030146/1‑93 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      <  2:2031167/93‑1 (MQ=255)
                                    aaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGCGCGCCCCGATTCAGGCAcgtcg                                                     >  2:1893916/1‑93 (MQ=255)
                                    aaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcg                                                     >  1:3584008/1‑93 (MQ=255)
                                    aaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcg                                                     >  1:6221107/1‑93 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcgt                                                    >  2:4060737/1‑93 (MQ=255)
                                      aaTGCCATGAGCTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACgtcgtc                                                   <  1:4322069/93‑1 (MQ=255)
                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACgtcgtc                                                   <  2:2544386/93‑1 (MQ=255)
                                         gCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAct                                                >  2:3507987/1‑93 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              >  1:4565923/1‑93 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              <  2:5324433/93‑1 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              <  2:6291750/93‑1 (MQ=255)
                                            aTGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTa                                             <  1:1937219/93‑1 (MQ=255)
                                            aTGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTa                                             >  1:1722243/1‑93 (MQ=255)
                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            >  1:851513/1‑93 (MQ=255)
                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            >  2:5593746/1‑93 (MQ=255)
                                                 cTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCGGATTCAGGCACGTCGTCACTCTACCcgc                                        <  1:1091479/93‑1 (MQ=255)
                                                 cTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgc                                        <  2:1409129/93‑1 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       <  2:2472410/93‑1 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       >  2:2197485/1‑93 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       >  1:2554280/1‑93 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGGCACTCTACCCGCGt                                      >  2:1001146/1‑93 (MQ=255)
                                                    gCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTa                                     <  1:1655938/93‑1 (MQ=255)
                                                     cACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTAc                                    <  1:2088244/93‑1 (MQ=255)
                                                      aCTGCTGCGGGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   <  2:1876529/93‑1 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGGACGTCGTCACTCTACCCGCGGACa                                   >  2:85531/1‑93 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   >  1:2390681/1‑93 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   <  2:4073222/93‑1 (MQ=255)
                                                       ctgctgCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACAt                                  <  2:294145/93‑1 (MQ=255)
                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                <  2:1343584/93‑1 (MQ=255)
                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               <  2:4410235/93‑1 (MQ=255)
                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               <  1:1449859/93‑1 (MQ=255)
                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               >  1:904526/1‑93 (MQ=255)
                                                           tgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATc                              <  2:421432/93‑1 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             <  2:2113663/93‑1 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             >  1:4823802/1‑93 (MQ=255)
                                                             cgcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcg                            <  2:3365840/93‑1 (MQ=255)
                                                              gcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcgg                           <  2:4452109/93‑1 (MQ=255)
                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         <  2:5092024/93‑1 (MQ=255)
                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         >  1:4503907/1‑93 (MQ=255)
                                                                  tGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGt                       >  1:4055587/1‑93 (MQ=255)
                                                                    cgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAg                     >  2:5141110/1‑93 (MQ=255)
                                                                      cTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaa                   <  2:3789462/93‑1 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  >  1:4423714/1‑93 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGGAGaaaa                 >  1:54138/1‑93 (MQ=255)
                                                                         aaCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAg                >  2:2144470/1‑93 (MQ=255)
                                                                         aaCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAg                <  2:925683/93‑1 (MQ=255)
                                                                          aCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGc               <  2:3315506/93‑1 (MQ=255)
                                                                          aCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGc               <  2:3864877/93‑1 (MQ=255)
                                                                            aGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGt             <  1:5370977/93‑1 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  1:2900629/93‑1 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  2:5265719/93‑1 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  1:2681630/93‑1 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  2:2734068/93‑1 (MQ=255)
                                                                                  gACACGGCGGCTCTGGGAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcg       <  2:3302468/93‑1 (MQ=255)
                                                                                     acGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcggca    >  2:2622586/1‑92 (MQ=255)
                                                                                       ggcggcTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGGGGCGAt  >  2:1511141/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838

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