Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 414,832 G→T 7.2% H308N (CAT→AAT)  sbcC ← ATP dependent, structure specific DNA nuclease ‑ SbcC subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913414,8320GT 7.2% 410.2 / 13.0 123H308N (CAT→AAT) sbcCATP dependent, structure specific DNA nuclease ‑ SbcC subunit
Reads supporting (aligned to +/- strand):  ref base G (63/51);  new base T (3/6);  total (66/57)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 126 total aligned reads displayed.

GCGGTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTT  >  NC_000913/414683‑414981
                                                                                                                                                     |                                                                                                                                                     
       tgtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgc                                                                                                                                                <  1:102575/150‑1 (MQ=255)
       tgtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgc                                                                                                                                                <  1:853649/150‑1 (MQ=255)
       tgtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAgc                                                                                                                                                >  2:1437467/1‑150 (MQ=255)
        gtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgca                                                                                                                                               <  2:170100/150‑2 (MQ=255)
        gtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAgcg                                                                                                                                               >  1:1314917/1‑150 (MQ=255)
        gtTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAgcg                                                                                                                                               <  1:967189/150‑1 (MQ=255)
           ccTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGcgccg                                                                                                                                            >  2:394245/1‑150 (MQ=255)
            cTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCGGGCGAATATGCGCCAGcgccgc                                                                                                                                           >  2:1440483/1‑150 (MQ=255)
             tGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgcagattg                                                                                                                                          <  2:278343/150‑7 (MQ=255)
             tGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAgcagattg                                                                                                                                          <  1:1335036/150‑7 (MQ=255)
              gTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCgcgc                                                                                                                                         <  1:1515495/150‑1 (MQ=255)
              gTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCgcgc                                                                                                                                         <  1:815071/150‑1 (MQ=255)
                aaCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgcag                                                                                                                                              <  1:320953/143‑3 (MQ=255)
                aaCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGCATATTCGCCAgcag                                                                                                                                              >  2:1077378/1‑141 (MQ=255)
                aaCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                       >  1:778441/1‑150 (MQ=255)
                aaCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                       >  1:221084/1‑150 (MQ=255)
                 aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgt                                                                                                                                      <  1:17883/150‑1 (MQ=255)
                 aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCcgt                                                                                                                                      >  1:304723/1‑150 (MQ=255)
                    aGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGtt                                                                                                                                   >  2:96401/1‑150 (MQ=255)
                       tATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTg                                                                                                                                >  2:1416991/1‑150 (MQ=255)
                        aTTCAGGCTTTGTTGCTGCTGCTGTAATTCTGATGACTGCTTCGCCGCGTGGTGGAGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGc                                                                                                                               >  2:1214825/1‑150 (MQ=255)
                           cAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGAt                                                                                                                            >  1:430429/1‑150 (MQ=255)
                                  tgttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTcc                                                                                                                     >  1:1551771/1‑150 (MQ=255)
                                  tgttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCAATGCGTTcc                                                                                                                     >  1:1109222/1‑150 (MQ=255)
                                   gttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTccc                                                                                                                    <  1:96401/150‑1 (MQ=255)
                                   gttgCTGCGGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTccc                                                                                                                    <  2:804544/150‑1 (MQ=255)
                                     tgctgctgcTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAg                                                                                                                  >  2:420588/1‑150 (MQ=255)
                                     tgctgctgcTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAg                                                                                                                  >  1:1285451/1‑150 (MQ=255)
                                        tgctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtgt                                                                                                               >  1:1719/1‑150 (MQ=255)
                                        tgctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtgt                                                                                                               >  2:574524/1‑150 (MQ=255)
                                          ctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTgg                                                                                                             >  1:1579847/1‑150 (MQ=255)
                                             ctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACg                                                                                                          >  1:650400/1‑150 (MQ=255)
                                             ctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACg                                                                                                          >  2:544056/1‑150 (MQ=255)
                                               gTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGaa                                                                                                        >  1:1201963/1‑150 (MQ=255)
                                               gTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGaa                                                                                                        >  1:437323/1‑150 (MQ=255)
                                                 aaTTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGa                                                                                                      >  2:1476837/1‑150 (MQ=255)
                                                 aaTTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGa                                                                                                      >  1:262596/1‑150 (MQ=255)
                                                      tgctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCg                                                                                                 <  2:487405/150‑1 (MQ=255)
                                                       gctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGt                                                                                                >  1:36558/1‑150 (MQ=255)
                                                       gctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGt                                                                                                >  1:792021/1‑150 (MQ=255)
                                                          gACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGcc                                                                                             <  1:1369900/150‑1 (MQ=255)
                                                            cTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCgg                                                                                           <  1:520442/150‑1 (MQ=255)
                                                              gCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                       <  2:1114293/104‑1 (MQ=255)
                                                              gCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                       >  1:1114293/1‑104 (MQ=255)
                                                              gCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGtt                                                                                         >  1:1257835/1‑150 (MQ=255)
                                                                   gccgcGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGcc                                                                                    >  2:533659/1‑150 (MQ=255)
                                                                   gccgcGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGcc                                                                                    >  1:182808/1‑150 (MQ=255)
                                                                   gccgcGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGcc                                                                                    >  2:1130311/1‑150 (MQ=255)
                                                                   gccgcGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGGTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGcc                                                                                    >  2:16331/1‑150 (MQ=255)
                                                                     cgcgTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAg                                                                                  >  1:330955/1‑150 (MQ=255)
                                                                     cgcgTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCCGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAg                                                                                  <  2:304723/150‑1 (MQ=255)
                                                                         tggtggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTa                                                                              <  2:306713/150‑1 (MQ=255)
                                                                           gtggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAg                                                                            >  1:1025158/1‑150 (MQ=255)
                                                                             ggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAgcg                                                                          <  1:1163010/150‑1 (MQ=255)
                                                                               cGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGcgcc                                                                        >  2:202143/1‑150 (MQ=255)
                                                                                     cTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGt                                                                  >  1:1582678/1‑150 (MQ=255)
                                                                                     cTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGt                                                                  >  1:1222110/1‑150 (MQ=255)
                                                                                     cTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGt                                                                  <  2:142814/150‑1 (MQ=255)
                                                                                       cgcgcgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTg                                                                <  2:27986/150‑1 (MQ=255)
                                                                                        gcgcgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGa                                                               <  2:430429/150‑1 (MQ=255)
                                                                                         cgcgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAg                                                              <  1:59322/150‑1 (MQ=255)
                                                                                         cgcgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTAATTCTTCAATATGCTGTCGAATATTCGCCACCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAg                                                              >  2:822286/1‑150 (MQ=255)
                                                                                          tggAAGAGGCATTTTGCTCTTTAAGAGAGTATTTAGTTCTTCAATCTGCTGGCGAATATGCGAGAGCGCGGCGCTGTGTTCTGCGATGCGTTCCCAGTGTTTACTAATATTTCGTGCCGGTTGTGCCAGACTAAGCTGCGCCTGTTGAgg                                                             <  2:36558/148‑1 (MQ=255)
                                                                                           cgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGt                                                            >  1:345238/1‑150 (MQ=255)
                                                                                               gcgcCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTgcg                                                        >  2:606464/1‑150 (MQ=255)
                                                                                                     ttGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACTAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCttttt                                                  <  2:857843/150‑1 (MQ=255)
                                                                                                      tgtgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGGGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttc                                                 <  2:1352403/150‑1 (MQ=255)
                                                                                                       gttcgctGTATGCGAGTATTTAGTTCTTCAATCTGCTGGCGATTATGCTCCTGCGCCGCGCTGTGTTCTGCGATGTGTTCCCAGTGTGGAGGAATATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGAGCTTTttct                                                <  2:189658/145‑1 (MQ=255)
                                                                                                       gtgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttct                                                <  1:879749/150‑1 (MQ=255)
                                                                                                       gtgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttct                                                <  1:394245/150‑1 (MQ=255)
                                                                                                           tctGTAAGTGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttctt                                            >  2:1560404/1‑150 (MQ=255)
                                                                                                              gTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                         >  2:1187280/1‑150 (MQ=255)
                                                                                                              gTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                         <  1:574524/150‑1 (MQ=255)
                                                                                                               tAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGc                                        <  1:1416991/150‑1 (MQ=255)
                                                                                                                aaGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCt                                       <  2:1201963/150‑1 (MQ=255)
                                                                                                                      gTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGcc                                 >  2:1331412/1‑150 (MQ=255)
                                                                                                                       tATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCt                                >  2:1622903/1‑150 (MQ=255)
                                                                                                                         tttACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGt                              >  1:1084182/1‑150 (MQ=255)
                                                                                                                         tttACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGt                              <  1:1041504/150‑1 (MQ=255)
                                                                                                                           tACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTg                            <  2:1257835/150‑1 (MQ=255)
                                                                                                                            aCTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGc                           <  1:420588/150‑1 (MQ=255)
                                                                                                                             cttcttCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCa                          >  1:1265054/1‑150 (MQ=255)
                                                                                                                             cttcttCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCa                          <  1:48937/150‑1 (MQ=255)
                                                                                                                               tcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAg                        <  1:1331412/150‑1 (MQ=255)
                                                                                                                                cttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                       >  2:1034249/1‑150 (MQ=255)
                                                                                                                                cttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                       <  2:858865/150‑1 (MQ=255)
                                                                                                                                cttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                       >  2:460150/1‑150 (MQ=255)
                                                                                                                                cttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCAATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAgg                       <  2:1109222/150‑1 (MQ=255)
                                                                                                                                 ttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGc                      <  2:650400/150‑1 (MQ=255)
                                                                                                                                   cAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCct                    >  2:187188/1‑150 (MQ=255)
                                                                                                                                    aaTCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctg                   >  2:336180/1‑150 (MQ=255)
                                                                                                                                      tCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgct                 <  1:1130311/150‑1 (MQ=255)
                                                                                                                                       ctgctgGCGAATATGGGCCAGCGCCGCGCTGGGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGGTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTTCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgctg                <  2:849097/150‑1 (MQ=255)
                                                                                                                                       ctgctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgctg                <  1:1011533/150‑1 (MQ=255)
                                                                                                                                       ctgctgGCGAATATGCGCCAGCGCCGCGCTATTTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGACAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACTAATGTTGCAAGGCctgctg                >  2:371493/1‑150 (MQ=255)
                                                                                                                                           tgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCg                                          >  1:1202457/1‑120 (MQ=255)
                                                                                                                                           tgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCg                                          <  2:1202457/120‑1 (MQ=255)
                                                                                                                                                   atGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctggc    >  2:1153635/1‑150 (MQ=255)
                                                                                                                                                    tGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctggct   >  1:775686/1‑150 (MQ=255)
                                                                                                                                                    tGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctggct   >  1:781299/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGGTCGTGTTCCTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTT  >  NC_000913/414683‑414981

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