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,7490AG98.3% 548.0 / ‑1.2 177A369A (GCA→GCGmetEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (76/98);  total (78/99)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 193 total aligned reads displayed.

TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838
                                                                                          |                                                                                         
tGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACggcggc                                                                                         <  1:4001797/93‑1 (MQ=255)
 gCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGct                                                                                        >  2:3054039/1‑93 (MQ=255)
   aGAAGTGAAAAGCAGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGctct                                                                                      >  1:243594/1‑93 (MQ=255)
      aGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGc                                                                                   <  2:3955306/93‑1 (MQ=255)
       gTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCa                                                                                  >  1:2968389/1‑93 (MQ=255)
        tGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCag                                                                                 <  1:1108877/93‑1 (MQ=255)
         gAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCaga                                                                                <  1:1287537/93‑1 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              <  1:3819161/93‑1 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              <  1:2475405/93‑1 (MQ=255)
           aaaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGt                                                                              >  1:4417092/1‑93 (MQ=255)
            aaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTg                                                                             <  2:503914/93‑1 (MQ=255)
            aaGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTg                                                                             <  2:2647105/93‑1 (MQ=255)
               cTGGTTTGCCTTCGCCCTACAAAAGTGCCATGAGCTGGCACTGCTGCGCGATGCTCTGAACAGTGGTGACACGGCGGCTCTGGCGGAGTGGAg                                                                          <  2:4207828/93‑1 (MQ=255)
               cTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAg                                                                          <  1:291360/93‑1 (MQ=255)
                tGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAgc                                                                         <  2:1263666/93‑1 (MQ=255)
                 ggTTTGCCTTCGCGCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGTTGACACGGCGGCTCTGGCAGAGTGGAgcg                                                                        <  2:4081514/93‑1 (MQ=255)
                    ttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTTGAGCGccc                                                                     <  2:3997675/93‑1 (MQ=255)
                    ttGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGccc                                                                     <  2:616699/93‑1 (MQ=255)
                      tccTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCg                                                                   <  2:4376786/92‑1 (MQ=255)
                       ccTTCGCCCTCCAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGa                                                                  <  1:1994169/93‑1 (MQ=255)
                        cTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAt                                                                 <  1:349201/93‑1 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  1:3026376/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  1:4178088/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  2:3206200/1‑93 (MQ=255)
                         ttCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGAtt                                                                >  2:1610680/1‑93 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              <  2:2928597/93‑1 (MQ=255)
                           cGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCa                                                              <  2:4586463/93‑1 (MQ=255)
                            gCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGATCTGGCAGAGTGGAGCGCCCCGATTCAg                                                             >  2:1529198/1‑93 (MQ=255)
                              ccTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGc                                                           >  2:3519465/1‑93 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAc                                                         >  2:4159047/1‑93 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAc                                                         >  2:2609876/1‑93 (MQ=255)
                                tACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAc                                                         <  1:918235/93‑1 (MQ=255)
                                 aCAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcg                                                        >  2:700115/1‑93 (MQ=255)
                                 aCAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcg                                                        >  1:3390735/1‑93 (MQ=255)
                                  cAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgt                                                       >  1:3265830/1‑93 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      >  2:3178962/1‑93 (MQ=255)
                                   aaaaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtc                                                      >  1:4091767/1‑93 (MQ=255)
                                    agaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcg                                                     >  2:3291541/3‑93 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcgt                                                    <  1:2094788/93‑1 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCAcgtcgt                                                    <  2:1017982/93‑1 (MQ=255)
                                     aaaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAATGGAGCGCCCCGATTCAGGCAcgtcgt                                                    <  1:794000/93‑1 (MQ=255)
                                      aaTGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACgtcgtc                                                   <  1:998669/93‑1 (MQ=255)
                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAc                                                 >  2:1310921/1‑93 (MQ=255)
                                        tGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCAc                                                 >  2:3660195/1‑93 (MQ=255)
                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActc                                               <  2:2231391/93‑1 (MQ=255)
                                          ccATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActc                                               <  1:1079779/93‑1 (MQ=255)
                                           cATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCActct                                              >  2:940450/1‑93 (MQ=255)
                                             tGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            <  1:627791/93‑1 (MQ=255)
                                             tGAACGGGCACTGCTGCGCGATGCGGTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAc                                            <  2:605400/93‑1 (MQ=255)
                                              gAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTAcc                                           <  2:1793679/93‑1 (MQ=255)
                                               aaCTGGCACTGCTGCGCGATGCGATGAACAGTGGTGACACGGCGGCTCTGGCAGTGTGGAGCGCCCCGATTCTGGCACGTCGTCACTCTAccg                                          >  2:239321/1‑92 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       >  2:1153666/1‑93 (MQ=255)
                                                  tGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCcgcg                                       >  1:3896449/1‑93 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGt                                      <  2:4049977/93‑1 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGt                                      <  1:1726551/93‑1 (MQ=255)
                                                   ggCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGt                                      <  2:1176425/93‑1 (MQ=255)
                                                    gCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTa                                     <  2:3810798/93‑1 (MQ=255)
                                                      aCTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACa                                   >  2:1070311/1‑93 (MQ=255)
                                                       ctgctgCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACAt                                  <  2:3399980/93‑1 (MQ=255)
                                                       ctgctgCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACAt                                  >  1:2930586/1‑93 (MQ=255)
                                                       ctgctgCGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACAt                                  >  2:1718626/1‑93 (MQ=255)
                                                        tgctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATa                                 <  1:2723566/93‑1 (MQ=255)
                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                >  1:3911130/1‑93 (MQ=255)
                                                         gctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATaa                                <  2:1342819/93‑1 (MQ=255)
                                                          ctgcGCGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAAt                               <  2:1936433/93‑1 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             <  1:3084677/93‑1 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             >  2:2276322/1‑93 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             <  1:3601992/93‑1 (MQ=255)
                                                            gcgcGATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATcc                             <  2:2533487/93‑1 (MQ=255)
                                                             cgcgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcg                            <  1:3933077/93‑1 (MQ=255)
                                                               cgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggc                          <  1:103271/93‑1 (MQ=255)
                                                               cgATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggc                          <  1:4687448/93‑1 (MQ=255)
                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         >  2:2828255/1‑93 (MQ=255)
                                                                gATGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcg                         >  1:3228674/1‑93 (MQ=255)
                                                                 aTGCGCTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCcggcgg                        <  1:1617527/93‑1 (MQ=255)
                                                                  tGCGCTGCACAGTGGTGACACGGCCGCTCTGCCAGAGTCGAGCGCCCAGATTCAGGCACGTCGTCACTCTCCCCGCGTACATAATCCGGCGGt                       <  2:2114469/93‑1 (MQ=255)
                                                                   gcgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTa                      >  2:3689325/1‑93 (MQ=255)
                                                                   gcgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTa                      >  1:2517718/1‑93 (MQ=255)
                                                                   gcgcTGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTa                      <  1:1464603/93‑1 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGCGAGCCCCGCTTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  >  1:1871305/1‑93 (MQ=255)
                                                                       tGAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaa                  <  1:3679973/93‑1 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 <  1:2631993/93‑1 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 <  1:2530933/93‑1 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 >  1:727274/1‑93 (MQ=255)
                                                                        gAACAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCACCGATTCAGGCACGTCGACACTCTACCCGCGTACATAATCCGGCGGTAGaaaa                 <  1:4212655/93‑1 (MQ=255)
                                                                         aaCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAg                >  1:4047379/1‑93 (MQ=255)
                                                                          aCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGc               >  1:547736/1‑93 (MQ=255)
                                                                          aCAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGc               <  1:909257/93‑1 (MQ=255)
                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCg              >  2:4429899/1‑93 (MQ=255)
                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCg              <  2:4715727/93‑1 (MQ=255)
                                                                           cAGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCg              >  1:1170269/1‑93 (MQ=255)
                                                                            aGTGGTGACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGt             <  2:1411639/93‑1 (MQ=255)
                                                                             gtggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTc            <  1:4562792/93‑1 (MQ=255)
                                                                              tggtgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCt           <  1:3295287/93‑1 (MQ=255)
                                                                               gttgACACGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTg          <  2:4547099/91‑1 (MQ=255)
                                                                                 tgACACGGCGGCTCTGGCAGAGTGGGGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGCCTggc        >  2:3550728/1‑93 (MQ=255)
                                                                                   acacGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTggcgg      >  1:2222498/1‑93 (MQ=255)
                                                                                     acGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGgcggcg    <  1:3448062/93‑1 (MQ=255)
                                                                                     acGGCGGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGgcggcg    >  1:3649990/1‑93 (MQ=255)
                                                                                       ggcggcTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAt  >  2:829811/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGCAGAAGTGAAAAGCTGGTTTGCCTTCGCCCTACAAAAATGCCATGAACTGGCACTGCTGCGCGATGCGCTGAACAGTGGTGACACGGCAGCTCTGGCAGAGTGGAGCGCCCCGATTCAGGCACGTCGTCACTCTACCCGCGTACATAATCCGGCGGTAGAAAAGCGTCTGGCGGCGAT  >  NC_000913/4011659‑4011838

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