Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,107,299 A→G V289A (GTG→GCG)  ymdD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,107,2990AG99.3% 943.5 / ‑3.5 295V289A (GTG→GCG) ymdDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (154/139);  total (156/139)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 306 total aligned reads displayed.

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAG  >  NC_000913/1107209‑1107387
                                                                                          |                                                                                        
tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                        <  2:1712286/93‑1 (MQ=255)
tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                        <  1:2214933/93‑1 (MQ=255)
    ccAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCa                                                                                    <  1:3481767/93‑1 (MQ=255)
     cAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAc                                                                                   <  2:516349/93‑1 (MQ=255)
      agatagatAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAcc                                                                                  <  2:8715825/93‑1 (MQ=255)
       gatagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCa                                                                                 <  2:6963427/93‑1 (MQ=255)
         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                               >  1:6243724/1‑93 (MQ=255)
         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                               <  1:5606901/93‑1 (MQ=255)
          agataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATcc                                                                              >  1:5298947/1‑93 (MQ=255)
          agataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATcc                                                                              <  1:5196524/93‑1 (MQ=255)
             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                           <  1:2111798/93‑1 (MQ=255)
             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                           >  1:2825460/1‑93 (MQ=255)
             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                           <  2:5902971/93‑1 (MQ=255)
             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                           <  1:8229348/93‑1 (MQ=255)
              aaaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAg                                                                          >  1:4260624/1‑93 (MQ=255)
              aaaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAg                                                                          <  1:4965309/93‑1 (MQ=255)
               aaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGa                                                                         >  2:6637819/1‑93 (MQ=255)
                aCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGAc                                                                        <  1:7238852/93‑1 (MQ=255)
                  aGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCg                                                                      <  2:1923013/93‑1 (MQ=255)
                   gCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGa                                                                     <  1:4766979/93‑1 (MQ=255)
                     gATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAgg                                                                   <  2:9241654/93‑1 (MQ=255)
                      aTGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGa                                                                  <  2:6480749/93‑1 (MQ=255)
                       tGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAc                                                                 >  2:9444800/1‑93 (MQ=255)
                       tGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAc                                                                 <  1:6693645/93‑1 (MQ=255)
                         ggTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGGGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGCCCa                                                               <  1:6805685/92‑1 (MQ=255)
                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                             >  1:1013625/1‑93 (MQ=255)
                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                             >  2:1425384/1‑93 (MQ=255)
                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                             >  1:7304782/1‑93 (MQ=255)
                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                             <  2:3105380/93‑1 (MQ=255)
                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                             >  1:652513/1‑93 (MQ=255)
                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                             <  1:1645301/93‑1 (MQ=255)
                            tAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATgg                                                            <  2:3135547/93‑1 (MQ=255)
                             aaCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGt                                                           >  2:180831/1‑93 (MQ=255)
                             aaCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGt                                                           <  1:9795077/93‑1 (MQ=255)
                             aaCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGt                                                           >  1:4912211/1‑93 (MQ=255)
                              aCAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTg                                                          >  1:569862/1‑93 (MQ=255)
                               cAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGa                                                         <  2:4754944/93‑1 (MQ=255)
                                aaaaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGAt                                                        <  1:2334059/93‑1 (MQ=255)
                                aaaaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGAt                                                        <  1:9681441/93‑1 (MQ=255)
                                 aaaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATc                                                       >  2:6555813/1‑93 (MQ=255)
                                  aaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCa                                                      >  2:7098066/1‑93 (MQ=255)
                                   aTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCAc                                                     >  1:9375082/1‑93 (MQ=255)
                                        tCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACt                                                <  2:2387499/93‑1 (MQ=255)
                                         cACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTc                                               <  2:4709620/93‑1 (MQ=255)
                                          aCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCg                                              >  1:5070588/1‑93 (MQ=255)
                                          aCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCg                                              <  1:9921408/93‑1 (MQ=255)
                                            ccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGTACCATGGTGATCACCGACTCGGt                                            >  2:5543344/1‑93 (MQ=255)
                                            ccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGt                                            <  2:5831345/93‑1 (MQ=255)
                                            ccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGt                                            <  2:2969981/93‑1 (MQ=255)
                                            ccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGt                                            <  2:6398956/93‑1 (MQ=255)
                                             cgcgCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGtt                                           <  2:10021638/93‑1 (MQ=255)
                                               cgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTc                                         >  2:3404114/1‑93 (MQ=255)
                                               cgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTc                                         <  1:2023653/93‑1 (MQ=255)
                                                 cTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGt                                       >  2:261561/1‑93 (MQ=255)
                                                   gACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTAc                                     <  2:4934371/93‑1 (MQ=255)
                                                   gACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGCACCATGGTGATCAACGACTCGGTTTCGTAc                                     >  2:2174569/1‑93 (MQ=255)
                                                   gACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTAACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTAc                                     >  2:1569574/1‑93 (MQ=255)
                                                     cTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                   >  1:442909/1‑93 (MQ=255)
                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGCCTCGGTTTCGTACATc                                  >  2:3800800/1‑93 (MQ=255)
                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                  >  2:6720942/1‑93 (MQ=255)
                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                  >  2:9490631/1‑93 (MQ=255)
                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                  >  2:7379740/1‑93 (MQ=255)
                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                  >  2:3952106/1‑93 (MQ=255)
                                                       ggAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATcc                                 >  1:5983284/1‑93 (MQ=255)
                                                        gAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCa                                <  1:1281931/93‑1 (MQ=255)
                                                        gAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCa                                <  2:4597691/93‑1 (MQ=255)
                                                        gAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCa                                >  2:5422850/1‑93 (MQ=255)
                                                         aaGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAg                               <  1:2614088/93‑1 (MQ=255)
                                                          aGTTAAGCACACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAgg                              <  1:6979642/93‑1 (MQ=255)
                                                          aGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAgg                              >  2:5779796/1‑93 (MQ=255)
                                                          aGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAgg                              <  1:91996/93‑1 (MQ=255)
                                                           gTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCCTCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCCTACATCCAGGc                             >  1:4172301/1‑93 (MQ=255)
                                                            ttAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                            <  1:1155737/93‑1 (MQ=255)
                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCTGGCAt                           <  2:9108484/93‑1 (MQ=255)
                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                           <  2:4100907/93‑1 (MQ=255)
                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                           <  2:3310777/93‑1 (MQ=255)
                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                           <  1:9523863/93‑1 (MQ=255)
                                                              aaGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATc                          >  1:3788274/1‑93 (MQ=255)
                                                                   aaCGGTGGCCAAAGGAGAAGACCGCATTCCCCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                     <  1:5816981/93‑1 (MQ=255)
                                                                   aaCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                     <  2:8108210/93‑1 (MQ=255)
                                                                   aaCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                     <  2:7569357/93‑1 (MQ=255)
                                                                     cGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATGGCCACTg                   >  2:8803741/1‑93 (MQ=255)
                                                                      ggTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGc                  >  2:2079077/1‑93 (MQ=255)
                                                                        tGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCa                <  2:9682033/93‑1 (MQ=255)
                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc            <  1:4048664/93‑1 (MQ=255)
                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc            <  1:9051008/93‑1 (MQ=255)
                                                                             aaaGGAGAAGCCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCAACTGCCATAgcg           <  2:7607582/93‑1 (MQ=255)
                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg           >  2:5080062/1‑93 (MQ=255)
                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg           >  2:7423160/1‑93 (MQ=255)
                                                                              aaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcgc          >  2:2759889/1‑93 (MQ=255)
                                                                                 gagaAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg       <  1:7972298/93‑1 (MQ=255)
                                                                                 gagaAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg       >  1:8303913/1‑93 (MQ=255)
                                                                                 gagaAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGACCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg       <  2:8672300/93‑1 (MQ=255)
                                                                                  agaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGt      >  1:301725/1‑93 (MQ=255)
                                                                                  agaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGt      >  1:7382722/1‑93 (MQ=255)
                                                                                  agaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGt      >  1:7595536/1‑93 (MQ=255)
                                                                                  agaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGt      >  1:8795446/1‑93 (MQ=255)
                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg  >  2:8876058/1‑93 (MQ=255)
                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg  >  2:3576911/1‑93 (MQ=255)
                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg  <  2:4424887/93‑1 (MQ=255)
                                                                                          |                                                                                        
TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAG  >  NC_000913/1107209‑1107387

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