Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,316,658 T→C 5.6% L161P (CTG→CCG)  rcsB → DNA‑binding transcriptional activator RcsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,316,6580TC 5.6% 216.2 / 2.1 71L161P (CTG→CCG) rcsBDNA‑binding transcriptional activator RcsB
Reads supporting (aligned to +/- strand):  ref base T (40/27);  new base C (2/2);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGACCGATCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAA  >  NC_000913/2316512‑2316807
                                                                                                                                                  |                                                                                                                                                     
accgaccgATCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGtt                                                                                                                                                    >  1:1348616/1‑150 (MQ=255)
   gaccgaTCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGc                                                                                                                                                 >  1:1476827/1‑150 (MQ=255)
     ccgaTCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCgg                                                                                                                                               >  1:1004577/1‑150 (MQ=255)
       gaTCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCCGTTTGCGGaa                                                                                                                                             >  1:8349/1‑150 (MQ=255)
       gaTCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCCGTTTGCGGaa                                                                                                                                             >  1:1330095/1‑150 (MQ=255)
          cTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGCCAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGc                                                                                                                                          >  1:5910/1‑150 (MQ=255)
           tGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCt                                                                                                                                         <  2:819615/150‑1 (MQ=255)
             ccGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTc                                                                                                                                       >  2:1242812/1‑150 (MQ=255)
                 aaGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGCGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTgg                                                                                                                                   <  1:639822/150‑1 (MQ=255)
                      ctcGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGAcc                                                                                                                              <  2:341729/150‑1 (MQ=255)
                       tcGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCg                                                                                                                             <  1:1314129/150‑1 (MQ=255)
                        cgccgcGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCga                                                                                                                            <  2:967483/150‑1 (MQ=255)
                        cgccgcGCGGGAGAAAGGGAAGAAATTTACCGCGGACAGCGTTGGTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTCCGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCga                                                                                                                            <  1:1438173/150‑1 (MQ=255)
                            gcgcTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATc                                                                                                                        <  1:1579597/150‑1 (MQ=255)
                            gcgcTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATc                                                                                                                        >  1:312457/1‑150 (MQ=255)
                                tGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTa                                                                                                                    <  2:1300690/150‑1 (MQ=255)
                                tGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTa                                                                                                                    >  2:1493675/1‑150 (MQ=255)
                                    gAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTaaaaa                                                                                                                >  2:1147881/1‑150 (MQ=255)
                                    gAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGGTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTaaaaa                                                                                                                <  1:579004/150‑1 (MQ=255)
                                      aaGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGc                                                                                                              >  2:257610/1‑150 (MQ=255)
                                      aaGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGc                                                                                                              <  1:1201328/150‑1 (MQ=255)
                                          gaagaaATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGaa                                                                                                          >  2:458743/1‑150 (MQ=255)
                                          caagaaATTTACCCCGCAAAGCGTTTCTCGCATGTTGGAAAAAATCATTGCTGGTGGTTACGGTGAAAAGCGTCTCTCGCCAAAAGAGAGTCAAGTTCTGAGCATGTTTGCGGAAGGCTTCCTGGTGACCGAGATCCCTAAAAAGCTGaa                                                                                                          >  2:144806/2‑150 (MQ=255)
                                                   tACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTAt                                                                                                 >  2:1211328/1‑150 (MQ=255)
                                                      cccGGAAAGCGTTTCTCGACTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGATAGTGAAGTTCTGCGCATGTTTGAGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAAACGCAGTATTaa                                                                                              >  1:215194/1‑150 (MQ=255)
                                                                    ctcGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGcc                                                                                <  2:863035/150‑1 (MQ=255)
                                                                    ctcGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGcc                                                                                <  2:1581864/150‑1 (MQ=255)
                                                                           gTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAgaagaa                                                                         >  2:456487/1‑150 (MQ=255)
                                                                           gTTGGAAAAAATCAGTGCTGGTGGTTACGGAGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGATATCGCTAAAAAGCTGAACCGCAGTATTAAATCCATCAGTAGACAgaagaa                                                                         <  2:1088102/150‑1 (MQ=255)
                                                                              ggAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATc                                                                      >  1:607351/1‑150 (MQ=255)
                                                                              ggAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATc                                                                      >  2:1046841/1‑150 (MQ=255)
                                                                              ggAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATc                                                                      >  1:258801/1‑150 (MQ=255)
                                                                               gAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCt                                                                     >  1:1625997/1‑150 (MQ=255)
                                                                               gAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCt                                                                     >  1:1299255/1‑150 (MQ=255)
                                                                               gAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCt                                                                     >  2:1470544/1‑150 (MQ=255)
                                                                                aaaaaaTCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTg                                                                    >  2:720651/1‑150 (MQ=255)
                                                                                              ggtggtTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTg                                                      >  1:799540/1‑150 (MQ=255)
                                                                                              ggtggtTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTg                                                      <  2:795046/150‑1 (MQ=255)
                                                                                                tggtTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTggg                                                    >  2:587513/1‑150 (MQ=255)
                                                                                                tggtTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCCGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTggg                                                    <  2:8349/150‑1 (MQ=255)
                                                                                                tggtTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCCGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTggg                                                    <  2:1330095/150‑1 (MQ=255)
                                                                                                  gtTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGgtg                                                  <  1:1156297/150‑1 (MQ=255)
                                                                                                     aCGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCg                                               <  1:249611/150‑1 (MQ=255)
                                                                                                     aCGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCg                                               <  2:258801/150‑1 (MQ=255)
                                                                                                      cGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCga                                              >  1:1179567/1‑150 (MQ=255)
                                                                                                       ggTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAAt                                                                       <  1:170000/124‑1 (MQ=255)
                                                                                                       ggTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAAt                                                                       >  2:170000/1‑124 (MQ=255)
                                                                                                       ggTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCgag                                             <  2:748360/150‑1 (MQ=255)
                                                                                                            cAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGa                                        >  1:1243172/1‑150 (MQ=255)
                                                                                                            cAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGa                                        >  1:547495/1‑150 (MQ=255)
                                                                                                              aGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGata                                      >  1:343251/1‑150 (MQ=255)
                                                                                                              aGCGTCTATCGCCAAAACAGAGTGAAGTTCTCCGCCTGTTAGCGGAAGTCTTCCTGTAGACAAAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTCCGATGATGAAGCTGGGTGTCGAGAACGata                                      >  2:548289/1‑150 (MQ=255)
                                                                                                                   ctctcGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGcc                                 <  1:798551/150‑1 (MQ=255)
                                                                                                                   ctctcGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGcc                                 >  2:1264266/1‑150 (MQ=255)
                                                                                                                    tctcGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCACTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGccc                                >  1:705321/1‑150 (MQ=255)
                                                                                                                     ctcGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCct                               >  2:328159/1‑150 (MQ=255)
                                                                                                                       cGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCctgc                             >  1:1482212/1‑150 (MQ=255)
                                                                                                                       cGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCctgc                             >  2:956770/1‑150 (MQ=255)
                                                                                                                        gCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCctgct                            >  2:597267/1‑150 (MQ=255)
                                                                                                                        gCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCctgct                            >  2:953120/1‑150 (MQ=255)
                                                                                                                        gCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCctgct                            >  2:460034/1‑150 (MQ=255)
                                                                                                                             aaGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAAtt                       >  1:278056/1‑150 (MQ=255)
                                                                                                                                  agTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTAtctc                  <  1:1211328/150‑1 (MQ=255)
                                                                                                                                   gTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTAtctct                 >  1:181694/1‑150 (MQ=255)
                                                                                                                                           cTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTg         >  2:865345/1‑150 (MQ=255)
                                                                                                                                             gcgcCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGAc       <  2:343251/150‑1 (MQ=255)
                                                                                                                                              cgcCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGAcc      <  1:956770/150‑1 (MQ=255)
                                                                                                                                              cgcCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGAcc      <  2:1179567/150‑1 (MQ=255)
                                                                                                                                              cgcCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGAcc      <  1:1264266/150‑1 (MQ=255)
                                                                                                                                                 cTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGTGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTa   <  2:1348616/150‑1 (MQ=255)
                                                                                                                                                 cTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTa   <  2:5910/150‑1 (MQ=255)
                                                                                                                                                  tGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTaa  >  2:535423/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACCGACCGATCTGCCGAAAGCTCTCGCCGCGCTGCAGAAAGGGAAGAAATTTACCCCGGAAAGCGTTTCTCGCCTGTTGGAAAAAATCAGTGCTGGTGGTTACGGTGACAAGCGTCTCTCGCCAAAAGAGAGTGAAGTTCTGCGCCTGTTTGCGGAAGGCTTCCTGGTGACCGAGATCGCTAAAAAGCTGAACCGCAGTATTAAAACCATCAGTAGCCAGAAGAAATCTGCGATGATGAAGCTGGGTGTCGAGAACGATATCGCCCTGCTGAATTATCTCTCTTCAGTGACCTTAA  >  NC_000913/2316512‑2316807

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: ‑