Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 323,006 T→C G223G (GGT→GGC ykgF → ferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913323,0060TC98.6% 272.1 / NA 73G223G (GGT→GGCykgFferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/44);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCTATGAGGGGCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACC  >  NC_000913/322858‑323153
                                                                                                                                                    |                                                                                                                                                   
gCTATGAGGGGCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCaa                                                                                                                                                   <  1:162998/151‑1 (MQ=255)
    tGAGGGGCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGAATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATgc                                                                                                                                                >  1:217631/1‑150 (MQ=255)
      aGGGGCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATgcgc                                                                                                                                              <  1:234470/150‑1 (MQ=255)
       ggggCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGa                                                                                                                                            <  2:809565/151‑1 (MQ=255)
          gCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAAtg                                                                                                                                         <  2:19427/151‑1 (MQ=255)
                 aCGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCa                                                                                                                                  <  1:981010/151‑1 (MQ=255)
                 aCGCCTGAAGCGATGACCTTATTAAGCTGGAAAAAAATCCGAGTAGATTTCCTCAGTGATGAAATACGTATTAGCGGCTGTAATTTCTCGGTGGCAGAGACCGGTTCTGTATGCTTCGTGTGCAATGAAGGCAGTCCTCGAATGTGTata                                                                                                                                   <  1:423412/150‑3 (MQ=255)
                  cGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCAc                                                                                                                                 <  1:1069171/151‑1 (MQ=255)
                  cGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCAc                                                                                                                                 <  2:1260447/151‑1 (MQ=255)
                   gCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACg                                                                                                                                <  1:550251/151‑1 (MQ=255)
                   gCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACg                                                                                                                                <  1:696923/151‑1 (MQ=255)
                    cccgAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGGATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGATCGGGATGCCTGGAGACCAATGAAGGTAATGCGCGAATGGGGACCACGc                                                                                                                               >  1:328357/4‑151 (MQ=255)
                      tGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTg                                                                                                                             >  2:353219/1‑151 (MQ=255)
                       gAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCt                                                                                                                              <  2:182641/149‑1 (MQ=255)
                       gAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTg                                                                                                                             <  1:1274124/150‑1 (MQ=255)
                       gAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTg                                                                                                                             <  2:454839/150‑1 (MQ=255)
                        aaGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGcc                                                                                                                           >  2:897456/1‑151 (MQ=255)
                            gATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACACGCTGCCTaaa                                                                                                                       >  2:508768/1‑150 (MQ=255)
                                 ccTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAAc                                                                                                                     >  2:907224/1‑148 (MQ=255)
                                 ccTTATTCATCCGGCAAAAAATCCACGAAGATTTCCTCAGTTCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACg                                                                                                                    <  1:177564/149‑1 (MQ=255)
                                  cTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCat                                                                                                                 >  2:713509/1‑151 (MQ=255)
                                    tATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCatat                                                                                                               >  2:472689/1‑151 (MQ=255)
                                    tATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCat                                                                                                                 <  2:66116/149‑1 (MQ=255)
                                     aTTCATCCGGCAAATAATCCGCGAAGATTTCCTCAGTGCTTAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGCTTCGTTCTGTCTGGTGACCAATGAAGGCGATGCGCGAATGTGTCCCACGCTGCCTGAAACGCagat                                                                                                               <  1:178717/150‑1 (MQ=255)
                                     aTTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCatat                                                                                                               <  1:400800/150‑1 (MQ=255)
                                      ttCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTg                                                                                                             >  2:57064/1‑151 (MQ=255)
                                        cATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCa                                                                                                           <  1:383462/151‑1 (MQ=255)
                                         aTCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGc                                                                                                            >  1:156417/1‑149 (MQ=255)
                                                aaaaaaTCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATggg                                                                                                   <  2:724501/151‑1 (MQ=255)
                                                  aaaaTCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGa                                                                                                  >  2:399677/1‑150 (MQ=255)
                                                   aaaTCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAAt                                                                                                >  1:1187547/1‑151 (MQ=255)
                                                   aaaTCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAAt                                                                                                <  1:373900/151‑1 (MQ=255)
                                                       ccGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGa                                                                                             <  2:1173234/150‑1 (MQ=255)
                                                       ccGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAg                                                                                            <  2:818267/151‑1 (MQ=255)
                                                       ccGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAg                                                                                            <  1:947848/151‑1 (MQ=255)
                                                        cgcgAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGc                                                                                           <  1:1165267/151‑1 (MQ=255)
                                                         gcgAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCg                                                                                          <  2:581682/151‑1 (MQ=255)
                                                           gAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGt                                                                                         >  2:434597/1‑150 (MQ=255)
                                                              gATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTg                                                                                     >  1:1200131/1‑151 (MQ=255)
                                                              gATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGGGGCAGAGACCGGTTCGGTATGCCTGGGGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTAtt                                                                                      >  2:411483/1‑150 (MQ=255)
                                                               aTTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGc                                                                                    >  1:986928/1‑151 (MQ=255)
                                                                        gTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGtt                                                                           <  1:1190150/151‑1 (MQ=255)
                                                                             gAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGc                                                                        <  1:547152/149‑1 (MQ=255)
                                                                               aaTAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAg                                                                    <  2:228448/151‑1 (MQ=255)
                                                                                 tAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGt                                                                  >  2:1293020/1‑151 (MQ=255)
                                                                                  aGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGt                                                                  >  2:729365/1‑150 (MQ=255)
                                                                                   ggTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGt                                                                  >  1:12638/1‑149 (MQ=255)
                                                                                   ggTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTa                                                                 <  1:1198325/150‑1 (MQ=255)
                                                                                        tACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGt                                                            >  2:1099607/1‑150 (MQ=255)
                                                                                        tACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGt                                                            >  2:287102/1‑150 (MQ=255)
                                                                                          ccGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTAtt                                                         <  2:802562/151‑1 (MQ=255)
                                                                                               tGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATc                                                     >  2:903150/1‑150 (MQ=255)
                                                                                                gTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATc                                                     <  1:244777/149‑1 (MQ=255)
                                                                                                   aTTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGGCTAAAACGCATATTTCAGTGATGGGAATGGAGCGTATTGAACCCGCGTTTGCCGAGGTAGATGTATAGATCaaaag                                                <  1:255300/151‑5 (MQ=255)
                                                                                                   aTTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCAt                                                <  2:852073/151‑1 (MQ=255)
                                                                                                     ttCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATg                                               <  1:1213222/150‑1 (MQ=255)
                                                                                                       cgcgGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTg                                            <  1:298840/151‑1 (MQ=255)
                                                                                                             ggCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGgcgcgc                                      <  2:428034/151‑1 (MQ=255)
                                                                                                             ggCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGgcgcg                                       <  2:728830/150‑1 (MQ=255)
                                                                                                              gCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCa                                     <  1:890825/151‑1 (MQ=255)
                                                                                                                 gagaCCGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTg                                  <  2:895245/151‑1 (MQ=255)
                                                                                                                     ccGGTTCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGt                              >  2:421532/1‑151 (MQ=255)
                                                                                                                          tCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGt                          <  1:450644/150‑1 (MQ=255)
                                                                                                                          tCGGTATGCCTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTg                         >  2:1028391/1‑151 (MQ=255)
                                                                                                                                  ccTGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTg                 <  2:898945/151‑1 (MQ=255)
                                                                                                                                    tGGTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGAc               >  2:860961/1‑151 (MQ=255)
                                                                                                                                      gTGACCAATGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACgg             <  2:81356/151‑1 (MQ=255)
                                                                                                                                             aTGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATaca       >  2:690401/1‑150 (MQ=255)
                                                                                                                                              tGAAGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATacaac     >  2:822172/1‑151 (MQ=255)
                                                                                                                                                aaGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAacac   >  1:433533/1‑151 (MQ=255)
                                                                                                                                                aaGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAacac   <  1:1292980/151‑1 (MQ=255)
                                                                                                                                                aaGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAacac   <  1:1063097/151‑1 (MQ=255)
                                                                                                                                                 aGGCAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACAcc  >  1:649946/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTATGAGGGGCCGGAAACGCCTGAAGCGATGACCTTATTCATCCGGCAAAAAATCCGCGAAGATTTCCTCAGTGCTGAAATAGGTATTACCGGCTGTAATTTCGCGGTGGCAGAGACCGGTTCGGTATGCCTGGTGACCAATGAAGGTAATGCGCGAATGTGTACCACGCTGCCTAAAACGCATATTGCAGTGATGGGAATGGAGCGTATTGCCCCCACGTTTGCCGAGGTAGATGTATTGATCACCATGCTGGCGCGCAGTGCCGTTGGTGCACGTTTGACGGGATACAACACC  >  NC_000913/322858‑323153

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