Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,541,456 A→G 23.7% I75T (ATT→ACT)  cysW ← sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,541,4560AG23.7% 245.2 / 83.3 143I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (59/50);  new base G (17/17);  total (76/67)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 148 total aligned reads displayed.

ATGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGA  >  NC_000913/2541357‑2541555
                                                                                                   |                                                                                                   
aTGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCaa                                                                                                     <  1:2948343/100‑1 (MQ=255)
 tGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAAt                                                                                                    >  2:794804/1‑100 (MQ=255)
  gTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAGGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATc                                                                                                   >  2:2038055/1‑100 (MQ=255)
    ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAg                                                                                                 <  1:3369816/100‑1 (MQ=255)
       gTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAgcgc                                                                                              >  1:3296667/1‑100 (MQ=255)
       gTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAgcgc                                                                                              >  2:1497278/1‑100 (MQ=255)
            gTCAGCAGTAACGGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCGAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatca                                                                                         <  2:322877/100‑1 (MQ=255)
             tCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcat                                                                                        >  2:387696/1‑100 (MQ=255)
              cagcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcatc                                                                                       >  1:1048320/1‑100 (MQ=255)
                 cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCAcc                                                                                    >  1:2818521/1‑100 (MQ=255)
                  agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCg                                                                                   <  2:3011117/100‑1 (MQ=255)
                  agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCg                                                                                   >  2:604999/1‑100 (MQ=255)
                   gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGt                                                                                  <  1:4023052/100‑1 (MQ=255)
                    tAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTc                                                                                 <  1:248628/100‑1 (MQ=255)
                      aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                                               <  1:3872802/100‑1 (MQ=255)
                      aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                                               <  2:617311/100‑1 (MQ=255)
                      aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                                               <  1:1497290/100‑1 (MQ=255)
                       cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGc                                                                              >  1:4099583/1‑100 (MQ=255)
                        tGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGcc                                                                             <  1:2965526/100‑1 (MQ=255)
                         gtcGTCCAGGGAAGTCAATCCGCGTTACCAGCCAGCGCAGCAGAATGCCGAACTCCAGGTTTACCGGTACGGCAGTCAGCGCGATCAGCACCGTCAGc                                                                              <  2:34319/96‑1 (MQ=255)
                         ggCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGcc                                                                             <  1:2969687/99‑1 (MQ=255)
                         ggCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                                            >  2:813302/1‑100 (MQ=255)
                            gTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGAt                                                                         <  1:1647620/100‑1 (MQ=255)
                            gTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGAt                                                                         >  1:2295053/1‑100 (MQ=255)
                             tCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATg                                                                        <  2:29123/100‑1 (MQ=255)
                              ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATgg                                                                       <  1:3948891/100‑1 (MQ=255)
                              ccAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATgg                                                                       <  1:158144/100‑1 (MQ=255)
                               cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                      >  2:1895485/1‑100 (MQ=255)
                               cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                      >  2:1136178/1‑100 (MQ=255)
                               cAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGCTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGc                                                                      >  2:1727653/1‑100 (MQ=255)
                                 gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCCATCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                                                    >  2:3430770/1‑100 (MQ=255)
                                    aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                                 <  2:1918157/100‑1 (MQ=255)
                                      gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                                               >  1:3957185/1‑100 (MQ=255)
                                      gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                                               >  1:3635161/1‑100 (MQ=255)
                                        tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                             <  1:3609023/100‑1 (MQ=255)
                                        tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                             >  2:2798263/1‑100 (MQ=255)
                                         aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                                            <  1:4055984/100‑1 (MQ=255)
                                         aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                                            >  2:2977031/1‑100 (MQ=255)
                                           aGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAGCACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTc                                                          >  2:3636373/1‑100 (MQ=255)
                                                gTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAt                                                     >  1:1613863/1‑100 (MQ=255)
                                                 tCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATc                                                    >  1:2880413/1‑100 (MQ=255)
                                                   aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                                  >  1:3013600/1‑100 (MQ=255)
                                                   aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                                  >  1:501763/1‑100 (MQ=255)
                                                   aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                                  <  2:4049459/100‑1 (MQ=255)
                                                   aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCgg                                                  <  1:3421898/100‑1 (MQ=255)
                                                   aCCAGCCAGGCCAGCAGAATGCAGAACACCAGGTTTACCGGTAAGGCAGTCAGCGCAATCATCACCGTCCGCCCGCTGGCGTGCAGCATGTCCGGATccg                                                  >  2:16544/1‑98 (MQ=255)
                                                     cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                                <  1:995630/100‑1 (MQ=255)
                                                     cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                                <  2:2768358/100‑1 (MQ=255)
                                                       gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                               <  2:3908948/99‑1 (MQ=255)
                                                        ccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGa                                             >  1:2559077/1‑100 (MQ=255)
                                                         caggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAt                                            >  1:3020602/1‑100 (MQ=255)
                                                          aggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAtt                                           >  1:2186175/1‑100 (MQ=255)
                                                          aggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAtt                                           >  2:532085/1‑100 (MQ=255)
                                                              cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                       <  2:1180442/100‑1 (MQ=255)
                                                              cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                       >  1:1054449/1‑100 (MQ=255)
                                                               agcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                      <  2:2374215/100‑1 (MQ=255)
                                                                gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa                                     >  1:3270995/1‑100 (MQ=255)
                                                                gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa                                     <  1:1748555/100‑1 (MQ=255)
                                                                gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa                                     >  1:1266863/1‑100 (MQ=255)
                                                                gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGACAGATTCTGTan                                     <  1:2716643/100‑2 (MQ=255)
                                                                  agAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaaaa                                   >  1:180105/1‑100 (MQ=255)
                                                                  agAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaaaa                                   <  2:1198943/100‑1 (MQ=255)
                                                                   gAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAc                                  >  2:113909/1‑100 (MQ=255)
                                                                   gAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAc                                  >  1:1130598/1‑100 (MQ=255)
                                                                     aTGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCg                                <  2:3013581/100‑1 (MQ=255)
                                                                     aTGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCg                                <  1:1856869/100‑1 (MQ=255)
                                                                     aTGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCg                                <  2:1989346/100‑1 (MQ=255)
                                                                      tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCgg                               <  1:347390/100‑1 (MQ=255)
                                                                       gCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                              <  2:3653459/100‑1 (MQ=255)
                                                                       gCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGc                              >  1:1149196/1‑100 (MQ=255)
                                                                        ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                             >  1:3155315/1‑100 (MQ=255)
                                                                        ccGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                             >  2:546089/1‑100 (MQ=255)
                                                                           aaCACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCa                          >  1:1054366/1‑100 (MQ=255)
                                                                            acacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAg                         <  1:2924510/100‑1 (MQ=255)
                                                                            acacCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                            >  2:170709/1‑97 (MQ=255)
                                                                             cacCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGc                        >  1:851900/1‑100 (MQ=255)
                                                                              atctggTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                       <  2:4086771/96‑1 (MQ=255)
                                                                              acCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                       >  2:2587319/1‑100 (MQ=255)
                                                                               ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                      >  1:3887415/1‑100 (MQ=255)
                                                                               ccAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                       >  2:2427388/1‑99 (MQ=255)
                                                                               ccAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                      <  1:2038068/100‑1 (MQ=255)
                                                                                cAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcccc                     <  1:605005/100‑1 (MQ=255)
                                                                                 aGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAGTCTGTAGAACCGGCATCAGc                        >  2:2859035/1‑96 (MQ=255)
                                                                                 aGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCt                    >  1:2187760/1‑100 (MQ=255)
                                                                                 aGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCt                    >  1:413292/1‑100 (MQ=255)
                                                                                       aCCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGa              <  2:334405/100‑1 (MQ=255)
                                                                                        ccGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGa              >  2:1529226/1‑99 (MQ=255)
                                                                                         cGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAAt            <  2:3272060/100‑1 (MQ=255)
                                                                                          ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg           <  1:25006/100‑1 (MQ=255)
                                                                                          ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg           <  1:155409/100‑1 (MQ=255)
                                                                                          ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg           <  1:874064/100‑1 (MQ=255)
                                                                                          ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg           >  2:299390/1‑100 (MQ=255)
                                                                                          ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGAGGGCGTGCAGCATGTCCGGAGCGGCCAGATTCTGGAAAACCGGCATCAGCCCCTTGCTGAATg           >  2:2903911/1‑100 (MQ=255)
                                                                                           gTACGGCAATCAGCGCGGTCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg           <  2:2238069/99‑1 (MQ=255)
                                                                                             aCGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc         <  2:1779046/99‑1 (MQ=255)
                                                                                              cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCCTGCTGAATGCCTg       <  2:2556836/100‑1 (MQ=255)
                                                                                              cGGCAATAAGCACGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCGGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg       >  2:2171553/1‑100 (MQ=255)
                                                                                              cGGCAACCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTg       <  1:3582807/100‑1 (MQ=255)
                                                                                                gCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCa     >  2:3309529/1‑100 (MQ=255)
                                                                                                   aTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGa  <  1:3228813/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATGTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACGA  >  NC_000913/2541357‑2541555

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