Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,930 G→C 5.2% S500C (TCC→TGC)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,9300GC 5.2% 418.0 / 12.5 135S500C (TCC→TGC) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base G (53/73);  new base C (4/3);  total (57/78)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

Only 100 of 135 total aligned reads displayed.

TAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960785‑2961076
                                                                                                                                                 |                                                                                                                                                  
tAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTc                                                                                                                                                >  2:530382/1‑150 (MQ=255)
 aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATCCGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                               <  2:16052/150‑1 (MQ=255)
   aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                             <  1:1740789/150‑1 (MQ=255)
   aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                             <  2:352933/150‑1 (MQ=255)
   aaCCGATTGCCACTCGCCCTGAGCGCTCTCCATCTCGTAGCCCAACAACATACGCGTCAGGCCAAATCGACAGGTGTGTTGTCCGGTGTCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                             <  1:221492/150‑1 (MQ=255)
    aCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTAGTCCGCTGGCGGGGAGTTCCAAATAGATAACGATGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTg                                                                                                                                            >  2:643348/1‑150 (MQ=255)
      cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTTAc                                                                                                                                          >  2:1462267/1‑150 (MQ=255)
      cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAc                                                                                                                                          <  1:242206/150‑1 (MQ=255)
      cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAc                                                                                                                                          <  1:650321/150‑1 (MQ=255)
      cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggg                                                                                                                                          <  1:1726293/150‑3 (MQ=255)
       gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAcc                                                                                                                                         <  1:924979/150‑1 (MQ=255)
       gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgggg                                                                                                                                         <  2:472351/150‑4 (MQ=255)
       gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgggg                                                                                                                                         <  1:1439513/150‑4 (MQ=255)
          tGCCACTCACCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGGTGACCCAc                                                                                                                                      >  1:155043/1‑150 (MQ=255)
           gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                     <  2:1225969/150‑1 (MQ=255)
            ccACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTAGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCCCTg                                                                                                                                    >  2:432630/1‑150 (MQ=255)
             cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACATCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAGTCGTTGACCCACTgg                                                                                                                                   <  1:645697/150‑1 (MQ=255)
             cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGGTGACCCACTgg                                                                                                                                   >  1:1529521/1‑150 (MQ=255)
               cTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACCACATACCCGTCAGGCCAAAACGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCg                                                                                                                                 <  1:1883981/150‑1 (MQ=255)
               cTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCg                                                                                                                                 <  1:884790/150‑1 (MQ=255)
                 cGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCCCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTa                                                                                                                               >  1:1537192/1‑150 (MQ=255)
                  gCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGAGTGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTaa                                                                                                                              <  2:859893/150‑1 (MQ=255)
                   cccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                             <  1:1276951/150‑1 (MQ=255)
                    ccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCCGGATTCGTTGACCCACTGGCGTaact                                                                                                                            >  1:1581157/1‑148 (MQ=255)
                     cTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATa                                                                                                                           <  2:1983297/150‑1 (MQ=255)
                           gcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAg                                                                                                                     >  2:1745205/1‑150 (MQ=255)
                            ctctCCCTCGCGTAGCCCAACAACCTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGCATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                                    >  1:1584539/1‑150 (MQ=255)
                            ctctCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCGCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                                    <  1:413369/150‑1 (MQ=255)
                            ctctCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                                    >  2:1728580/1‑150 (MQ=255)
                                     gcgTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTAGAACTCTGTCGTCGACTTAAAAATAAAAAAAAATATCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGGAACTACCGCAGCCCTtcgtcg                                                                                                           >  2:657505/1‑150 (MQ=255)
                                       gTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGt                                                                                                         <  2:913974/150‑1 (MQ=255)
                                         aGCCCAACAACATCCGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                       <  2:1167818/150‑1 (MQ=255)
                                         aGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACCGCAGCCCTTCTTCGGTGa                                                                                                       >  2:1892276/1‑150 (MQ=255)
                                           cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                     <  1:43542/150‑1 (MQ=255)
                                           cccAACAACATACGCGTCAGGCCAAATCCCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                     <  2:1537192/150‑1 (MQ=255)
                                            ccAACACCATCCGCGTCAGGCCAACTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGt                                                                                                    <  1:274707/150‑1 (MQ=255)
                                             caacaaCATACGCGTCAGGGCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                   <  2:1549568/150‑1 (MQ=255)
                                             caacaaCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCTGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                   <  1:385423/150‑1 (MQ=255)
                                             caacaaCATACGCGTCAGGCCAAAACGCCAGGTGTGTTGTCCGGTGGCGTAGAGTTCCAGATATCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCATTCTTCGGTGATGTc                                                                                                   >  2:179021/1‑150 (MQ=255)
                                                caacaTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCAATCTCGCGAAAGTTGTCGTCATCTATGCCCCAACGAATGCCCGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCaaa                                                                                                >  2:163129/1‑150 (MQ=255)
                                                    aTACGCGTCAGGCCAAATCGCCAGGTGTGTTGACCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgc                                                                                            >  1:1098793/1‑150 (MQ=255)
                                                     tACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGCTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGGAAATAACGCCGCCCGTCGTCGGTGATGTCAAACcgcg                                                                                           >  2:285110/1‑150 (MQ=255)
                                                       cgcgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgcc                                                                                         >  2:1431333/1‑150 (MQ=255)
                                                        gcgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccg                                                                                        >  1:1974816/1‑150 (MQ=255)
                                                         cgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccgc                                                                                       >  2:1702724/1‑150 (MQ=255)
                                                         cgTCAGGCCAAATCGCCAGGTGTGTTGACCGGTGGCGGGGAGTTCCAGCTAGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccgc                                                                                       >  2:42697/1‑150 (MQ=255)
                                                          gTCAGGCCATATCGCCAGGTGTGTTGGCCGGTGGCGGGGAGTTCCAGCTCGCGAACTTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCgccgcc                                                                                      <  1:561227/150‑1 (MQ=255)
                                                            cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAGCCGCGCCGCCAg                                                                                    <  2:1114678/150‑1 (MQ=255)
                                                            cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAg                                                                                    <  1:1745205/150‑1 (MQ=255)
                                                                ccAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAcc                                                                                <  2:1841339/150‑1 (MQ=255)
                                                                  aaaTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCgg                                                                              >  2:1845251/1‑150 (MQ=255)
                                                                     tCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcac                                                                           <  1:230612/150‑1 (MQ=255)
                                                                      cGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcaca                                                                          >  1:1212600/1‑150 (MQ=255)
                                                                         cAGGTGTGTTGTGCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCTGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  2:1333/150‑1 (MQ=255)
                                                                         cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATATCGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       >  2:599241/1‑150 (MQ=255)
                                                                         cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  1:530382/150‑1 (MQ=255)
                                                                         cAGGTGTGTTGTCCGGTGGCGGGGAGGGCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  2:475479/150‑1 (MQ=255)
                                                                           ggTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAg                                                                     <  2:835648/150‑1 (MQ=255)
                                                                             tgtgtTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCCGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCa                                                                   >  2:768537/1‑150 (MQ=255)
                                                                                  tgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACCGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACCTCCAGCAACGcc                                                              >  1:239050/1‑150 (MQ=255)
                                                                                      cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGca                                                          <  1:330948/150‑1 (MQ=255)
                                                                                      cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGca                                                          <  1:153271/150‑1 (MQ=255)
                                                                                          ggCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATc                                                      <  2:1339593/150‑1 (MQ=255)
                                                                                                gAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                                <  1:835943/150‑1 (MQ=255)
                                                                                                 aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                               <  1:1133355/150‑1 (MQ=255)
                                                                                                  ggtCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGaca                                              <  1:643348/148‑1 (MQ=255)
                                                                                                        gCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAAc                                        <  1:1347221/150‑1 (MQ=255)
                                                                                                         cTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACg                                       >  2:1675300/1‑150 (MQ=255)
                                                                                                           cgcgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGAc                                     <  1:1056585/150‑1 (MQ=255)
                                                                                                             cgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACAGGCGTAAAAAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAACGGCACATCCCGCACCGCCAGCACATCCTCTGACCCCAAACGACTg                                   >  2:1319682/1‑150 (MQ=255)
                                                                                                                aCGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCTGCAACGCCAGCACATCCTGTGACACAAAACGACTGTCa                                <  2:178165/150‑1 (MQ=255)
                                                                                                                 cGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGGCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAg                               <  2:1355614/150‑1 (MQ=255)
                                                                                                                 cGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAg                               <  1:856813/150‑1 (MQ=255)
                                                                                                                   ttGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGc                             <  1:1061142/150‑1 (MQ=255)
                                                                                                                       cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCGCCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTCTCAGGCAggg                         >  1:693685/1‑148 (MQ=255)
                                                                                                                       cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTg                         >  1:1460279/1‑150 (MQ=255)
                                                                                                                         tcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGACAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACCTCCTCTGACACACAACGACTGGCAGGCAGTGAt                       >  2:79986/1‑150 (MQ=255)
                                                                                                                          cATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa                      >  2:37462/1‑150 (MQ=255)
                                                                                                                          cATCAATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACAGCGCCGCCAGCACCGGCACATCCAGCAACGACAGCACATCCTCTGACACAAAACGACTGTCACGCAGTGaaa                      >  1:69780/1‑148 (MQ=255)
                                                                                                                          aaTCTATGCCCCACCGAATGCCGGATGCGTGGACCCACTGGCGTAAATAAGGCACCCCTTCTTCGTTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa                      <  2:1769449/149‑1 (MQ=255)
                                                                                                                            tCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAAc                    <  2:1858900/150‑1 (MQ=255)
                                                                                                                             cTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCTCATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACa                   >  2:295337/1‑150 (MQ=255)
                                                                                                                              tATGCCCCAACGAATGCCGTATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCTCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAg                  <  2:159814/150‑1 (MQ=255)
                                                                                                                               aTGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg                 <  2:656630/150‑1 (MQ=255)
                                                                                                                                tGCCCCAACGAATGCCGGATTCGTTGACCCACTAGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGc                <  2:1735841/150‑1 (MQ=255)
                                                                                                                                 gCCCCCACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt               <  1:1803775/150‑1 (MQ=255)
                                                                                                                                 gCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACAGTCAGGCAGTGATAACAGGCt               <  2:1212600/150‑1 (MQ=255)
                                                                                                                                  ggtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              >  2:1634378/4‑150 (MQ=255)
                                                                                                                                  ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              <  2:606643/150‑1 (MQ=255)
                                                                                                                                  ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              <  2:614731/150‑1 (MQ=255)
                                                                                                                                   gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGCCACAAAACGACTGTCGGGCAGTGAGAACAGGCTGa             >  1:1791911/3‑150 (MQ=255)
                                                                                                                                   gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa             >  1:802295/3‑150 (MQ=255)
                                                                                                                                   cccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCCCATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt            >  2:1381978/1‑150 (MQ=255)
                                                                                                                                    ccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt            >  1:77400/1‑150 (MQ=255)
                                                                                                                                     cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATa           <  1:1634394/150‑1 (MQ=255)
                                                                                                                                         gAATGCCGGATTCGTCGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGTCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACg       >  2:16792/1‑150 (MQ=255)
                                                                                                                                           aTGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCt     <  2:1619079/150‑1 (MQ=255)
                                                                                                                                           aTGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGGGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCt     <  1:1728580/150‑1 (MQ=255)
                                                                                                                                              ccGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  <  2:1529521/150‑1 (MQ=255)
                                                                                                                                              ccGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  <  2:1116844/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960785‑2961076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: