Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 578,102 C→T Q164* (CAA→TAA)  rrrD → DLP12 prophage; putative lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,1020CT99.3% 533.9 / NA 149Q164* (CAA→TAA) rrrDDLP12 prophage; putative lysozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (72/76);  minor base G (1/0);  total (73/76)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

GAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGA  >  NC_000913/577956‑578248
                                                                                                                                                  |                                                                                                                                                  
gAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTaa                                                                                                                                                  >  1:601602/1‑149 (MQ=255)
   cATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTg                                                                                                                                             >  1:1711912/1‑151 (MQ=255)
    aTCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGa                                                                                                                                            <  1:1164327/151‑1 (MQ=255)
      cGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAg                                                                                                                                           <  4:1649552/150‑1 (MQ=255)
      cGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAg                                                                                                                                           <  3:245587/150‑1 (MQ=255)
      cGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGc                                                                                                                                          >  3:304963/1‑151 (MQ=255)
       gACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCa                                                                                                                                         <  2:274474/151‑1 (MQ=255)
         cGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCaga                                                                                                                                       <  1:1024719/151‑1 (MQ=255)
         cGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCaga                                                                                                                                       <  1:296999/151‑1 (MQ=255)
           cttcttCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGt                                                                                                                                     >  2:355562/1‑151 (MQ=255)
             tcttcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTa                                                                                                                                    <  4:1872304/150‑1 (MQ=255)
               ttcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAAcc                                                                                                                                 >  1:1825277/1‑151 (MQ=255)
                tcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAAcc                                                                                                                                 <  4:81787/150‑1 (MQ=255)
                tcGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAAcc                                                                                                                                 <  2:56726/150‑1 (MQ=255)
                  gCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACcgcg                                                                                                                              >  1:1500944/1‑151 (MQ=255)
                  gCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACcgcg                                                                                                                              <  2:53879/151‑1 (MQ=255)
                   cAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACcgcg                                                                                                                              >  1:995246/1‑150 (MQ=255)
                   cAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACcgcg                                                                                                                              >  2:465261/1‑150 (MQ=255)
                   cAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACcgcg                                                                                                                              >  1:1537581/1‑150 (MQ=255)
                         aaaCCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTatat                                                                                                                       >  4:12588/1‑151 (MQ=255)
                          aaCCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATc                                                                                                                      >  1:1381904/1‑151 (MQ=255)
                          aaCCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGGCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCCCGATTATATc                                                                                                                      >  1:467638/1‑151 (MQ=255)
                              aGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGct                                                                                                                  >  1:1607544/1‑151 (MQ=255)
                              aGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGct                                                                                                                  <  1:896115/151‑1 (MQ=255)
                               gggCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGctc                                                                                                                 <  2:224126/151‑1 (MQ=255)
                                   gATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGAt                                                                                                             <  1:539359/151‑1 (MQ=255)
                                     tatCAAAGGCGCATGTGATCAGCTGCGTCGCTGGAAATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTACCCGCGATTATATCCGCTCTGAtt                                                                                                            <  3:414601/150‑1 (MQ=255)
                                      atCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACAGCGATTATATCCGCTCTGATTAt                                                                                                          <  2:262625/151‑1 (MQ=255)
                                        cAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCt                                                                                                        >  1:1004804/1‑151 (MQ=255)
                                         aaaGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCt                                                                                                        >  1:1287414/1‑150 (MQ=255)
                                            ggCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGATCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGcat                                                                                                    <  2:143698/151‑1 (MQ=255)
                                              cgcATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGcatca                                                                                                  <  1:693333/151‑1 (MQ=255)
                                                cATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCatcatc                                                                                                <  1:1064158/151‑1 (MQ=255)
                                                 aTGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCg                                                                                               >  4:1710587/1‑151 (MQ=255)
                                                 aTGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCg                                                                                               <  4:4111/151‑1 (MQ=255)
                                                   gtgATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGt                                                                                              <  4:1361051/150‑1 (MQ=255)
                                                   gtgATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTc                                                                                             <  2:361899/151‑1 (MQ=255)
                                                   gtgATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTc                                                                                             >  4:353574/1‑151 (MQ=255)
                                                     gATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCa                                                                                            <  4:24818/150‑1 (MQ=255)
                                                     gATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAg                                                                                           <  4:1602214/151‑1 (MQ=255)
                                                      aTCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGc                                                                                          <  1:1583043/151‑1 (MQ=255)
                                                      aTCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGc                                                                                          <  1:130558/151‑1 (MQ=255)
                                                        cAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCt                                                                                        >  4:747992/1‑151 (MQ=255)
                                                         aGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTg                                                                                       >  4:1877267/1‑151 (MQ=255)
                                                               gTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTc                                                                                     <  1:1209957/147‑1 (MQ=255)
                                                                tCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTTAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATggg                                                                                >  1:802596/1‑151 (MQ=255)
                                                                tCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATggg                                                                                >  1:206552/1‑151 (MQ=255)
                                                                 cGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCAGGGGc                                                                               <  3:456350/151‑1 (MQ=255)
                                                                  gCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTTAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCg                                                                              >  4:232329/1‑151 (MQ=255)
                                                                  gCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCg                                                                              <  1:914749/151‑1 (MQ=255)
                                                                       aCATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCaa                                                                         <  1:1031644/151‑1 (MQ=255)
                                                                         aTACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGCGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGACAGCCTGTCATGGGCGGTCAAt                                                                        >  1:1249472/1‑150 (MQ=255)
                                                                           aCGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCa                                                                      <  4:1304857/150‑1 (MQ=255)
                                                                           aCGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCa                                                                      <  4:778660/150‑1 (MQ=255)
                                                                           aCGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCAt                                                                     <  1:748882/151‑1 (MQ=255)
                                                                           aCGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCAt                                                                     <  4:1847952/151‑1 (MQ=255)
                                                                              cTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTAc                                                                  >  4:1398100/1‑151 (MQ=255)
                                                                                ggCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTAcc                                                                 >  4:1406205/1‑150 (MQ=255)
                                                                                  cGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTg                                                              <  4:1678919/151‑1 (MQ=255)
                                                                                   ggTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGa                                                             <  1:844851/151‑1 (MQ=255)
                                                                                    gTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCAACGTCAGCCTGTCATGGGCGGTCAATCATTACCGt                                                               <  1:1415551/148‑1 (MQ=255)
                                                                                         cAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGc                                                       >  4:1796730/1‑151 (MQ=255)
                                                                                         cAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGc                                                       >  1:544181/1‑151 (MQ=255)
                                                                                          aaTGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCa                                                      <  4:1361105/151‑1 (MQ=255)
                                                                                             ggAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAAt                                                    <  4:1540700/150‑1 (MQ=255)
                                                                                              gAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCg                                                  >  3:350775/1‑151 (MQ=255)
                                                                                              gAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCg                                                  <  1:1408075/151‑1 (MQ=255)
                                                                                                   gCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTAc                                             <  1:1303301/151‑1 (MQ=255)
                                                                                                       gatgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACaaa                                          <  4:1105760/150‑1 (MQ=255)
                                                                                                       gatgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACaa                                           >  4:480795/1‑149 (MQ=255)
                                                                                                         tgaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTc                                       >  1:444275/1‑151 (MQ=255)
                                                                                                          gaCTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTCAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTcc                                      >  1:1073351/1‑151 (MQ=255)
                                                                                                             tcgtcgTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAgc                                   >  4:1194781/1‑151 (MQ=255)
                                                                                                              cgtcgtGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAgc                                   <  1:586723/150‑1 (MQ=255)
                                                                                                                   tGAGCTTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGcgcg                                >  1:137129/1‑148 (MQ=255)
                                                                                                                   tGAGATTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGAc                              >  1:1198137/1‑150 (MQ=255)
                                                                                                                    gagaTTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACa                             >  4:198371/1‑150 (MQ=255)
                                                                                                                      gaTTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACaaa                           >  4:1144308/1‑150 (MQ=255)
                                                                                                                      gaTTGAGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACaaa                           >  4:1071408/1‑150 (MQ=255)
                                                                                                                           aGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGcc                     <  2:187615/151‑1 (MQ=255)
                                                                                                                           aGCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGcc                     <  1:285906/151‑1 (MQ=255)
                                                                                                                            gCGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGcc                     >  1:1928162/1‑150 (MQ=255)
                                                                                                                             cGTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCag                   <  3:540331/151‑1 (MQ=255)
                                                                                                                              gTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGcc                     >  1:1181438/1‑148 (MQ=255)
                                                                                                                              gTGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCaga                  <  4:1489701/151‑1 (MQ=255)
                                                                                                                               tGAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCag                   >  2:403220/1‑149 (MQ=255)
                                                                                                                                gAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCagag                 <  3:344200/150‑1 (MQ=255)
                                                                                                                                gAAGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAgaga                <  4:673249/151‑1 (MQ=255)
                                                                                                                                  aGTCTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAAc              >  4:941880/1‑151 (MQ=255)
                                                                                                                                     cTGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTg            >  4:1452591/1‑150 (MQ=255)
                                                                                                                                      tGTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGa           >  2:211338/1‑150 (MQ=255)
                                                                                                                                       gTTTGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGaa          <  1:995285/150‑1 (MQ=255)
                                                                                                                                        tttGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGaa          <  3:212172/149‑1 (MQ=255)
                                                                                                                                        tttGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGc        >  4:367038/1‑151 (MQ=255)
                                                                                                                                        tttGTGGGGGTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGGCAGCCTGTCATGGGCGGTCAATCATTCCCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAAGGCCAGAGAACTGa           <  3:108609/148‑1 (MQ=255)
                                                                                                                                          tgtgGGGGTAACAGGGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATTGGCGGTCAATCATTACCTGGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGc        <  3:533548/149‑1 (MQ=255)
                                                                                                                                             gggggTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc    <  4:929763/150‑1 (MQ=255)
                                                                                                                                             gggggTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGc    <  4:184315/150‑1 (MQ=255)
                                                                                                                                             gggggTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCg   <  3:465245/151‑1 (MQ=255)
                                                                                                                                             gggggTAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCg   <  4:800843/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAACATCGACGCTTCTTCGCAAAATAAACCAGGGCGATATCAAAGGCGCATGTGATCAGCTGCGTCGCTGGACATACGCTGGCGGTAAGCAATGGAAAGGCCTGATGACTCGTCGTGAGATTGAGCGTGAAGTCTGTTTGTGGGGGCAACAGTGAGCAGAGTAACCGCGATTATATCCGCTCTGATTATCTGCATCATCGTCAGCCTGTCATGGGCGGTCAATCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGA  >  NC_000913/577956‑578248

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑