Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,880,440 T→A 5.1% D261V (GAC→GTC)  recF ← gap repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,880,4400TA 5.1% 284.4 / 4.5 78D261V (GAC→GTC) recFgap repair protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base A (2/2);  total (40/38)
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.34e-01

CAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAA  >  NC_000913/3880297‑3880582
                                                                                                                                               |                                                                                                                                              
caAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgc                                                                                                                                            >  2:1113459/1‑148 (MQ=255)
     tCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGaa                                                                                                                                       <  1:1113459/148‑1 (MQ=255)
              aaaTCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGt                                                                                                                              >  3:567379/1‑148 (MQ=255)
                atcatcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATg                                                                                                                                     <  4:257160/139‑1 (MQ=255)
                atcatcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATg                                                                                                                                     >  3:257160/1‑139 (MQ=255)
                    tcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATg                                                                                                                                     >  3:755341/1‑135 (MQ=255)
                    tcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATg                                                                                                                                     <  4:755341/135‑1 (MQ=255)
                                     ccgccgCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCg                                                                                                             >  3:320880/1‑142 (MQ=255)
                                     ccgccgCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCg                                                                                                             <  4:320880/142‑1 (MQ=255)
                                     ccgccgCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGcgcg                                                                                                       >  3:34516/1‑148 (MQ=255)
                                             cgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGt                                                                                               <  4:34516/148‑1 (MQ=255)
                                                  tCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCt                                                                                          <  4:567379/148‑1 (MQ=255)
                                                     cGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtg                                                                                                                      <  2:195181/117‑1 (MQ=255)
                                                     cGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGACGCGAATGCGTAAGTCCGCTTtg                                                                                                                      >  1:195181/1‑117 (MQ=255)
                                                           aGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATcgc                                                                                 >  2:1030360/1‑148 (MQ=255)
                                                            ggAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 <  3:912148/99‑1 (MQ=255)
                                                            ggAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 >  4:912148/1‑99 (MQ=255)
                                                             gAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 <  4:276820/98‑1 (MQ=255)
                                                             gAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 >  3:276820/1‑98 (MQ=255)
                                                              aaCTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAg                                                                                                 <  1:722937/129‑1 (MQ=255)
                                                              aaCTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAg                                                                                                 >  2:722937/1‑129 (MQ=255)
                                                              aaCTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGc                                                                                           >  1:856894/1‑135 (MQ=255)
                                                              aaCTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGc                                                                                           <  2:856894/135‑1 (MQ=255)
                                                                ctctCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCa                                                                            >  1:348336/1‑148 (MQ=255)
                                                                              gACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAg                                                              >  1:6576/1‑148 (MQ=255)
                                                                                       cgcACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGc                                                     <  2:348336/148‑1 (MQ=255)
                                                                                       cgcACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGc                                                     <  1:1030360/148‑1 (MQ=255)
                                                                                         cacaTCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAAtt                                                                       >  4:283519/1‑128 (MQ=255)
                                                                                         cacaTCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAAtt                                                                       <  3:283519/128‑1 (MQ=255)
                                                                                                aCAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgt                                                                                                     <  1:992935/91‑1 (MQ=255)
                                                                                                aCAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgt                                                                                                     >  2:992935/1‑91 (MQ=255)
                                                                                                 cAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 <  2:160013/62‑1 (MQ=255)
                                                                                                 cAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 >  1:160013/1‑62 (MQ=255)
                                                                                                 cAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCAc                                                           >  1:1074756/1‑132 (MQ=255)
                                                                                                 cAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCAc                                                           <  2:1074756/132‑1 (MQ=255)
                                                                                                    cTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 <  2:588668/59‑1 (MQ=255)
                                                                                                    cTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTa                                                                                                                                 >  1:588668/1‑59 (MQ=255)
                                                                                                    cTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCa                                                                            >  3:426902/1‑112 (MQ=255)
                                                                                                    cTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCa                                                                            <  4:426902/112‑1 (MQ=255)
                                                                                                     ttAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAg                                                                                                 >  1:990803/1‑90 (MQ=255)
                                                                                                     ttAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAg                                                                                                 <  2:990803/90‑1 (MQ=255)
                                                                                                     ttAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCtt                                       >  3:638772/1‑148 (MQ=255)
                                                                                                        aGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTgc                                                                                                          <  4:399472/78‑1 (MQ=255)
                                                                                                        aGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTgc                                                                                                          >  3:399472/1‑78 (MQ=255)
                                                                                                         gCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTc                                                        >  1:350166/1‑127 (MQ=255)
                                                                                                         gCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTc                                                        <  2:350166/127‑1 (MQ=255)
                                                                                                           tGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTccc                                 >  1:608428/1‑148 (MQ=255)
                                                                                                           tGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTccc                                 <  2:608428/148‑1 (MQ=255)
                                                                                                              ccACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGc                              >  4:633699/1‑148 (MQ=255)
                                                                                                               cACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTc                                                                 >  1:766201/1‑112 (MQ=255)
                                                                                                               cACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTc                                                                 <  2:766201/112‑1 (MQ=255)
                                                                                                               cACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGcc                             >  1:17793/1‑148 (MQ=255)
                                                                                                                 cgcgATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCcgc                           >  1:24778/1‑148 (MQ=255)
                                                                                                                     aTAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaaa                                                                          <  3:223910/97‑1 (MQ=255)
                                                                                                                     aTAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaaa                                                                          >  4:223910/1‑97 (MQ=255)
                                                                                                                        aGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGaa                    <  2:6576/148‑1 (MQ=255)
                                                                                                                        aGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGaa                    <  3:633699/148‑1 (MQ=255)
                                                                                                                         ggTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGaaa                   <  2:24778/148‑1 (MQ=255)
                                                                                                                          gTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAag                  <  4:638772/148‑1 (MQ=255)
                                                                                                                           tATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTcc                                                                >  2:616671/1‑101 (MQ=255)
                                                                                                                           tATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTcc                                                                <  1:616671/101‑1 (MQ=255)
                                                                                                                           tATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTc                                               <  3:737048/118‑1 (MQ=255)
                                                                                                                           tATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAACTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTc                                               >  4:737048/1‑118 (MQ=255)
                                                                                                                              cTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTa                                                                      <  2:1033883/92‑1 (MQ=255)
                                                                                                                              cTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTa                                                                      >  1:1033883/1‑92 (MQ=255)
                                                                                                                                tCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTg                                             >  4:654140/1‑115 (MQ=255)
                                                                                                                                tCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTg                                             <  3:654140/115‑1 (MQ=255)
                                                                                                                                tCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAg            >  3:452717/1‑148 (MQ=255)
                                                                                                                                     cGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCAc                                                           >  4:11502/1‑96 (MQ=255)
                                                                                                                                     cGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCAc                                                           <  3:11502/96‑1 (MQ=255)
                                                                                                                                      ttcgcgCCGACGGTGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGaa              <  3:106703/138‑1 (MQ=255)
                                                                                                                                      ttcgcgCCGACGGTGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCagag      >  4:850044/3‑148 (MQ=255)
                                                                                                                                        cgcgccGACGGTGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTg                                             <  3:98823/103‑1 (MQ=255)
                                                                                                                                        cgcgCCGACGGTGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTg                                             >  4:98823/1‑107 (MQ=255)
                                                                                                                                         gcACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGaaa   <  4:452717/148‑1 (MQ=255)
                                                                                                                                          cACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGaaaa  >  3:693557/1‑148 (MQ=255)
                                                                                                                                              gTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTtctc                                   <  2:987397/111‑1 (MQ=255)
                                                                                                                                              gTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTtctc                                   >  1:987397/1‑111 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAA  >  NC_000913/3880297‑3880582

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