Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,761,620 T→G 9.7% K641Q (AAG→CAG)  yfjK ← radiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,761,6200TG 9.7% 472.1 / 40.1 144K641Q (AAG→CAG) yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base T (62/68);  new base G (7/7);  total (69/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 147 total aligned reads displayed.

TCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATA  >  NC_000913/2761472‑2761767
                                                                                                                                                    |                                                                                                                                                   
   gAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTAc                                                                                                                                                  >  4:405467/1‑149 (MQ=255)
    aaaaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACa                                                                                                                                                 >  3:69050/1‑149 (MQ=255)
    aaaaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACa                                                                                                                                                 <  1:145650/149‑1 (MQ=255)
      aaTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGc                                                                                                                                               <  2:578803/149‑1 (MQ=255)
       aTTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAg                                                                                                                                                >  4:554682/1‑147 (MQ=255)
        ttCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgt                                                                                                                                             <  3:504457/149‑1 (MQ=255)
           cAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtggt                                                                                                                                          <  3:74976/149‑1 (MQ=255)
            aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                         >  2:263935/1‑149 (MQ=255)
            aGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggtg                                                                                                                                         >  4:651145/1‑149 (MQ=255)
             gCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGgtgt                                                                                                                                        >  1:523527/1‑149 (MQ=255)
              cAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGtt                                                                                                                                       <  4:69050/149‑1 (MQ=255)
               aaGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGttt                                                                                                                                      <  3:186904/149‑1 (MQ=255)
                aGTCTACCAGTCGGTGTAACGTCTGAATGGGCATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTa                                                                                                                                     >  1:390915/1‑149 (MQ=255)
                 gTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTaa                                                                                                                                    >  3:31664/1‑149 (MQ=255)
                  tCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAAc                                                                                                                                   <  1:373533/149‑1 (MQ=255)
                    tACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTa                                                                                                                                 >  4:671356/1‑149 (MQ=255)
                      ccAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACa                                                                                                                               >  3:159277/1‑149 (MQ=255)
                      ccAGTCGGTGTAACGTCTGAATGGGAATGCCAATCCCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACa                                                                                                                               <  4:450254/149‑1 (MQ=255)
                       cAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTGTAAAGCgtg                                                                                                                                            <  4:263720/135‑1 (MQ=255)
                        aGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtc                                                                                                                             <  4:330854/149‑1 (MQ=255)
                              tgtAACGTCTGAAAGGGAATGCCAATCTCTTCTAAGGAATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                            <  4:640155/128‑1 (MQ=255)
                                tAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCa                                                                                                                     >  1:59943/1‑149 (MQ=255)
                                 aaCGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCaa                                                                                                                    <  4:69838/149‑1 (MQ=255)
                                    gTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCt                                                                                                                 >  3:344078/1‑149 (MQ=255)
                                     tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTACACCACGCTGTAAAGCgtggta                                                                                                                                         >  4:126934/1‑123 (MQ=255)
                                     tCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTACACCACGCTGTAAAGCgtggta                                                                                                                                         <  3:126934/124‑2 (MQ=255)
                                      cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                               <  3:420492/149‑1 (MQ=255)
                                      cTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                               <  4:521124/149‑1 (MQ=255)
                                       tGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACAGGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTATTTAACTACATCCTCCATCAAGCGGAg                                                                                                              >  4:445585/1‑149 (MQ=255)
                                        gAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGc                                                                                                             <  1:126924/149‑1 (MQ=255)
                                         aaTGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCa                                                                                                            >  2:554876/1‑149 (MQ=255)
                                                   ccAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTATGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtag                                                                                                  >  2:88341/1‑149 (MQ=255)
                                                    cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcctc                                                                                                                          <  2:131199/124‑1 (MQ=255)
                                                    cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACAtcctc                                                                                                                          >  1:131199/1‑124 (MQ=255)
                                                    cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                 >  1:431072/1‑149 (MQ=255)
                                                    cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                 >  3:263074/1‑149 (MQ=255)
                                                    cAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtagg                                                                                                 <  3:23634/149‑1 (MQ=255)
                                                      aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggt                                                                                                                                          >  4:272590/1‑106 (MQ=255)
                                                      aTCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtggt                                                                                                                                          <  3:272590/106‑1 (MQ=255)
                                                       tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTACACCACGCTGTAAAGCgtggtatac                                                                                                                                      <  2:171622/109‑5 (MQ=255)
                                                       tctcTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTACACCACGCTGTAAAGCgtggtatac                                                                                                                                      >  1:171622/1‑105 (MQ=255)
                                                         tcttctAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTTTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGtata                                                                                            >  4:412573/1‑149 (MQ=255)
                                                              tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                       >  2:335222/1‑149 (MQ=255)
                                                              tAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGa                                                                                       >  3:540897/1‑149 (MQ=255)
                                                               aaGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGaa                                                                                      <  4:31664/149‑1 (MQ=255)
                                                                aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                     <  4:540897/149‑1 (MQ=255)
                                                                aGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                     <  4:345762/149‑1 (MQ=255)
                                                                      tGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTg                                                                               <  2:431072/149‑1 (MQ=255)
                                                                        acacCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGcc                                                                                         <  2:448476/137‑1 (MQ=255)
                                                                             tggtggTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGa                                                                        <  3:671356/149‑1 (MQ=255)
                                                                               gtggtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGc                                                                      <  2:390915/149‑1 (MQ=255)
                                                                                  gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                          <  3:269922/142‑1 (MQ=255)
                                                                                  gtAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                          >  4:269922/1‑142 (MQ=255)
                                                                                       tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAg                                                                                           <  3:640540/120‑1 (MQ=255)
                                                                                       tGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaaa                                                                <  1:61316/147‑1 (MQ=255)
                                                                                        ggAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                            <  2:209574/70‑1 (MQ=255)
                                                                                          aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                           <  4:328937/149‑1 (MQ=255)
                                                                                          aaaCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                           <  3:651145/149‑1 (MQ=255)
                                                                                            aCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                               <  4:102878/95‑1 (MQ=255)
                                                                                            aCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGa                                                                                                               >  3:102878/1‑95 (MQ=255)
                                                                                            aCTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAttgt                                                         >  3:681887/1‑149 (MQ=255)
                                                                                              tttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCAcc                                                                                                          >  1:673938/1‑98 (MQ=255)
                                                                                               ttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                      <  2:488941/149‑1 (MQ=255)
                                                                                               ttCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGa                                                      >  4:143257/1‑149 (MQ=255)
                                                                                                       ccAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAg                                                                                     >  1:449374/1‑110 (MQ=255)
                                                                                                         aaTTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCgtg                                                                                                                                            <  1:194509/53‑1 (MQ=255)
                                                                                                           tttGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                          >  3:377146/1‑117 (MQ=255)
                                                                                                              gACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAagag                                       <  2:193520/149‑1 (MQ=255)
                                                                                                                cATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTATTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                     <  3:570206/149‑1 (MQ=255)
                                                                                                                 aTGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTAc                                    >  4:350427/1‑149 (MQ=255)
                                                                                                                        tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCTAAGGTATTGGTTATGAGCTTCaaa                                                                <  2:65560/114‑1 (MQ=255)
                                                                                                                        tAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCaaa                                                                >  1:65560/1‑114 (MQ=255)
                                                                                                                         aaTCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTcc                            <  1:263935/149‑1 (MQ=255)
                                                                                                                           tCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTcccc                          <  3:143257/149‑1 (MQ=255)
                                                                                                                             gACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCtt                        >  1:305713/1‑149 (MQ=255)
                                                                                                                                    ctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGCGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAAt                 <  3:445585/149‑1 (MQ=255)
                                                                                                                                    ctgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAAGGTATAGCCGAAGGTATTTGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAGAGAGTACGTTTCTCCCCTTCCGAAAt                 <  4:590123/149‑1 (MQ=255)
                                                                                                                                     tgcGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATc                >  1:296259/1‑149 (MQ=255)
                                                                                                                                         tataCCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACt            <  4:159277/149‑1 (MQ=255)
                                                                                                                                          acaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTAt                                                                          >  2:631767/3‑86 (MQ=255)
                                                                                                                                          acaCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTg           >  1:342361/3‑149 (MQ=255)
                                                                                                                                            aCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                     >  1:503929/1‑121 (MQ=255)
                                                                                                                                            aCCACGCTGTAAAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTa                                     <  2:503929/121‑1 (MQ=255)
                                                                                                                                             ccACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTAc                                    <  4:314456/121‑1 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCAtt                                                           <  3:318317/97‑1 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAaga                                        <  3:576666/116‑1 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt       >  1:223642/1‑149 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt       >  1:200616/1‑149 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt       >  1:120835/1‑149 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt       >  4:220829/1‑149 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt       <  3:412573/149‑1 (MQ=255)
                                                                                                                                              cACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGt       >  2:373789/1‑149 (MQ=255)
                                                                                                                                               aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTtaggtag                                                                                                  <  3:279228/57‑1 (MQ=255)
                                                                                                                                               aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg      <  2:305713/149‑1 (MQ=255)
                                                                                                                                               aCGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCATAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAGTCATTGTAGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTg      >  4:119896/1‑149 (MQ=255)
                                                                                                                                                cGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGc     <  3:416811/149‑1 (MQ=255)
                                                                                                                                                 gCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCa    >  4:352485/1‑149 (MQ=255)
                                                                                                                                                 gCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCa    <  2:223642/149‑1 (MQ=255)
                                                                                                                                                   tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATa  >  4:475233/1‑149 (MQ=255)
                                                                                                                                                   tttACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATa  >  2:147926/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGGAAAATTCCAGCAAGTCTACCAGTCGGTGTAACGTCTGAATGGGAATGCCAATCTCTTCTAAGGCATTGACACCTGGTGGTAGATGGAAACTTTCGAACGCCAATTTGACATGTGTGTAATCGACTTTGCTGCGTATACCACGCTTTACAGCGTGGTGTTTAACTACATCCTCCATCAAGCTGAGCACCTTAGGTAGGGTATAGCCGAAGGTATTGGTTATGAGCTTCAAATCATTGTTGAGGGCAATGGAAAGAGTACGTTTCTCCCCTTCAGAAATCCACTGTCGTGCATA  >  NC_000913/2761472‑2761767

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