Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 451,404 A→G F456L (TTC→CTC)  ampG ← AmpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913451,4040AG99.3% 922.2 / ‑3.3 293F456L (TTC→CTC) ampGAmpG
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (124/167);  total (125/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 320 total aligned reads displayed.

TAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAGG  >  NC_000913/451315‑451493
                                                                                          |                                                                                         
gaGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGGAAAGAGAGGt                                                                                         <  2:7556010/92‑1 (MQ=255)
    gcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTc                                                                                    >  2:8158136/1‑93 (MQ=255)
    gcCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGAGCGTc                                                                                    >  2:7781811/1‑93 (MQ=255)
      cAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCaa                                                                                  >  1:5035357/1‑93 (MQ=255)
       aGATAATCCAGCAAACCACCAAGCACGACGCCAGAACGGGCTACTAAACCCCCGGCTTCCAGCAGACCAGGCCGGAAAGAGAGGCGCGCCaaa                                                                                 >  2:2363572/1‑93 (MQ=255)
       aGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCaaa                                                                                 >  1:711237/1‑93 (MQ=255)
       aGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCaaa                                                                                 >  1:6567349/1‑93 (MQ=255)
         aTAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATc                                                                               <  2:3619261/93‑1 (MQ=255)
          tAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATcc                                                                              <  2:2253192/93‑1 (MQ=255)
          tAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATcc                                                                              <  2:8243272/93‑1 (MQ=255)
            aTCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAACTCCAg                                                                            >  1:9592186/1‑93 (MQ=255)
             tCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgc                                                                           >  1:1189892/1‑93 (MQ=255)
               cAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAgcgc                                                                         >  1:2787094/1‑93 (MQ=255)
                  cAAACCACCAAGCACGACGCCAGCAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTc                                                                      <  2:7208468/93‑1 (MQ=255)
                  cAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTc                                                                      <  2:8180580/93‑1 (MQ=255)
                    aaCCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCa                                                                    >  2:2781163/1‑93 (MQ=255)
                      ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                  >  1:7972810/1‑93 (MQ=255)
                      ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                  >  1:5341097/1‑93 (MQ=255)
                      ccaccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATc                                                                  <  1:4544528/93‑1 (MQ=255)
                       caccaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCg                                                                 >  1:6117548/1‑93 (MQ=255)
                        accaAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGt                                                                >  1:1009960/1‑93 (MQ=255)
                            aGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAAt                                                            >  2:7954166/1‑93 (MQ=255)
                            aGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAAt                                                            >  2:6136882/1‑93 (MQ=255)
                             gCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATa                                                           >  2:5968038/1‑93 (MQ=255)
                             gCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATa                                                           <  2:9660135/93‑1 (MQ=255)
                               acgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGc                                                         <  1:2888626/93‑1 (MQ=255)
                                cgacgCCAGAAAGGGCGACTAAACCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                                        <  1:7661670/93‑1 (MQ=255)
                                cgacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCa                                                        <  2:6452238/93‑1 (MQ=255)
                                 gacgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAg                                                       <  1:6789747/93‑1 (MQ=255)
                                  acgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGt                                                      >  2:3727972/1‑93 (MQ=255)
                                  acgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGt                                                      >  2:3977157/1‑93 (MQ=255)
                                  acgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGGGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAga                                                      >  1:486449/1‑92 (MQ=255)
                                   cgCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGAGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTa                                                     >  2:6066531/1‑93 (MQ=255)
                                     ccAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAc                                                   <  2:6339854/93‑1 (MQ=255)
                                     ccAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAc                                                   >  1:4251965/1‑93 (MQ=255)
                                      cAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAAcc                                                  <  1:9153619/93‑1 (MQ=255)
                                       aGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACca                                                 >  2:6890597/1‑93 (MQ=255)
                                       aGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACca                                                 <  1:4720788/93‑1 (MQ=255)
                                       aGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAAACCAAAACGACCATCGTCAATAGCAGTAACca                                                 <  1:8940090/93‑1 (MQ=255)
                                          aaGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCacac                                              <  2:4666799/93‑1 (MQ=255)
                                           aGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACg                                             >  1:4614550/1‑93 (MQ=255)
                                           aGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACg                                             <  1:2098592/93‑1 (MQ=255)
                                                 aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                       >  2:5423767/1‑93 (MQ=255)
                                                 aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                       >  1:8546070/1‑93 (MQ=255)
                                                 aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                       <  2:1764330/93‑1 (MQ=255)
                                                 aCTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAAc                                       >  1:2319755/1‑93 (MQ=255)
                                                  cTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                                      <  1:9510971/93‑1 (MQ=255)
                                                  cTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                                      <  2:5773628/93‑1 (MQ=255)
                                                  cTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACa                                      <  1:2518998/93‑1 (MQ=255)
                                                   tAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAg                                     <  2:9028214/93‑1 (MQ=255)
                                                   tAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAg                                     >  2:187900/1‑93 (MQ=255)
                                                     aaaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGc                                   >  2:2849578/1‑93 (MQ=255)
                                                      aaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCt                                  <  2:8747087/93‑1 (MQ=255)
                                                      aaCCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACACGCt                                  >  1:4415224/1‑93 (MQ=255)
                                                        ccccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGa                                <  1:7480296/93‑1 (MQ=255)
                                                        agccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCAATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGa                                <  1:5981685/91‑1 (MQ=255)
                                                          ccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcg                              >  2:4959233/1‑93 (MQ=255)
                                                          ccGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcg                              >  2:7265815/1‑93 (MQ=255)
                                                           cGACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgc                             <  1:9902865/93‑1 (MQ=255)
                                                            gACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc                            <  2:2638682/93‑1 (MQ=255)
                                                            gACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc                            <  1:5616518/93‑1 (MQ=255)
                                                            gACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc                            >  1:8733745/1‑93 (MQ=255)
                                                            gACTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcc                            >  1:8176034/1‑93 (MQ=255)
                                                             aCTTCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATTGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgccc                           >  1:930202/1‑93 (MQ=255)
                                                               ttCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccgc                         <  1:3547326/93‑1 (MQ=255)
                                                               ttCCAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGAcgcccgc                         <  1:898276/93‑1 (MQ=255)
                                                                 ccAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCccgccg                       >  2:4384082/1‑93 (MQ=255)
                                                                 ccAGCAGAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCccgccg                       >  2:6525076/1‑93 (MQ=255)
                                                                  cagcagAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCcgccgc                      >  2:56991/1‑93 (MQ=255)
                                                                  cagcagAGCAGGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCcgccgc                      >  1:7762126/1‑93 (MQ=255)
                                                                   agcagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCgccgcc                     >  2:4556789/1‑93 (MQ=255)
                                                                     gagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAg                   <  1:7303224/92‑1 (MQ=255)
                                                                     cagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGCCGCCCGCCGCCAg                   >  1:5250014/1‑93 (MQ=255)
                                                                     cagagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAg                   <  2:3178007/93‑1 (MQ=255)
                                                                      agagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgt                  <  1:6451440/93‑1 (MQ=255)
                                                                       gagcagGCAGGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAgtg                 >  1:7813448/1‑93 (MQ=255)
                                                                         gcaggcagGAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTc               <  1:1329327/93‑1 (MQ=255)
                                                                          gaggcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTcc              <  2:9590490/92‑1 (MQ=255)
                                                                            ggcaggAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcac            >  2:2992353/1‑93 (MQ=255)
                                                                             gcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  2:9786951/93‑1 (MQ=255)
                                                                             gcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  2:9349475/93‑1 (MQ=255)
                                                                             gcaggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  1:5157665/93‑1 (MQ=255)
                                                                             gcagaAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  2:84565/93‑1 (MQ=255)
                                                                             gaaggAAAGAGAGGTCCGTCAAACCCAGCGCGTCTATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCcaca           <  2:7504003/91‑1 (MQ=255)
                                                                              caggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACAt          <  1:6901703/93‑1 (MQ=255)
                                                                               aggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATg         <  1:1764339/93‑1 (MQ=255)
                                                                               aggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATg         <  2:1192986/93‑1 (MQ=255)
                                                                               aggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATg         >  2:8656604/1‑93 (MQ=255)
                                                                               aggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATg         <  2:886341/93‑1 (MQ=255)
                                                                                 ggAAAGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACCCGGCCAACAGTCTGCCGCCCCCCGCCAGGGTCCACATgg        >  1:8335825/1‑93 (MQ=255)
                                                                                    aaGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaa     <  2:3376530/93‑1 (MQ=255)
                                                                                    aaGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaa     >  1:3887130/1‑93 (MQ=255)
                                                                                    aaGAGAGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaa     >  1:4017982/1‑93 (MQ=255)
                                                                                     agagagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaag    >  2:1049540/1‑92 (MQ=255)
                                                                                     agagagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaaa    <  1:5410639/93‑1 (MQ=255)
                                                                                     agagagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaaa    <  2:9172587/93‑1 (MQ=255)
                                                                                     agagagGCGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCaaa    >  1:6621585/1‑93 (MQ=255)
                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAg   <  2:1633418/93‑1 (MQ=255)
                                                                                      gagaGGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAg   <  2:5249098/93‑1 (MQ=255)
                                                                                       agagGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAgg  <  1:1604505/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAGCGCCAGATAATCCAGCAAACCACCAAGCACGACGCCAGAAAGGGCGACTAAAACCCCGACTTCCAGCAGAGCAGGCAGGAAAGAGAAGTGCGTCAAATCCAGCGCGTCCATCGTCAATAGCAGTAACCACACGGCCAACAGGCTGACGCCCGCCGCCAGTGTCCACATGGCAAAGG  >  NC_000913/451315‑451493

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