Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,439 G→A A148A (GCC→GCT hemY ← putative protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,4390GA99.1% 423.9 / NA 108A148A (GCC→GCThemYputative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (49/58);  total (49/59)
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

Only 100 of 111 total aligned reads displayed.

TACGCCTGTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAA  >  NC_000913/3987293‑3987586
                                                                                                                                                  |                                                                                                                                                   
taCGCCTGTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGcc                                                                                                                                                  >  2:233062/1‑150 (MQ=255)
 aCGCCTGTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTgccagcca                                                                                                                                                 <  2:1488623/150‑1 (MQ=255)
 aCGCCTGTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTgccagcca                                                                                                                                                 <  2:1239262/150‑1 (MQ=255)
      tGTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCg                                                                                                                                           >  2:363464/1‑151 (MQ=255)
       gTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCg                                                                                                                                           <  1:230695/150‑1 (MQ=255)
         tCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCt                                                                                                                                         <  4:294283/150‑1 (MQ=255)
          ccGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTg                                                                                                                                        <  2:1558308/150‑1 (MQ=255)
              cAGACGTAATCCTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTgcgc                                                                                                                                     <  1:329580/149‑1 (MQ=255)
                gACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTc                                                                                                                                 >  2:1148236/1‑151 (MQ=255)
                  cGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTcc                                                                                                                                >  3:446211/1‑150 (MQ=255)
                           tCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATgttggt                                                                                                                      >  2:241968/1‑151 (MQ=255)
                           tCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATgttggt                                                                                                                      >  2:601223/1‑151 (MQ=255)
                           tCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATgttggt                                                                                                                      >  4:1153206/1‑151 (MQ=255)
                           tCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATgttggt                                                                                                                      >  3:49198/1‑151 (MQ=255)
                            ccGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATgttggtt                                                                                                                     >  2:1301451/1‑151 (MQ=255)
                               gATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGgc                                                                                                                  >  1:32488/1‑151 (MQ=255)
                                aTGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGgcg                                                                                                                 >  1:427043/1‑151 (MQ=255)
                                aTGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGgcg                                                                                                                 >  2:608897/1‑151 (MQ=255)
                                  ggcgtggcgtAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGcgcg                                                                                                               >  4:577711/1‑151 (MQ=255)
                                    cgttgcttAACTTCCAACAGCTTATCCACGCCGTGGCTTGCAGCATGGTTTTCATTACGGGCTCGTTTCTTTCTTTTACGGGTGAATTCTCCCCGAATGGTTTCCTTTCCGGCCATCTCCTCTCCCCCTACCAGCCGTTGGTCGgcgc                                                                                                                <  2:1598844/148‑1 (MQ=255)
                                    cgttgcgtAACTTCCAGCAGCTTATCCACGCCGTGGCGTGGAGCATGTTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGTCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGcgcg                                                                                                               <  4:1537530/149‑1 (MQ=255)
                                         cgtAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTtca                                                                                                        <  4:1175491/151‑1 (MQ=255)
                                          gtAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTtcat                                                                                                       >  2:569013/1‑151 (MQ=255)
                                          gtAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTtcat                                                                                                       <  4:382087/151‑1 (MQ=255)
                                            aaCTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTtcatca                                                                                                     <  3:428972/151‑1 (MQ=255)
                                             aCTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATcac                                                                                                    >  2:1018728/1‑151 (MQ=255)
                                                tCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATcacc                                                                                                   <  4:1315173/149‑1 (MQ=255)
                                                 ccAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATcaccac                                                                                                 <  2:395470/150‑1 (MQ=255)
                                                 ccAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACg                                                                                                >  4:933798/1‑151 (MQ=255)
                                                   agcagcTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgtt                                                                                              <  3:212429/151‑1 (MQ=255)
                                                     cagcTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgttgt                                                                                            >  3:369471/1‑151 (MQ=255)
                                                      agcTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgttgtt                                                                                           <  3:427061/151‑1 (MQ=255)
                                                        cTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTgc                                                                                         >  2:1554884/1‑151 (MQ=255)
                                                           aTCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGcgcg                                                                                      <  4:1516489/151‑1 (MQ=255)
                                                           aTCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTAAGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGcgcg                                                                                      <  4:837448/151‑1 (MQ=255)
                                                             ccACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGc                                                                                    >  1:229523/1‑151 (MQ=255)
                                                             ccACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGc                                                                                    <  4:106575/151‑1 (MQ=255)
                                                              cACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCt                                                                                   >  2:916623/1‑151 (MQ=255)
                                                                 gCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCg                                                                                >  2:274241/1‑151 (MQ=255)
                                                                 gCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCg                                                                                >  3:37410/1‑151 (MQ=255)
                                                                  ccGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCgg                                                                               >  2:322784/1‑151 (MQ=255)
                                                                       gCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGcc                                                                          <  1:284537/151‑1 (MQ=255)
                                                                       gCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGcc                                                                          >  2:1401678/1‑151 (MQ=255)
                                                                        cGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCa                                                                         <  4:1293535/151‑1 (MQ=255)
                                                                        cGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCa                                                                         <  2:1353984/151‑1 (MQ=255)
                                                                           gcagcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTa                                                                      <  3:401536/151‑1 (MQ=255)
                                                                            cagcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGt                                                                       >  2:10024/1‑149 (MQ=255)
                                                                              gcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGAt                                                                   <  4:610707/151‑1 (MQ=255)
                                                                               caTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGAt                                                                   >  3:239878/1‑150 (MQ=255)
                                                                               caTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATa                                                                  <  1:369454/151‑1 (MQ=255)
                                                                                aTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAg                                                                 <  2:926489/151‑1 (MQ=255)
                                                                                  ggTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGtt                                                               <  2:573105/151‑1 (MQ=255)
                                                                                  ggTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGtt                                                               <  4:1619463/151‑1 (MQ=255)
                                                                                  ggTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGt                                                                <  4:1062701/150‑1 (MQ=255)
                                                                                   cttttCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTc                                                              >  4:744954/2‑151 (MQ=255)
                                                                                    ttttCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTc                                                              >  4:1376201/1‑150 (MQ=255)
                                                                                        cATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTcacca                                                          <  2:1209449/150‑1 (MQ=255)
                                                                                           tACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCg                                                       <  3:265836/150‑1 (MQ=255)
                                                                                             cGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgtt                                                    <  2:454563/151‑1 (MQ=255)
                                                                                                 ccAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgttgtt                                                 >  1:334165/1‑150 (MQ=255)
                                                                                                 ccAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgttgtt                                                 <  3:334182/150‑1 (MQ=255)
                                                                                                  cAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTcc                                               >  2:1258154/1‑151 (MQ=255)
                                                                                                  cAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTcc                                               <  2:1684399/151‑1 (MQ=255)
                                                                                                   aGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTc                                                <  4:170072/149‑1 (MQ=255)
                                                                                                   aGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcg                                              <  4:1019884/151‑1 (MQ=255)
                                                                                                    gTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcgc                                             <  4:612295/151‑1 (MQ=255)
                                                                                                     ttGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcgcg                                            >  4:1114328/1‑151 (MQ=255)
                                                                                                       gCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTg                                          <  4:1149695/151‑1 (MQ=255)
                                                                                                          gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGAt                                        <  1:465088/150‑1 (MQ=255)
                                                                                                          gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATc                                       >  2:307953/1‑151 (MQ=255)
                                                                                                          gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATc                                       >  4:469864/1‑151 (MQ=255)
                                                                                                          gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATc                                       >  3:448490/1‑151 (MQ=255)
                                                                                                           aCGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCg                                      >  4:1320038/1‑151 (MQ=255)
                                                                                                             gTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGc                                    <  3:112813/151‑1 (MQ=255)
                                                                                                              tACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGc                                    >  2:181948/1‑150 (MQ=255)
                                                                                                              tACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGc                                    >  2:892830/1‑150 (MQ=255)
                                                                                                               aCGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCa                                   <  2:576913/150‑1 (MQ=255)
                                                                                                                cgcgGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAtt                                 <  1:253354/151‑1 (MQ=255)
                                                                                                                cgcgGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGAGGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAtt                                 <  2:700403/151‑1 (MQ=255)
                                                                                                                  cgGGTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAtttt                               <  4:47878/151‑1 (MQ=255)
                                                                                                                   gggTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttt                              >  2:1032642/1‑151 (MQ=255)
                                                                                                                   gggTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttt                              <  4:1090/151‑1 (MQ=255)
                                                                                                                   gggTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttt                              >  4:686984/1‑151 (MQ=255)
                                                                                                                    ggTGATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTc                             <  4:384440/151‑1 (MQ=255)
                                                                                                                       gATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGcc                          >  1:164132/1‑151 (MQ=255)
                                                                                                                       gATTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGcc                          >  2:1047313/1‑151 (MQ=255)
                                                                                                                        aTTTCTACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCa                         <  2:916743/151‑1 (MQ=255)
                                                                                                                             tACCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGc                    >  3:296457/1‑151 (MQ=255)
                                                                                                                              aCCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCTCGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCa                      <  1:85186/148‑1 (MQ=255)
                                                                                                                              aCCGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCt                   <  2:1068201/151‑1 (MQ=255)
                                                                                                                               ccGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCtt                  <  4:242354/151‑1 (MQ=255)
                                                                                                                               ccGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCtt                  <  4:364006/151‑1 (MQ=255)
                                                                                                                                cGGAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCttt                 <  4:833771/151‑1 (MQ=255)
                                                                                                                                 agAATGGTGTCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCtttt                >  2:389706/2‑151 (MQ=255)
                                                                                                                                       gtgtCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACtt          <  1:392475/151‑1 (MQ=255)
                                                                                                                                       gtgtCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACtt          <  4:1214475/151‑1 (MQ=255)
                                                                                                                                        tgtCGTTGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTg         >  3:78758/1‑151 (MQ=255)
                                                                                                                                              tGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCa                      >  4:636765/1‑132 (MQ=255)
                                                                                                                                              tGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCa                      <  2:635906/132‑1 (MQ=255)
                                                                                                                                              tGCCAGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGAt    <  4:1238587/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACGCCTGTTCCGCCAGACGTAATACTTCCGGATGGCGTGGCGTAACTTCCAGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAA  >  NC_000913/3987293‑3987586

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