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.2% 1173.6 / ‑2.3 370A369A (GCA→GCGmetEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (188/179);  total (189/181)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 399 total aligned reads displayed.

TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838
                                                                                          |                                                                                         
 gCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGct                                                                                        >  2:6345197/1‑93 (MQ=255)
   aGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctct                                                                                      <  1:901422/93‑1 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGGTCTg                                                                                     >  1:4114718/1‑93 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTg                                                                                     <  1:104101/93‑1 (MQ=255)
    gAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTg                                                                                     <  2:7046706/93‑1 (MQ=255)
       gTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCa                                                                                  <  2:5036352/93‑1 (MQ=255)
         gAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCaga                                                                                <  2:7498589/93‑1 (MQ=255)
          aaaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCagag                                                                               <  2:8231529/93‑1 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCagagg                                                                              >  2:2434129/1‑92 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCagact                                                                              >  1:1228663/1‑91 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              >  1:9560914/1‑93 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              <  1:499883/93‑1 (MQ=255)
            aaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTg                                                                             >  1:9624068/1‑93 (MQ=255)
               cTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAg                                                                          <  1:2166964/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         <  1:3975223/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         >  1:9965701/1‑93 (MQ=255)
                 ggTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgcg                                                                        <  1:7731089/93‑1 (MQ=255)
                 ggTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgcg                                                                        >  2:1752801/1‑93 (MQ=255)
                 ggTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgcc                                                                        >  1:7398692/1‑92 (MQ=255)
                   tttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcc                                                                      >  1:4681120/1‑93 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    >  1:8361686/1‑93 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    >  2:6594293/1‑93 (MQ=255)
                     tGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGcccc                                                                    <  2:8664798/93‑1 (MQ=255)
                       ccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGa                                                                  <  1:7117424/93‑1 (MQ=255)
                       ccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGa                                                                  <  2:1207916/93‑1 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                <  2:5582276/93‑1 (MQ=255)
                          tggCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATACGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTc                                                               <  2:89789/91‑1 (MQ=255)
                          tCGCCCTACAAAAATGCCATGAGCTGGCACTGCTCCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTc                                                               <  1:7891385/93‑1 (MQ=255)
                          tCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTc                                                               >  2:8579284/1‑93 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGGTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              <  1:847504/93‑1 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              >  2:4716333/1‑93 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              >  2:8477752/1‑93 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              <  1:3630744/93‑1 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             <  2:1425161/93‑1 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             >  2:7950989/1‑93 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             >  2:2165556/1‑93 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             <  1:3023644/93‑1 (MQ=255)
                             cccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGCCACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAgg                                                            <  2:1488971/93‑1 (MQ=255)
                             cccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAgg                                                            >  1:279515/1‑93 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          >  1:5343428/1‑93 (MQ=255)
                               cTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCa                                                          >  1:9379728/1‑93 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAc                                                         >  2:2579643/1‑93 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGGGGAGCGCCCCGATTCAGGCAc                                                         >  1:3148115/1‑93 (MQ=255)
                                  cAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgt                                                       <  2:4317544/93‑1 (MQ=255)
                                  cAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgt                                                       <  2:6464553/93‑1 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      <  1:2855892/93‑1 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      >  1:5591476/1‑93 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcgt                                                    >  1:3997967/1‑93 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcgt                                                    >  1:4281785/1‑93 (MQ=255)
                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACgtcgtc                                                   >  1:1825402/1‑93 (MQ=255)
                                       aTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGGACGTCGGCa                                                  >  1:6909478/1‑93 (MQ=255)
                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAc                                                 <  1:9146259/93‑1 (MQ=255)
                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAc                                                 <  1:6335691/93‑1 (MQ=255)
                                         gCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAct                                                <  2:2980106/93‑1 (MQ=255)
                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActc                                               <  2:8383009/93‑1 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGGGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              >  2:1021129/1‑93 (MQ=255)
                                            aTGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTa                                             >  1:6115181/1‑93 (MQ=255)
                                            aTGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTa                                             <  2:309242/93‑1 (MQ=255)
                                            aTGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTa                                             >  2:38509/1‑93 (MQ=255)
                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            >  2:6083546/1‑93 (MQ=255)
                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            >  1:8663078/1‑93 (MQ=255)
                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            >  2:3564794/1‑93 (MQ=255)
                                                aCTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcg                                         >  2:4231237/1‑93 (MQ=255)
                                                 cTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgc                                        >  2:4017883/1‑93 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGCTTCAGGCACGTCGTCACTCTACCcgcg                                       >  1:5967037/1‑93 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGt                                      <  2:1921506/93‑1 (MQ=255)
                                                    gCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGCGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGtc                                     >  1:1146599/1‑92 (MQ=255)
                                                    gCACTGCTGCGCGATGCGCTGAACAGGGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTa                                     <  2:7607988/93‑1 (MQ=255)
                                                    gCACTGCCGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTa                                     >  2:8145777/1‑93 (MQ=255)
                                                     cACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTAc                                    >  2:9301528/1‑93 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   >  2:6250380/1‑93 (MQ=255)
                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                <  2:1678566/93‑1 (MQ=255)
                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                >  1:5275544/1‑93 (MQ=255)
                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               >  1:8554934/1‑93 (MQ=255)
                                                           tgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATc                              >  1:9852698/1‑93 (MQ=255)
                                                            gcggGAGGCGCTGAACGGTGGTGACACGGCGGCTGTGGCAGAGGGGGGCCCCCCGATTCGGGCCCGTCGCCCCTCTCCCCGCGTACATAATcc                             <  2:3367569/93‑1 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             >  1:9235082/1‑93 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGACGGCTCTGGCAGAGTGGAGCGCGCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             <  1:1276676/93‑1 (MQ=255)
                                                                 gtGCGCTGAACGGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATGCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcgg                        <  2:1307788/92‑1 (MQ=255)
                                                                    cgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAg                     <  1:8831118/93‑1 (MQ=255)
                                                                      cTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaa                   >  2:9390876/1‑93 (MQ=255)
                                                                      cTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaa                   >  2:9778530/1‑93 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  <  1:8042027/93‑1 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATACGGCGGTAGaaaa                 >  2:5176749/1‑93 (MQ=255)
                                                                         aaCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAg                <  1:5162527/93‑1 (MQ=255)
                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTTCATAATCCGGCGGTAGAAAAGCg              >  1:8157244/1‑93 (MQ=255)
                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCg              >  1:1549085/1‑93 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            >  2:9604937/1‑93 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  1:461357/93‑1 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGGACATAATCCGGCGGTAGAAAAGCGTc            >  1:4505846/1‑93 (MQ=255)
                                                                              tggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGtat           <  2:5802455/93‑3 (MQ=255)
                                                                               ggtgACACGGCGGCTCTGGCATAGTGGAGCGCCCCTATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTg          >  1:1762999/1‑93 (MQ=255)
                                                                                 tgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        <  1:2350819/93‑1 (MQ=255)
                                                                                 tgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        <  1:10030082/93‑1 (MQ=255)
                                                                                 tgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggc        <  2:2729850/93‑1 (MQ=255)
                                                                                  gACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcg       >  2:4504085/1‑93 (MQ=255)
                                                                                   acacGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcgg      >  1:6467862/1‑93 (MQ=255)
                                                                                     acGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGgcggcg    >  1:32714/1‑93 (MQ=255)
                                                                                      cggcggCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTCCCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGa   >  1:435279/1‑93 (MQ=255)
                                                                                       ggcggcTCTGGGAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAt  <  1:9658361/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838

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