Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,520,878 A→G V71A (GTG→GCG)  b1451 ← putative outer membrane receptor for iron transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,8780AG99.4% 1096.4 / ‑3.2 325V71A (GTG→GCG) b1451putative outer membrane receptor for iron transport
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (169/154);  total (171/154)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 342 total aligned reads displayed.

ATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATCC  >  NC_000913/1520788‑1520967
                                                                                          |                                                                                         
 taccgTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGgcgcgc                                                                                        >  1:3009697/1‑93 (MQ=255)
   ccgTAAGTGGAGCGGGATCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCg                                                                                      <  1:2975232/93‑1 (MQ=255)
   ccgTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCg                                                                                      >  2:7885765/1‑93 (MQ=255)
      tAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTc                                                                                   <  1:1514888/93‑1 (MQ=255)
       aaGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCa                                                                                  <  2:3960081/93‑1 (MQ=255)
          tGGAGCGGGTGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGACGCTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTg                                                                               >  2:2667176/1‑93 (MQ=255)
           ggAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGa                                                                              <  2:6745385/93‑1 (MQ=255)
            gAGCGGGAGCCAAATCCGCGAATCGACAGCTGTTACTCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAACCCAGGCGCGCCGGTCCGTGAt                                                                             <  2:7803667/93‑1 (MQ=255)
            gAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGAt                                                                             <  1:770098/93‑1 (MQ=255)
              gCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTc                                                                           <  2:2187300/93‑1 (MQ=255)
                gggAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCgg                                                                         <  2:6238688/93‑1 (MQ=255)
                 ggAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGa                                                                        >  1:1658071/1‑93 (MQ=255)
                  gAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGAc                                                                       >  2:4393537/1‑93 (MQ=255)
                  gAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGAc                                                                       >  1:5394427/1‑93 (MQ=255)
                    gCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACaa                                                                     <  1:8533476/93‑1 (MQ=255)
                       aaaTCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGtt                                                                  >  1:3561312/1‑93 (MQ=255)
                       aaaTCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGtt                                                                  >  2:6544209/1‑93 (MQ=255)
                        aaTCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTa                                                                 >  2:2070148/1‑93 (MQ=255)
                        aaTCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTa                                                                 >  2:6122593/1‑93 (MQ=255)
                         aTCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTaa                                                                >  2:69703/1‑93 (MQ=255)
                          tCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAAt                                                               <  2:4654561/93‑1 (MQ=255)
                          tCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAAt                                                               <  1:2447173/93‑1 (MQ=255)
                          cccGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGCTCAGTGATTCGGACAAGTTAAt                                                               <  2:5526326/92‑1 (MQ=255)
                           ccGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATg                                                              <  1:3883349/93‑1 (MQ=255)
                           ccGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATg                                                              >  1:5542845/1‑93 (MQ=255)
                            cgcgAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAATTTAATgc                                                             >  1:1077612/1‑93 (MQ=255)
                             gcgAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATgcg                                                            >  2:3890161/1‑93 (MQ=255)
                             gcgAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATgcg                                                            >  1:8286217/1‑93 (MQ=255)
                             gcgAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATgcg                                                            >  1:607591/1‑93 (MQ=255)
                              cgAATCGACAGCTGTGAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATgcgc                                                           <  1:7830676/93‑1 (MQ=255)
                              cgAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATgcgc                                                           <  1:4229556/93‑1 (MQ=255)
                                 aTCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGgt                                                        <  2:1968398/93‑1 (MQ=255)
                                  tCGACAGCTGTAAATCTTGAGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGgtg                                                       >  1:2785247/1‑93 (MQ=255)
                                   cGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGgtgt                                                      <  2:4106329/93‑1 (MQ=255)
                                   cGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGgtgt                                                      <  1:4742955/93‑1 (MQ=255)
                                   cGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGgtgt                                                      <  2:8009184/93‑1 (MQ=255)
                                   cGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGgtgt                                                      >  1:7305961/1‑93 (MQ=255)
                                    gACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGtt                                                     <  2:7270248/93‑1 (MQ=255)
                                      cAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGc                                                   >  1:592028/1‑93 (MQ=255)
                                      cAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGc                                                   >  2:8307658/1‑93 (MQ=255)
                                       aGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGcc                                                  >  1:513453/1‑93 (MQ=255)
                                       aGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGCGTTGcc                                                  <  1:8736821/93‑1 (MQ=255)
                                        gCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCa                                                 <  2:80022/93‑1 (MQ=255)
                                         cTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAg                                                <  2:2455537/93‑1 (MQ=255)
                                           gTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGgc                                              >  2:3817820/1‑93 (MQ=255)
                                            tAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGgcg                                             >  2:5680703/1‑93 (MQ=255)
                                             aaaTCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGgcgc                                            >  1:6433134/1‑93 (MQ=255)
                                             aaaTCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGgcgc                                            >  2:1321515/1‑93 (MQ=255)
                                             aaaTCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGgcgc                                            <  1:2019152/93‑1 (MQ=255)
                                               aTCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGgcgcgt                                          >  2:3691606/1‑91 (MQ=255)
                                                tCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGtc                                         <  1:1002578/93‑1 (MQ=255)
                                                tCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGtc                                         >  2:8550540/1‑93 (MQ=255)
                                                tCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGtc                                         >  2:1487416/1‑93 (MQ=255)
                                                  ttGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGtctc                                       <  1:1838419/93‑1 (MQ=255)
                                                  ttGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGtctc                                       >  1:2258428/1‑93 (MQ=255)
                                                  ttGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGtctc                                       <  1:6337537/93‑1 (MQ=255)
                                                   tGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGGTTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTctcc                                      >  1:3318695/1‑93 (MQ=255)
                                                   tGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTctcc                                      <  1:5299164/93‑1 (MQ=255)
                                                   tGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTctcc                                      <  2:748227/93‑1 (MQ=255)
                                                      gcATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCg                                   >  2:7353201/1‑93 (MQ=255)
                                                      gcATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCg                                   >  1:411595/1‑93 (MQ=255)
                                                         tAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCa                                >  1:4468552/1‑93 (MQ=255)
                                                         tAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCa                                <  2:5240072/93‑1 (MQ=255)
                                                           gTTCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATc                              >  2:2294017/1‑93 (MQ=255)
                                                            ttCTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATcc                             <  2:110355/93‑1 (MQ=255)
                                                              cTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATccac                           >  2:1734027/1‑93 (MQ=255)
                                                              cTGCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATccac                           >  2:2534443/1‑93 (MQ=255)
                                                                gCCGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATccacca                         >  1:568794/1‑93 (MQ=255)
                                                                 ccGGTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCcaccac                        <  2:711885/93‑1 (MQ=255)
                                                                   ggTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGc                      >  1:6901007/1‑93 (MQ=255)
                                                                    gTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCt                     >  1:1006064/1‑93 (MQ=255)
                                                                    gTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCt                     <  1:7261183/93‑1 (MQ=255)
                                                                    gTTTTGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCt                     <  1:6472801/93‑1 (MQ=255)
                                                                     ttttGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCtt                    <  2:3140409/93‑1 (MQ=255)
                                                                     ttttGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCtt                    <  2:1766926/93‑1 (MQ=255)
                                                                     ttttGTACCTGCAAACCAGGCGCGCCGGTCAGAGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCta                    >  1:4031091/1‑92 (MQ=255)
                                                                      tttGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTa                   <  2:4498506/93‑1 (MQ=255)
                                                                       ttGTGCCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTAc                  <  1:4570285/93‑1 (MQ=255)
                                                                       ttGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTAc                  >  1:7962707/1‑93 (MQ=255)
                                                                       ttGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTAc                  >  1:4796634/1‑93 (MQ=255)
                                                                        tGTCCCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTAct                 <  1:8082684/93‑1 (MQ=255)
                                                                        tGTACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTAct                 <  2:5001940/93‑1 (MQ=255)
                                                                          tACCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTActgc               <  2:5788309/93‑1 (MQ=255)
                                                                           aCCTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTActgct              >  1:6965560/1‑93 (MQ=255)
                                                                            ccTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTActgctg             >  1:5776335/1‑93 (MQ=255)
                                                                             cTGCAAACCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTActgctgc            >  2:8757624/1‑93 (MQ=255)
                                                                                cAAACCAGGCGCCCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTgg         <  2:8649137/93‑1 (MQ=255)
                                                                                 aaaCCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAAGGAGCGGAGTTGCCAGGCGCGTCTCCTCGCATACCAACACGCTTACTGCTGCTggg        >  2:5996389/1‑93 (MQ=255)
                                                                                  aaCCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTgggg       >  1:853964/1‑93 (MQ=255)
                                                                                  aaCCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTgggg       <  1:363632/93‑1 (MQ=255)
                                                                                  aaCCAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTgggg       >  1:1436482/1‑93 (MQ=255)
                                                                                  aaCCAGGCGCGCCGGTCAGTGATTCGGACAAGCTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTggag       >  2:4892335/1‑91 (MQ=255)
                                                                                  aaCCAGGCGCGCCGGTCAGTGATTAGGACAAGTTAATGCGCGGTGTTGCCAGGTGCGTCTCCTCGCCATATCCCACGCTTACTGCTGCTgggg       >  2:4311646/1‑93 (MQ=255)
                                                                                    ccAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTa     <  1:658273/93‑1 (MQ=255)
                                                                                     cAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTAt    <  2:5438543/93‑1 (MQ=255)
                                                                                     cAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTAt    >  2:7738402/1‑93 (MQ=255)
                                                                                     cAGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTAt    >  1:3201141/1‑93 (MQ=255)
                                                                                      aGGCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATc   >  2:1251698/1‑93 (MQ=255)
                                                                                       ggCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATcc  >  2:5336175/1‑93 (MQ=255)
                                                                                       ggCGCGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCGTCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATcc  <  2:3586765/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGACAAGTTAATGCGCGGTGTTGCCAGGCGCATCTCCTCGCCATCCACCACGCTTACTGCTGCTGGGGTATCC  >  NC_000913/1520788‑1520967

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