Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 785,828 C→T P66S (CCT→TCT)  aroG → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913785,8280CT99.3% 565.0 / NA 142P66S (CCT→TCT) aroG3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, phenylalanine repressible
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (76/65);  minor base G (0/1);  total (76/66)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TTCCTCCTGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATCCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATA  >  NC_000913/785682‑785975
                                                                                                                                                  |                                                                                                                                                   
ttCCTCCTGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGt                                                                                                                                                 >  1:833071/1‑151 (MQ=255)
ttCCTCCTGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGt                                                                                                                                                 >  1:1902518/1‑151 (MQ=255)
  cctcctGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTg                                                                                                                                                  <  1:1396888/148‑1 (MQ=255)
  cctcctGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTcg                                                                                                                                               <  4:768791/151‑1 (MQ=255)
    tcctGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTcgcg                                                                                                                                             >  3:24918/1‑151 (MQ=255)
     cctGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCgg                                                                                                                                            >  4:1783536/1‑151 (MQ=255)
     cctGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCgg                                                                                                                                            <  1:1841192/151‑1 (MQ=255)
       tGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGc                                                                                                                                           >  1:888382/1‑150 (MQ=255)
          cGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCaaaa                                                                                                                                       <  4:1295636/151‑1 (MQ=255)
           tcATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCaaa                                                                                                                                        >  4:530448/2‑149 (MQ=255)
           gCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAag                                                                                                                                      <  1:1047432/151‑1 (MQ=255)
            cATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCaaa                                                                                                                                        >  3:490233/1‑148 (MQ=255)
               tGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGt                                                                                                                                   <  4:1002152/150‑1 (MQ=255)
                gCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTa                                                                                                                                  <  2:448782/150‑1 (MQ=255)
                gCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTa                                                                                                                                  <  1:600478/150‑1 (MQ=255)
                gCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTAt                                                                                                                                 <  2:46777/151‑1 (MQ=255)
                    gAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCa                                                                                                                             <  4:1474767/151‑1 (MQ=255)
                       aaaTTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTc                                                                                                                          >  1:239440/1‑151 (MQ=255)
                        aaTTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCg                                                                                                                         >  4:1721319/1‑151 (MQ=255)
                        aaTTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCg                                                                                                                         >  1:239617/1‑151 (MQ=255)
                         aTTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGc                                                                                                                        >  2:391312/1‑151 (MQ=255)
                           tCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCtt                                                                                                                      >  1:312621/1‑151 (MQ=255)
                            ccccGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCtt                                                                                                                      <  4:819809/150‑1 (MQ=255)
                            ccccGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCtt                                                                                                                      <  4:697462/150‑1 (MQ=255)
                            ccccGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTg                                                                                                                     <  3:464452/151‑1 (MQ=255)
                               cGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTg                                                                                                                  >  1:1588261/1‑151 (MQ=255)
                                gCTACTGAAAATGCCGCGAATGCGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTTGCCCATGCTCAATTCATGATGCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGc                                                                                                                    <  1:618125/148‑1 (MQ=255)
                                 cTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTg                                                                                                                  >  1:1470077/1‑149 (MQ=255)
                                    cTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCt                                                                                                             >  4:763517/1‑151 (MQ=255)
                                     tGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTg                                                                                                            >  4:780224/1‑151 (MQ=255)
                                      gAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCt                                                                                                             <  1:1215351/149‑1 (MQ=255)
                                         aaTGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTg                                                                                                        <  1:1060279/151‑1 (MQ=255)
                                           tGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGa                                                                                                       >  1:225619/1‑150 (MQ=255)
                                           tGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGa                                                                                                       >  4:1611283/1‑150 (MQ=255)
                                             ccGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAaga                                                                                                    <  4:1579784/151‑1 (MQ=255)
                                               gcgAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGc                                                                                                  >  1:1899714/1‑151 (MQ=255)
                                                 gAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTg                                                                                                >  4:224004/1‑151 (MQ=255)
                                                  aaTACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTg                                                                                                >  3:9309/1‑150 (MQ=255)
                                                    tACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGaa                                                                                              >  4:1040911/1‑150 (MQ=255)
                                                     aCGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGaa                                                                                              >  1:1872421/1‑149 (MQ=255)
                                                      cGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGaa                                                                                              >  1:1641722/1‑148 (MQ=255)
                                                          tGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCCCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATg                                                                                         >  4:1450835/1‑149 (MQ=255)
                                                                tGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGa                                                                                  <  3:27414/150‑1 (MQ=255)
                                                                tGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGa                                                                                  >  2:27414/1‑150 (MQ=255)
                                                                 gCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGaa                                                                                 <  1:250552/150‑1 (MQ=255)
                                                                 gCCCGAAAAGCGATCCATAAGATCCTGAAAGGGAATGATGATCGCCTGTTGGTTGTGATTGCCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGaa                                                                                 >  4:139909/1‑150 (MQ=255)
                                                                  cccGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAAt                                                                               >  2:34196/1‑151 (MQ=255)
                                                                  cccGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAAt                                                                               >  2:97448/1‑151 (MQ=255)
                                                                   ccGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATc                                                                              <  2:393513/151‑1 (MQ=255)
                                                                     gAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGGTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGt                                                                            <  4:683491/151‑1 (MQ=255)
                                                                      aaaaGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGt                                                                            >  1:479639/1‑150 (MQ=255)
                                                                        aaGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAAt                                                                         >  3:111559/1‑151 (MQ=255)
                                                                         aGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATg                                                                        <  1:1703959/151‑1 (MQ=255)
                                                                         aGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATg                                                                        >  1:1737049/1‑151 (MQ=255)
                                                                           cGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATgcg                                                                      <  1:1260893/151‑1 (MQ=255)
                                                                            gATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATgcgc                                                                     >  4:1223426/1‑151 (MQ=255)
                                                                            gATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATgcgc                                                                     <  1:939822/151‑1 (MQ=255)
                                                                             aTCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGcgcg                                                                    >  2:249533/1‑151 (MQ=255)
                                                                             aTCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGcgcg                                                                    <  3:340822/151‑1 (MQ=255)
                                                                               ccATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTc                                                                  >  4:1318109/1‑151 (MQ=255)
                                                                                    aGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTAttt                                                             >  3:449253/1‑151 (MQ=255)
                                                                                    aGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTAttt                                                             <  4:833873/151‑1 (MQ=255)
                                                                                     gATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTAtttt                                                            >  1:1423361/1‑151 (MQ=255)
                                                                                     gATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTAtttt                                                            >  3:518469/1‑151 (MQ=255)
                                                                                     gATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTAtttt                                                            <  1:702649/151‑1 (MQ=255)
                                                                                        ccTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGa                                                          <  1:1285721/150‑1 (MQ=255)
                                                                                         cTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGaa                                                         >  1:1312537/1‑150 (MQ=255)
                                                                                           gAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGaaa                                                        >  1:226402/1‑149 (MQ=255)
                                                                                            aaaGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAACTGGAAATCGTAATGCGCGTCTATTTTGaaaatc                                                     <  4:189655/151‑3 (MQ=255)
                                                                                              aGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGcca                                                   <  3:21157/151‑2 (MQ=255)
                                                                                              aGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCcg                                                   >  4:1737178/1‑151 (MQ=255)
                                                                                              aGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCcg                                                   <  4:587308/151‑1 (MQ=255)
                                                                                              aGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCcg                                                   <  4:977521/151‑1 (MQ=255)
                                                                                                gTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCcgc                                                  <  3:269167/150‑1 (MQ=255)
                                                                                                 tAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGTGTCTATTTTGAAAAGCCGCGt                                                >  3:267577/1‑151 (MQ=255)
                                                                                                  aaTGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGt                                                >  1:468610/1‑150 (MQ=255)
                                                                                                  aaTGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTa                                               <  2:195319/151‑1 (MQ=255)
                                                                                                    tgatgaTCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTAcc                                             >  4:111516/1‑151 (MQ=255)
                                                                                                     gatgatCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCTCTGCGTTAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTAACa                                            <  1:1865940/151‑1 (MQ=255)
                                                                                                        gatCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACgg                                         >  4:1909629/1‑151 (MQ=255)
                                                                                                         atCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGt                                        <  1:1220057/151‑1 (MQ=255)
                                                                                                              cTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCt                                   >  1:721042/1‑151 (MQ=255)
                                                                                                                gttggttgTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTgg                                 <  2:9309/151‑1 (MQ=255)
                                                                                                                   ggttgTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGaaa                              <  4:239242/151‑1 (MQ=255)
                                                                                                                     ttgTGATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGGACCACGGTGGGCTGGAAAg                             >  1:938496/1‑150 (MQ=255)
                                                                                                                        tgATTGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCt                         >  3:64139/1‑151 (MQ=255)
                                                                                                                          aTTGGCCCATGCTAAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCAGTCGCTTGCTGGCGCTGCGTGAAGAGCTGAGAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGAGTACCACGGTGGGCTGGAAAGGGCt                         <  4:963041/149‑1 (MQ=255)
                                                                                                                            tGGCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGCAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGAtt                     >  3:393799/1‑151 (MQ=255)
                                                                                                                             ggCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGAtt                     >  3:398094/1‑150 (MQ=255)
                                                                                                                             ggCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTa                    <  4:1175515/151‑1 (MQ=255)
                                                                                                                             ggCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTa                    <  1:1025781/151‑1 (MQ=255)
                                                                                                                              gCCCATGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTa                    >  2:199527/1‑150 (MQ=255)
                                                                                                                                  aTGCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGAt               >  1:329713/1‑151 (MQ=255)
                                                                                                                                    gCTCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATcc             <  4:13766/151‑1 (MQ=255)
                                                                                                                                      tCAATTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGc           >  1:1457267/1‑151 (MQ=255)
                                                                                                                                        aaTTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCat         >  3:319583/1‑151 (MQ=255)
                                                                                                                                        aaTTCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCat         <  1:845335/151‑1 (MQ=255)
                                                                                                                                           tCATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATg      >  1:63248/1‑151 (MQ=255)
                                                                                                                                            cATGATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATgg     <  4:225691/151‑1 (MQ=255)
                                                                                                                                               gATTCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGat   <  1:1782836/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTCCTCCTGTCGCATTGCTGGAAAAATTCCCCGCTACTGAAAATGCCGCGAATACGGTTGCCCATGCCCGAAAAGCGATCCATAAGATCCTGAAAGGTAATGATGATCGCCTGTTGGTTGTGATTGGCCCATGCTCAATTCATGATCCTGTCGCGGCAAAAGAGTATGCCACTCGCTTGCTGGCGCTGCGTGAAGAGCTGAAAGATGAGCTGGAAATCGTAATGCGCGTCTATTTTGAAAAGCCGCGTACCACGGTGGGCTGGAAAGGGCTGATTAACGATCCGCATATGGATA  >  NC_000913/785682‑785975

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