Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,503,660 T→C T149A (ACG→GCG)  b1433 ← putative membrane transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,503,6600TC90.9% 606.3 / 42.4 230T149A (ACG→GCG) b1433putative membrane transport protein
Reads supporting (aligned to +/- strand):  ref base T (9/12);  new base C (108/101);  total (117/113)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 234 total aligned reads displayed.

ACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCGGAATAATGC  >  NC_000913/1503569‑1503752
                                                                                           |                                                                                            
acgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTa                                                                                             <  1:7864314/93‑1 (MQ=255)
 cgatcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCaa                                                                                            >  1:6639304/1‑93 (MQ=255)
   atcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAAt                                                                                          >  1:846058/1‑93 (MQ=255)
   atcCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAAt                                                                                          >  1:5659208/1‑93 (MQ=255)
    tcCCAATAATCATCGCGGCAATTACCACATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAAtt                                                                                         >  2:40558/1‑93 (MQ=255)
      ccAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGa                                                                                       >  1:2234722/1‑93 (MQ=255)
       cAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGAc                                                                                      <  2:7635246/93‑1 (MQ=255)
       cAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGAc                                                                                      <  2:6181323/93‑1 (MQ=255)
       cAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGAc                                                                                      >  1:177966/1‑93 (MQ=255)
        aataatCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGAcc                                                                                     >  2:8852056/1‑93 (MQ=255)
        aataatCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGAcc                                                                                     >  1:306249/1‑93 (MQ=255)
         ataatCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCg                                                                                    >  1:5476931/1‑93 (MQ=255)
          taatCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGt                                                                                   >  2:5706213/1‑93 (MQ=255)
           aatCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTc                                                                                  >  2:8326842/1‑93 (MQ=255)
           aatCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTc                                                                                  <  1:2764806/93‑1 (MQ=255)
             tcatcGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCa                                                                                <  1:9101666/93‑1 (MQ=255)
                tcGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGAc                                                                             >  2:8210356/1‑93 (MQ=255)
                tcGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGAc                                                                             <  2:807774/93‑1 (MQ=255)
                tcGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGAc                                                                             >  2:8809031/1‑93 (MQ=255)
                 cgcgGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACt                                                                            <  2:2052809/93‑1 (MQ=255)
                  gcgGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTg                                                                           >  1:5711405/1‑93 (MQ=255)
                   cgGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTgg                                                                          >  1:1277551/1‑93 (MQ=255)
                    ggCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGc                                                                         <  1:902433/93‑1 (MQ=255)
                    ggCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGc                                                                         <  1:1659107/93‑1 (MQ=255)
                    ggCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGc                                                                         <  2:8496681/93‑1 (MQ=255)
                        aTTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAAc                                                                     >  2:7272244/1‑93 (MQ=255)
                        aTTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAAc                                                                     >  2:2626325/1‑93 (MQ=255)
                        aTTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAAc                                                                     >  2:9601765/1‑93 (MQ=255)
                        aTTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAAc                                                                     >  1:4969037/1‑93 (MQ=255)
                         ttACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACg                                                                    >  2:2107456/1‑93 (MQ=255)
                          tACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGc                                                                   <  1:3735689/93‑1 (MQ=255)
                           aCCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGcc                                                                  >  1:1541633/1‑93 (MQ=255)
                             cGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTg                                                                <  2:7567156/93‑1 (MQ=255)
                                aTAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTaa                                                             >  1:3617576/1‑93 (MQ=255)
                                aTAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTaa                                                             >  2:2872148/1‑93 (MQ=255)
                                 tAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTaaa                                                            <  2:1786179/93‑1 (MQ=255)
                                  aggccggcgCAACGGCCTTGGTTGCCAGCCATTCCAGCAACATCCTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAAc                                                           <  1:1751938/89‑1 (MQ=255)
                                  aGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAAc                                                           >  2:6682675/1‑93 (MQ=255)
                                  aGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAAc                                                           <  1:9997245/93‑1 (MQ=255)
                                   gcgcGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAAcc                                                          <  1:8913864/93‑1 (MQ=255)
                                   gcgcGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAAcc                                                          >  2:3271519/1‑93 (MQ=255)
                                      gggcgCAAGGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCaaa                                                       <  1:2895291/92‑1 (MQ=255)
                                       ggcgCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAg                                                      <  1:607571/93‑1 (MQ=255)
                                        gcgcAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGc                                                     >  1:3546747/1‑93 (MQ=255)
                                        gcgcAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGc                                                     >  1:9724146/1‑93 (MQ=255)
                                         cgcAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCg                                                    >  2:3512311/1‑93 (MQ=255)
                                          gcAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGt                                                   <  2:2896607/93‑1 (MQ=255)
                                             aCGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTCAACCAAAGCGtaat                                                >  1:7031387/1‑93 (MQ=255)
                                             aCGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGtaat                                                >  1:1428317/1‑93 (MQ=255)
                                              cGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGtaata                                               <  2:6778423/93‑1 (MQ=255)
                                                gCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATaaa                                             <  2:3870105/93‑1 (MQ=255)
                                                gCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATaaa                                             >  1:6416750/1‑93 (MQ=255)
                                                 ccTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATaaaa                                            <  2:1486632/93‑1 (MQ=255)
                                                   ttGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATc                                          >  1:8314953/1‑93 (MQ=255)
                                                   ttGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATc                                          >  2:7607648/1‑93 (MQ=255)
                                                   ttGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATc                                          >  1:5939009/1‑93 (MQ=255)
                                                    tGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATcc                                         <  1:9472148/93‑1 (MQ=255)
                                                     ggTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATccc                                        <  2:6345904/93‑1 (MQ=255)
                                                     ggTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATccc                                        <  2:9079443/93‑1 (MQ=255)
                                                      gTTGCCAGCCATACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCg                                       >  2:8636399/1‑93 (MQ=255)
                                                           cagccaTACCAGCAACGTCCTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAg                                  <  1:9869321/93‑1 (MQ=255)
                                                             gccaTACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCa                                >  2:440958/1‑93 (MQ=255)
                                                             gccaTACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCa                                >  1:7431908/1‑93 (MQ=255)
                                                               caTACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCAtt                              >  2:9836672/1‑93 (MQ=255)
                                                               caTACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCAtt                              >  1:8782860/1‑93 (MQ=255)
                                                                atcccAGCAACATACTTCCACACAACGCAAATTGACCGTCGAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTg                             <  1:1759011/90‑1 (MQ=255)
                                                                 tACCAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTgc                            <  2:8779810/93‑1 (MQ=255)
                                                                   ccAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTgccg                          >  2:51974/1‑93 (MQ=255)
                                                                    cAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTgccgc                         <  1:6514339/93‑1 (MQ=255)
                                                                    cAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTgccgc                         >  2:7978016/1‑93 (MQ=255)
                                                                    cAGCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTgccgc                         <  1:4802146/93‑1 (MQ=255)
                                                                      gCAACATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGccgccg                       >  1:7260268/1‑93 (MQ=255)
                                                                        aaCATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCa                     >  2:6057663/1‑93 (MQ=255)
                                                                        aaCATACTTCCACACAACGCAAATTGACCGTCCAGACTGACAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCCCa                     >  2:6098228/1‑93 (MQ=255)
                                                                        aaCATACTTACAGACAACGCAAATTGGCAGTCCAGGCTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCa                     <  1:6101774/93‑1 (MQ=255)
                                                                         aCATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCaa                    >  1:5122251/1‑93 (MQ=255)
                                                                          cATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAg                   <  1:1236199/93‑1 (MQ=255)
                                                                          cATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAg                   >  1:5518073/1‑93 (MQ=255)
                                                                          cATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAg                   >  2:5063972/1‑93 (MQ=255)
                                                                          cATACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAg                   <  1:4309014/93‑1 (MQ=255)
                                                                           aTACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGc                  <  2:6563856/93‑1 (MQ=255)
                                                                            tACTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCg                 <  1:6121376/93‑1 (MQ=255)
                                                                             aCTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGa                <  1:1152138/93‑1 (MQ=255)
                                                                              cTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAg               <  2:581273/93‑1 (MQ=255)
                                                                              cTTCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAg               <  2:140327/93‑1 (MQ=255)
                                                                               ttCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAgg              >  1:9296211/1‑92 (MQ=255)
                                                                               ttCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGt              <  2:5413551/93‑1 (MQ=255)
                                                                               ttCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGt              >  1:4811459/1‑93 (MQ=255)
                                                                                tCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTg             <  1:2431682/93‑1 (MQ=255)
                                                                                tCCACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTg             <  1:4669189/93‑1 (MQ=255)
                                                                                 gcacacaACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAGCCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGc            <  2:1317270/92‑1 (MQ=255)
                                                                                 ccACACAACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGc            >  1:7738683/1‑93 (MQ=255)
                                                                                  cacacaACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCg           <  1:4625940/93‑1 (MQ=255)
                                                                                  cacacaACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCg           >  2:9968825/1‑93 (MQ=255)
                                                                                   acacaACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCgg          >  2:4984232/1‑93 (MQ=255)
                                                                                     acaACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCGGaa        >  1:2159402/1‑93 (MQ=255)
                                                                                      caACGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCGGaat       >  2:4504510/1‑93 (MQ=255)
                                                                                        aCGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCGGaataa     >  2:400242/1‑93 (MQ=255)
                                                                                        aCGCAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGGGCGGaataa     >  2:547711/1‑93 (MQ=255)
                                                                                           tAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCGGAATAATgc  <  1:5389839/93‑1 (MQ=255)
                                                                                           |                                                                                            
ACGATCCCAATAATCATCGCGGCAATTACCGCATAGCGCGGCGCAACGGCCTTGGTTGCCAGCCATACCAGCAACATACTTCCACACAACGTAAATTGACCGTCCAGACTGGCAAACGCCTGTAAACCAAAGCGTAATAAAATCCCGGCAAGCATTGCCGCCGCAAGCGAGTGCGGAATAATGC  >  NC_000913/1503569‑1503752

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