Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,623 C→A 10.5% V640L (GTA→TTA)  abpB ← CP4‑57 prophage; putative helicase YfjK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6230CA10.5% 171.5 / 12.2 68V640L (GTA→TTA) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base A (2/5);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

CCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTG  >  NC_000913/2761482‑2761758
                                                                                                                                             |                                                                                                                                       
ccAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACa                                                                                                                                        <  2:490168/143‑1 (MQ=255)
  aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGc                                                                                                                                      >  1:300145/1‑143 (MQ=255)
     caGTTTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACAACTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                   <  2:466985/142‑1 (MQ=255)
     aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtg                                                                                                                                   <  2:357601/143‑1 (MQ=255)
     aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtg                                                                                                                                   <  1:312171/143‑1 (MQ=255)
       gTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggt                                                                                                                                 <  1:502472/143‑1 (MQ=255)
          tgccAGTCGGTTGAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtata                                                                                                                              <  2:128819/141‑4 (MQ=255)
            ccAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggtatacg                                                                                                                            <  1:606887/143‑6 (MQ=255)
             ccggcGGTGTAACGTCTGAATGGGAATGCCAGTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                           <  2:817048/139‑1 (MQ=255)
             cAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATAGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                           <  1:897221/143‑1 (MQ=255)
             cAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGAACCAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                           <  1:768652/143‑1 (MQ=255)
                  ggTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTCACTACa                                                                                                                      >  1:715707/1‑143 (MQ=255)
                  ggTGTAACGTCTGAATGGGAATGACAATCTCTACTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTAGAACGCCAATTTGACATGTGAGTAATCGACTTTGCTGCGTATACAACGATTTACAGCGTGGCGTTTAACTaaa                                                                                                                      >  2:47891/1‑141 (MQ=255)
                    tgtAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                    <  1:605208/143‑1 (MQ=255)
                       aaCGTCTGAATGGGGATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACGTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcctc                                                                                                                 <  2:502229/143‑1 (MQ=255)
                            cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCa                                                                                                            <  2:49675/143‑1 (MQ=255)
                                 tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGGTGGCTCTAGCCAAACCTCAATTTGACAGTTGGGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTg                                                                                                       <  2:1023749/143‑1 (MQ=255)
                                 tGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTg                                                                                                       >  2:821543/1‑143 (MQ=255)
                                  gggAATGCCAATCTCTTCTAAGGCATTGACACCTGCTGTGAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                      <  1:816603/143‑1 (MQ=255)
                                  gggAAGGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                      <  1:637452/143‑1 (MQ=255)
                                         ccAAGCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCtt                                                                                               <  1:61038/143‑1 (MQ=255)
                                          cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTta                                                                                              <  1:351573/143‑1 (MQ=255)
                                            aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtagg                                                                                            <  1:597678/143‑1 (MQ=255)
                                             tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCCCCTtcggt                                                                                           >  1:956486/1‑143 (MQ=255)
                                                cttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                        >  2:672237/1‑143 (MQ=255)
                                                   ctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTACCTATACCACGCTTTTCAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCCTAGGTAGGGta                                                                                     >  2:369995/1‑143 (MQ=255)
                                                     aaGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACTCCTATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGGTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtata                                                                                   <  2:507116/143‑1 (MQ=255)
                                                     aaGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGAGTTTGCTGCGTATACCACGCTTTACAGCGCGGTGTTTAACTACATCGTCCATCAAGCTGAGCACCTTagggagggagtg                                                                                   >  2:200337/1‑139 (MQ=255)
                                                      aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTCGGGTACAg                                                                                  >  2:899813/1‑143 (MQ=255)
                                                      aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAACAGAGAACATGCGGTAAGGTATAg                                                                                  >  2:587151/1‑143 (MQ=255)
                                                                   tggtggTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTgg                                                                     >  1:702511/1‑143 (MQ=255)
                                                                       ggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAGAGCGTGGTGTTTAACTACATCCTCCCTCAAGCTGAGCAACCTAGGCACAGTACAGCGGAAGGTCTCGGTTac                                                                 >  2:460474/1‑142 (MQ=255)
                                                                         tAGGTGGAAACGTTCGAACGCCGATCGGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                               <  1:836257/143‑1 (MQ=255)
                                                                         tAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                               <  1:461804/143‑1 (MQ=255)
                                                                           gATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTGGCGGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCGGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGc                                                             >  2:799359/1‑143 (MQ=255)
                                                                            aTGGAAACTTTCGAACGCCAAGTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCt                                                            <  1:56250/143‑1 (MQ=255)
                                                                                 aaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTGTGAGCTTCaaa                                                       >  1:1003094/1‑143 (MQ=255)
                                                                                    tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCa                                                    >  1:716848/1‑143 (MQ=255)
                                                                                         aaCGCCAATTTTGCGTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGTTGTTTAACTACATCCTCCATCAAGCTTAGCACCTTAGGTAGGGTATAGCCGAAGTTATTGGTTATGAGCTTCAAATCAttgtt                                               <  2:894354/143‑1 (MQ=255)
                                                                                         aaCGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttgtt                                               <  1:526084/143‑1 (MQ=255)
                                                                                         aaCGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCATAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttgtt                                               >  1:250857/1‑143 (MQ=255)
                                                                                             ccAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAgg                                           <  1:865875/143‑1 (MQ=255)
                                                                                              cAATTTTCCATGTGTGTAATCGACTTTGCCGCTTATACATCGTGTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                          <  2:906911/143‑1 (MQ=255)
                                                                                              cAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAggg                                          >  1:519057/1‑143 (MQ=255)
                                                                                                aTTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCa                                        >  2:43070/1‑143 (MQ=255)
                                                                                                  ttGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCTTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAAt                                      <  2:185606/143‑1 (MQ=255)
                                                                                                       aTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGAGTTTAACTACATCCTCCATCACTCTGAGCCCCTTAGGTAGGGTATAGCCGAAGGTATTGGATATGAGCCGCAAATCATTTTTGCGGGCCATGGaaa                                 >  2:59947/1‑143 (MQ=255)
                                                                                                          tgtgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATCGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAggag                              >  2:546154/1‑143 (MQ=255)
                                                                                                            tgtAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                            >  2:858882/1‑143 (MQ=255)
                                                                                                              tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACg                          >  1:224265/1‑143 (MQ=255)
                                                                                                                aTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGATCTCCGCTTCTGATCTGTGTTGAGGGCAATGGAAAGAGTACGtt                        >  2:216902/1‑143 (MQ=255)
                                                                                                                 tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAGTGTTGAGGGCAATGGAAAGAGTACGttt                       >  2:248515/1‑143 (MQ=255)
                                                                                                                  cGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTACATGAAGCTGGGCACCTCAGGGCAGTTAGAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtc                      >  1:1048248/1‑143 (MQ=255)
                                                                                                                    aCTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCGTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTGAAATCATTGTTGAGGGCAATGGAAAGAGAACGGTtgtc                    >  2:574419/1‑143 (MQ=255)
                                                                                                                    aCTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtctc                    <  1:563610/143‑1 (MQ=255)
                                                                                                                    aCTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTtctc                    <  1:561539/143‑1 (MQ=255)
                                                                                                                           tgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCa             >  2:302160/1‑143 (MQ=255)
                                                                                                                             cGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCGATAGAAAGAGTACGTTTCTCCCTTCCAgg           >  1:338289/1‑142 (MQ=255)
                                                                                                                              aaacaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATCGTTGAGGGCAATGGAAAGAGTACGCTTCTACACTTCAGaa          >  2:662892/5‑143 (MQ=255)
                                                                                                                                atatctcGCTTTACAGCGTGTTGTTTAACTACATCCTCCATCAAGCTGAGCTCCTGAGGGAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAAt        <  2:511992/137‑1 (MQ=255)
                                                                                                                                ataCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCAGCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAAt        <  2:300145/143‑1 (MQ=255)
                                                                                                                                  aCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATcc      >  2:220578/1‑143 (MQ=255)
                                                                                                                                   ccACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa     >  1:50179/1‑143 (MQ=255)
                                                                                                                                   ccACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCa     >  1:334277/1‑143 (MQ=255)
                                                                                                                                    cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAACTCCAc    >  1:281942/1‑143 (MQ=255)
                                                                                                                                    cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCAc    >  2:33127/1‑143 (MQ=255)
                                                                                                                                     aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTATGGTTCTCCCCTTCAGAAAACCACt   >  2:20585/1‑143 (MQ=255)
                                                                                                                                      cGCTTTACAGCGTGGTGTTTAACTATATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTg  >  1:796974/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                       
CCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTG  >  NC_000913/2761482‑2761758

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