Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,130,739 A→G T182A (ACT→GCT)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,7390AG99.7% 971.0 / ‑5.4 309T182A (ACT→GCT) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (150/158);  total (151/158)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 315 total aligned reads displayed.

TATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATT  >  NC_000913/4130649‑4130828
                                                                                          |                                                                                         
tATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCt                                                                                         >  1:7859400/1‑93 (MQ=255)
 aTCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTc                                                                                        >  2:6796992/1‑93 (MQ=255)
  tCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCa                                                                                       >  1:6199854/1‑93 (MQ=255)
  tCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCa                                                                                       >  2:3286831/1‑93 (MQ=255)
   ccGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAg                                                                                      <  2:60973/93‑1 (MQ=255)
     ggAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGtt                                                                                    <  2:6336135/93‑1 (MQ=255)
      gAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttc                                                                                   >  2:7780497/1‑93 (MQ=255)
      gAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttc                                                                                   >  2:3190575/1‑93 (MQ=255)
       aaGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttct                                                                                  >  2:617461/1‑93 (MQ=255)
        aGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttctt                                                                                 >  2:882341/1‑93 (MQ=255)
         gTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttcttc                                                                                <  1:16342/93‑1 (MQ=255)
         attCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttcttc                                                                                <  1:6467608/92‑1 (MQ=255)
          ttCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttcttct                                                                               <  1:2789942/93‑1 (MQ=255)
          ttCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATATGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttcttct                                                                               >  2:5438555/1‑93 (MQ=255)
           tCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGttcttctt                                                                              <  1:4560191/93‑1 (MQ=255)
             aCCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCg                                                                            <  1:1787554/93‑1 (MQ=255)
             aCCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCg                                                                            <  1:4995333/93‑1 (MQ=255)
              cccGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGa                                                                           <  1:9827/93‑1 (MQ=255)
              cccGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGa                                                                           <  1:7432792/93‑1 (MQ=255)
              cccGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGa                                                                           <  2:1374171/93‑1 (MQ=255)
              cccGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGa                                                                           <  1:5520278/93‑1 (MQ=255)
              cccGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGa                                                                           <  1:3647878/93‑1 (MQ=255)
                cGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATg                                                                         <  1:5422902/93‑1 (MQ=255)
                 ggAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGt                                                                        <  1:67903/93‑1 (MQ=255)
                  gAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTc                                                                       >  2:6712420/1‑93 (MQ=255)
                   aaGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCg                                                                      >  2:763759/1‑93 (MQ=255)
                       aaaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGACAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAg                                                                  >  2:6623625/1‑93 (MQ=255)
                        aaaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGc                                                                 <  2:4130448/93‑1 (MQ=255)
                         aaaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCt                                                                >  2:4783094/1‑93 (MQ=255)
                          aaGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTa                                                               >  1:2158340/1‑93 (MQ=255)
                           aGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTAc                                                              <  2:3470614/93‑1 (MQ=255)
                           aGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTAc                                                              >  2:2119230/1‑93 (MQ=255)
                           aGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTAc                                                              <  2:3632692/93‑1 (MQ=255)
                            gcgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTAcc                                                             <  2:2578491/93‑1 (MQ=255)
                             cgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCt                                                            >  2:7444360/1‑93 (MQ=255)
                             cgcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCt                                                            <  1:3168957/93‑1 (MQ=255)
                              gcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTg                                                           >  1:5966955/1‑93 (MQ=255)
                              gcTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTg                                                           <  2:1573882/93‑1 (MQ=255)
                               cTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGc                                                          <  2:5516269/93‑1 (MQ=255)
                               cTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGc                                                          <  2:6266738/93‑1 (MQ=255)
                               cTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGc                                                          >  2:4976003/1‑93 (MQ=255)
                                 cAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGt                                                        >  1:6967562/1‑93 (MQ=255)
                                  aGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGtt                                                       <  2:4633797/93‑1 (MQ=255)
                                   ggCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGttt                                                      >  2:3804158/1‑93 (MQ=255)
                                   ggCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGttt                                                      <  1:7925756/93‑1 (MQ=255)
                                    gCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGACAGGTACCTGCGtttt                                                     >  1:6912048/1‑93 (MQ=255)
                                    gCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGtttt                                                     >  1:729601/1‑93 (MQ=255)
                                    gCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGtttt                                                     <  2:6517331/93‑1 (MQ=255)
                                      ggATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCg                                                   >  2:3540201/1‑93 (MQ=255)
                                       gATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGt                                                  >  1:3225578/1‑93 (MQ=255)
                                       gATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGt                                                  >  2:1923317/1‑93 (MQ=255)
                                       gATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGt                                                  <  1:2603608/93‑1 (MQ=255)
                                         tttGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGa                                                >  1:2711051/1‑93 (MQ=255)
                                         tttGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGa                                                >  1:5255592/1‑93 (MQ=255)
                                           tGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGAcc                                              >  1:2431535/1‑93 (MQ=255)
                                           tGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGAcc                                              >  2:5897865/1‑93 (MQ=255)
                                           tGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGAcc                                              >  1:4417810/1‑93 (MQ=255)
                                           tGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATCGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGAcc                                              <  1:5894729/93‑1 (MQ=255)
                                              ttAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCt                                           >  2:4861934/1‑93 (MQ=255)
                                                aaTCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCtgt                                         <  1:424354/93‑1 (MQ=255)
                                                aaTCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCtgt                                         <  2:2692576/93‑1 (MQ=255)
                                                aaTCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCtgt                                         <  1:1037901/93‑1 (MQ=255)
                                                  tCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTATCTGCGTTTTCGTGACCGCTgtgt                                       >  2:5677327/1‑93 (MQ=255)
                                                  tCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTgtgt                                       >  2:4198445/1‑93 (MQ=255)
                                                  tCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTgtgt                                       >  1:1280223/1‑93 (MQ=255)
                                                    tGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATCGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTAt                                     >  1:2079406/1‑93 (MQ=255)
                                                     gAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATCGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATc                                    <  2:1159110/93‑1 (MQ=255)
                                                       aaCGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATcgg                                  >  1:3424510/1‑93 (MQ=255)
                                                       aaCGCAAAGTGGATGCCGGAGCCAACCGCGCGATCGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATcgg                                  >  2:1292887/1‑93 (MQ=255)
                                                         cGCAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATcggcg                                >  1:7505724/1‑93 (MQ=255)
                                                           cAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTCCCTGCGTTTTCGTGACCGCTGTGTATCGGCggg                              >  1:7625428/1‑93 (MQ=255)
                                                           cAAAGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCggg                              >  2:7104446/1‑93 (MQ=255)
                                                            aaaGTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTAACCGCTGTGTATCGGCGGGc                             <  1:7121670/93‑1 (MQ=255)
                                                               gTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCgggggtt                          >  2:5134860/1‑89 (MQ=255)
                                                               gTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGAtt                          >  2:2918748/1‑93 (MQ=255)
                                                               gTGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCAtt                          <  1:3836819/93‑1 (MQ=255)
                                                                tGGATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTg                         <  2:6285212/93‑1 (MQ=255)
                                                                 ggATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGa                        <  2:1806464/93‑1 (MQ=255)
                                                                  gATGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCCTTGAt                       >  2:7735925/1‑93 (MQ=255)
                                                                   aTGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCTGCGGGCATTGAtg                      >  2:1329309/1‑93 (MQ=255)
                                                                   aTGCCGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAtg                      <  1:6297260/93‑1 (MQ=255)
                                                                       cGGAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTGCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGa                  <  2:7483248/93‑1 (MQ=255)
                                                                         gAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGGTGTGGaaa                >  1:373491/1‑93 (MQ=255)
                                                                         gAGCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGaaa                <  1:6416007/93‑1 (MQ=255)
                                                                           gCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGaaaat              >  2:2392136/1‑91 (MQ=255)
                                                                           gCCAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAatt              >  2:2811561/1‑93 (MQ=255)
                                                                            ccAACCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAatta             <  2:1653004/93‑1 (MQ=255)
                                                                            caaaCCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGGTTGGTGTGGTAatta             >  2:6427532/3‑93 (MQ=255)
                                                                               aCCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTc          >  2:6232814/1‑93 (MQ=255)
                                                                               aCCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTc          >  2:4003507/1‑93 (MQ=255)
                                                                               aCCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTc          <  1:2417913/93‑1 (MQ=255)
                                                                               aCCGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTc          >  1:5135734/1‑93 (MQ=255)
                                                                                ccGCGCGATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTcc         <  1:144589/93‑1 (MQ=255)
                                                                                  gcgcgATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTTCTGTGGAAATTATTTCgg       >  2:5204253/1‑93 (MQ=255)
                                                                                  gcgcgATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGGAATTATTCCgg       >  1:1955517/1‑93 (MQ=255)
                                                                                  gcgcgATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCgg       <  1:112600/93‑1 (MQ=255)
                                                                                  gcgcgATCGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCgg       <  1:7045214/93‑1 (MQ=255)
                                                                                    gcgATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGa     >  1:4744647/1‑93 (MQ=255)
                                                                                     cgATTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGaa    <  2:6953479/93‑1 (MQ=255)
                                                                                       aTTGCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGGGGAAATTATTTCTGGGAtt  >  2:4543604/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATCCGGAAGTTCACCCGGAAGCAAAAAGCGCTCAGGCGGATTTGCTTAATCTGAAACGCAAAGTGGATGCCGGAGCCAACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATT  >  NC_000913/4130649‑4130828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑