Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 11.7% K641Q (AAG→CAG)  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,6200TG11.7% 176.2 / 14.1 68K641Q (AAG→CAG) abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base T (34/26);  new base G (2/6);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

ATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTG  >  NC_000913/2761479‑2761758
                                                                                                                                             |                                                                                                                                          
aTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCttt                                                                                                                                           <  1:433085/143‑1 (MQ=255)
   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)
                                                                                                                                             |                                                                                                                                          
ATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTG  >  NC_000913/2761479‑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: