Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,186,925 A→G T504A (ACT→GCT)  yqiG → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,186,9250AG99.8% 1373.3 / ‑5.0 393T504A (ACT→GCT) yqiGputative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (210/182);  total (211/182)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 418 total aligned reads displayed.

TATCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTCCACTTCT  >  NC_000913/3186835‑3187014
                                                                                           |                                                                                         
  tCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgcttc                                                                                       >  1:1679106/1‑92 (MQ=255)
  tCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgctta                                                                                       >  2:2390432/1‑93 (MQ=255)
  tCAGAATGAAAAAGAGACTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgctta                                                                                       <  1:6834164/93‑1 (MQ=255)
   cAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTCAAATACATACGTAAGTTTGGCTCGTAAtgcttac                                                                                      <  2:6600950/93‑1 (MQ=255)
    aGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgcttact                                                                                     >  2:6372009/1‑93 (MQ=255)
     gAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTg                                                                                    <  1:8380540/93‑1 (MQ=255)
     gAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTg                                                                                    >  1:4301889/1‑93 (MQ=255)
       aTGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTggg                                                                                  >  2:3031452/1‑93 (MQ=255)
        tGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGa                                                                                 >  2:7713924/1‑93 (MQ=255)
        tGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGa                                                                                 >  1:3056604/1‑93 (MQ=255)
         gAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGAc                                                                                <  2:1594883/93‑1 (MQ=255)
           aaaaGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGc                                                                              <  1:1910862/93‑1 (MQ=255)
           aaaaGAGAGTGATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGc                                                                              <  2:1393439/93‑1 (MQ=255)
            aaaGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGcc                                                                             <  2:2165737/93‑1 (MQ=255)
             aaGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCca                                                                            >  2:7472859/1‑93 (MQ=255)
             aaGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCca                                                                            >  1:2133779/1‑93 (MQ=255)
              agagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcag                                                                           >  2:8395796/1‑93 (MQ=255)
              agagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcag                                                                           >  2:6145903/1‑93 (MQ=255)
              agagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcag                                                                           >  1:3777132/1‑93 (MQ=255)
              agagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcag                                                                           <  2:4818249/93‑1 (MQ=255)
              agagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcag                                                                           <  1:62536/93‑1 (MQ=255)
              agagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcac                                                                           >  1:6003557/1‑92 (MQ=255)
               gagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagc                                                                          <  1:8764055/93‑1 (MQ=255)
                agagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagca                                                                         <  2:2252054/93‑1 (MQ=255)
                agagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagca                                                                         <  1:4757104/93‑1 (MQ=255)
                agagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagca                                                                         <  1:5100458/93‑1 (MQ=255)
                agagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagca                                                                         <  2:3484637/93‑1 (MQ=255)
                   gTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTa                                                                      <  2:3408410/93‑1 (MQ=255)
                   gTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTa                                                                      <  2:200962/93‑1 (MQ=255)
                    ttATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTaa                                                                     <  2:7585417/93‑1 (MQ=255)
                       tatTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgt                                                                  >  2:2900112/1‑93 (MQ=255)
                        atTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgtg                                                                 >  2:8054853/1‑93 (MQ=255)
                         ttGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGa                                                                <  2:3965295/93‑1 (MQ=255)
                         ttGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGa                                                                >  2:7288153/1‑93 (MQ=255)
                          tGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGaa                                                               >  1:4561123/1‑93 (MQ=255)
                            tCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAatt                                                             >  2:4104254/1‑93 (MQ=255)
                            tCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAatt                                                             >  1:6322070/1‑93 (MQ=255)
                              aCGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAattat                                                           <  2:2544202/93‑1 (MQ=255)
                               cGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAattatt                                                          <  1:8672872/93‑1 (MQ=255)
                                gTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTc                                                         >  1:1577163/1‑93 (MQ=255)
                                gTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTc                                                         >  2:8778826/1‑93 (MQ=255)
                                 tttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCa                                                        <  2:5340450/93‑1 (MQ=255)
                                 tttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCa                                                        >  2:1042622/1‑93 (MQ=255)
                                 tttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCa                                                        >  2:2618965/1‑93 (MQ=255)
                                  ttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCAt                                                       <  1:3451121/93‑1 (MQ=255)
                                  ttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCAt                                                       <  2:5457377/93‑1 (MQ=255)
                                   tAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCAtt                                                      >  2:3692845/1‑93 (MQ=255)
                                   tAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCAtt                                                      <  2:2980366/93‑1 (MQ=255)
                                    aaCCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTa                                                     <  2:8179337/93‑1 (MQ=255)
                                     aCCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTAt                                                    >  1:1423084/1‑93 (MQ=255)
                                     aCCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGATCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTAt                                                    >  1:2564386/1‑93 (MQ=255)
                                      ccAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATc                                                   >  2:5655232/1‑93 (MQ=255)
                                       cAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCa                                                  >  1:75889/1‑93 (MQ=255)
                                        aGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCAc                                                 <  1:390184/93‑1 (MQ=255)
                                         gTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACt                                                >  1:7981375/1‑93 (MQ=255)
                                         gTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACt                                                <  2:8554122/93‑1 (MQ=255)
                                          tATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACtt                                               >  2:3874487/1‑93 (MQ=255)
                                           aTCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTa                                              <  1:7088145/93‑1 (MQ=255)
                                           aTCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTa                                              >  1:6116864/1‑93 (MQ=255)
                                           aTCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTa                                              >  2:5972552/1‑93 (MQ=255)
                                           aTCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTa                                              <  1:2400471/93‑1 (MQ=255)
                                             cTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGc                                            >  2:1256371/1‑93 (MQ=255)
                                               tGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCcg                                          <  2:4280871/93‑1 (MQ=255)
                                                 aaTCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCcgcg                                        >  1:2239376/1‑93 (MQ=255)
                                                  aTCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGa                                       <  1:1014717/93‑1 (MQ=255)
                                                   tCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAt                                      >  1:940547/1‑93 (MQ=255)
                                                   tCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAt                                      >  1:5376097/1‑93 (MQ=255)
                                                    cATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAtt                                     <  2:4498616/93‑1 (MQ=255)
                                                    cATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAtt                                     <  1:645059/93‑1 (MQ=255)
                                                     aTTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAttt                                    >  1:5734129/1‑93 (MQ=255)
                                                     aTTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAttt                                    >  2:1005839/1‑93 (MQ=255)
                                                     aTTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAttt                                    >  1:2000275/1‑93 (MQ=255)
                                                        aaGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGa                                 <  2:5247706/93‑1 (MQ=255)
                                                          ggTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGata                               <  2:5135545/93‑1 (MQ=255)
                                                          ggTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGata                               >  1:4473821/1‑93 (MQ=255)
                                                           gTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGatat                              <  1:29303/93‑1 (MQ=255)
                                                           gTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGatat                              >  1:2651614/1‑93 (MQ=255)
                                                           gTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGatat                              >  2:475896/1‑93 (MQ=255)
                                                           gTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGatat                              <  1:4933604/93‑1 (MQ=255)
                                                           gTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGatat                              >  1:5850983/1‑93 (MQ=255)
                                                            ttAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATc                             >  1:7000264/1‑93 (MQ=255)
                                                             tAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCg                            <  2:7853344/93‑1 (MQ=255)
                                                             tAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCg                            <  1:7658238/93‑1 (MQ=255)
                                                              aaaTACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCgg                           >  1:2618151/1‑93 (MQ=255)
                                                               aaTACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCggg                          <  2:5514775/93‑1 (MQ=255)
                                                               aaTACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCggg                          >  2:2743388/1‑93 (MQ=255)
                                                                 tacatacGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGTGcc                        >  2:7432988/1‑93 (MQ=255)
                                                                   catccGTAAGTTTGCTCGGTAAGGCTTCCGGGGACGCCAGCAGTAATTTGAATTTTTCATATCACTTAGCCGCGATTTTGATATCGGGCCAtt                     <  1:7316477/93‑1 (MQ=255)
                                                                      acGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTaa                   >  1:8043791/1‑93 (MQ=255)
                                                                      acGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTaa                   <  2:581718/93‑1 (MQ=255)
                                                                       cGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTaaa                  >  1:621225/1‑93 (MQ=255)
                                                                        gTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTaaaa                 >  1:4056294/1‑93 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc              >  2:6498735/1‑93 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc              >  1:3197291/1‑93 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc              >  1:7505036/1‑93 (MQ=255)
                                                                            gTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACg             >  1:7769780/1‑93 (MQ=255)
                                                                             tttGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGt            >  2:6930454/1‑93 (MQ=255)
                                                                              ttGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtc           <  1:6665626/93‑1 (MQ=255)
                                                                               tGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtct          <  2:1716702/93‑1 (MQ=255)
                                                                                 gCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTcc        >  2:4339229/1‑93 (MQ=255)
                                                                                           |                                                                                         
TATCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTCCACTTCT  >  NC_000913/3186835‑3187014

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