Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,171,417 G→A L850L (CTG→TTG)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,171,4170GA98.5% 202.3 / NA 66L850L (CTG→TTG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/30);  minor base C (0/1);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCGCT  >  NC_000913/1171170‑1171663
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
gCTTCACGACGTTTTTGGGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGCTAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCCTTCATCACCCGTTCCAATtc                                                                                                                                                                                                                                                      >  1:341246/1‑251 (MQ=255)
          gTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGAGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGGTGATGATGGAAATCATTCATCACCCGTTCCAATtcgcgctcgcgc                                                                                                                                                                                                                                            >  2:441654/1‑251 (MQ=255)
           tttttGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCTTCACCCGTTCCAATtcgcgctcgcgc                                                                                                                                                                                                                                            >  2:11592/1‑250 (MQ=255)
             tttGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTCCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATc                                                                                                                                                                                                                                         >  2:454759/1‑251 (MQ=255)
                gtgCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGGCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTg                                                                                                                                                                                                                                       <  2:487623/250‑1 (MQ=255)
                      cGGTAGTAATCGCTTTCGGATGCGGTGTCAGCAACCATGCATATGCATGATGATGCGAACGTCCGACGAGACCGCGTAACTGGTGCAGCTGCGCCGGACCGAAGTGATCCGCGCGTGCAATGATAATAGTGTTGGCTGTCGGGGTGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCCCGTTGAAACGATGATGATGGAAATAATACATAACCAGTTCCCATTCCCGCTCGCGCATGTGCCCGTg                                                                                                                                                                                                                                 >  2:436287/1‑250 (MQ=255)
                       ggTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGAc                                                                                                                                                                                                                               >  2:457377/1‑251 (MQ=255)
                           gTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGAt                                                                                                                                                                                                                           <  1:100749/251‑1 (MQ=255)
                                 gCTTTTGGATGCGGTGGCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGGGCAGCTGCGCCAGAACGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGAt                                                                                                                                                                                                                     <  2:391867/251‑1 (MQ=255)
                                  cTTTTTGAGGCGGTTTCAGCAACCATGCATTTGCCTGATGATGCGAACGTCCGACGCGCCCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATACTGGTTGTACAAACCAGCACGTTGAACCGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGc                                                                                                                                                                                                                        <  2:266250/247‑1 (MQ=255)
                                   ttttGGGTGCGGCGTCAGGAACAATGCATATGCATGATGATGCGATCGTCCGGCGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCCATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAAAGGTTGTACAAACCATCACGTTGAAACGTTGCTGCTGGAAAACATTCATCACCAGTTCCAATTCGCGCTCGCGCCTCTGCCCGTGACCGCTGGCGATcc                                                                                                                                                                                                                   >  2:473913/1‑251 (MQ=255)
                                    tttGGATGCGGTGTCAGCATTCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGGCGGGATGTCGATCCCGGTTGCGATAATGGTTGTACAAACCAGCACGTTGGAACGTTGATGATGGAAATCATGCATCTCCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATc                                                                                                                                                                                                                    <  1:83522/249‑1 (MQ=255)
                                      gggatggggggttaGCAACCAGGCATATGCCTGTTGATGCGAACGTCCGACGCGACCGCATAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAAGAGTGTTGGCTGTAGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGAGGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGACGGCGATCcgcg                                                                                                                                                                                                                <  2:487678/238‑1 (MQ=255)
                                         aTGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCtt                                                                                                                                                                                                             >  1:436136/1‑251 (MQ=255)
                                                 tCAGCAACCATGCCTATGCCTGATGATGCGAACGTCCGACGCAACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGCATTCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCCATGGCGATCCGCGCTTCTGGCAcc                                                                                                                                                                                                     >  2:238026/1‑251 (MQ=255)
                                                     cAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATGCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGtt                                                                                                                                                                                                 <  1:61656/251‑1 (MQ=255)
                                                       aCCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCt                                                                                                                                                                                               <  1:411713/251‑1 (MQ=255)
                                                            gCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAg                                                                                                                                                                                          <  1:454759/251‑1 (MQ=255)
                                                                     tgatgaTGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTAGCGCATCTGCCCGTGACCGATGGCGTTCCGCGCTTCTGGCACCAGTTCTGCAAGCCGTTcggc                                                                                                                                                                                 >  2:499353/1‑251 (MQ=255)
                                                                      tgtgatgCGAACGTCCGACGCGACCCAGTTACTGGTGCAGCGGCGCCAGCCCGAAGTGATCCGCGCGTTCAATGATAATATTGTTGGCTGTCGGGATGTCGAGCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTcggc                                                                                                                                                                                 <  2:285208/248‑1 (MQ=255)
                                                                      gatgatGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTcggcg                                                                                                                                                                                >  1:37009/1‑251 (MQ=255)
                                                                            gCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTc                                                                                                                                                                          <  1:348319/251‑1 (MQ=255)
                                                                            gCGAACGGCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTc                                                                                                                                                                          <  2:222153/251‑1 (MQ=255)
                                                                                    ccgaccggaacgcgTAACTGGTGCAGCTGCGCCAGAACGAAGTGATCCGCGCGGTCAATGATAATAGTGTTTGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGt                                                                                                                                                                   <  1:142297/240‑1 (MQ=255)
                                                                                    ccGACGCGACGTCGTAACTGGTGAAGCTGCGCCAGACAGAAGTGATCGGCGCGTGTAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATACATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGAGGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGtt                                                                                                                                                                  <  1:241734/251‑1 (MQ=255)
                                                                                      cgcgccggtctgTAACTGGGGCAGCTGCGCCAGCCCGAAGGGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGAGGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGtttt                                                                                                                                                                <  2:200989/240‑1 (MQ=255)
                                                                                           gACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCcac                                                                                                                                                            >  1:466457/1‑250 (MQ=255)
                                                                                            cccGCGTAACTAGTGCAGCTGCGCCAGCCCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGAGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCCcat                                                                                                                                                          >  2:242709/2‑251 (MQ=255)
                                                                                            aCCGCGTAACTGGTGCAGCTGCGCCAGACCGACGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAACTCATTCATCAACCGTTCCAATTCGCGCTCGCGCCTCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCCGCCGTTAGGCGGCTTTCTGAATGTTTTCCAcat                                                                                                                                                          >  2:361035/1‑251 (MQ=255)
                                                                                             ccGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCAcatc                                                                                                                                                         <  1:471823/251‑1 (MQ=255)
                                                                                                 gTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCACGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGtttt                                                                                                                                                                >  1:200989/1‑240 (MQ=255)
                                                                                                   aaCTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTa                                                                                                                                                   <  1:457377/251‑1 (MQ=255)
                                                                                                    aCTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTag                                                                                                                                                  <  1:502716/251‑1 (MQ=255)
                                                                                                                  gcCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTg                                                                                                                                    >  1:480764/1‑251 (MQ=255)
                                                                                                                      gACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGtcct                                                                                                                                <  2:108690/251‑1 (MQ=255)
                                                                                                                       aCCGAAGTGATCCGCGCGTTCAATGATAATAGGGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTTTACAATCCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGtcctc                                                                                                                               <  2:370656/251‑1 (MQ=255)
                                                                                                                          gAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCcgc                                                                                                                            >  1:469150/1‑251 (MQ=255)
                                                                                                                          gAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCACGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGGTAATAAACCTTTCCTCcgc                                                                                                                            >  2:371835/1‑251 (MQ=255)
                                                                                                                                aTCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGCTGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTg                                                                                                                                    >  2:502774/1‑237 (MQ=255)
                                                                                                                                aTCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTg                                                                                                                                    <  1:502774/237‑1 (MQ=255)
                                                                                                                                   cgcgcgTTCAATGATAAGAGTGTTGGCTGTCGGGAGGTCGATCCCGGTTTCGATAATGGTTTTACAAACCAGCACGTTGAAACGTTGATGATGTAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAAttt                                                                                                                   <  1:475466/251‑1 (MQ=255)
                                                                                                                                       cgTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACg                                                                                                               >  1:303627/1‑251 (MQ=255)
                                                                                                                                        gTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGt                                                                                                                                   <  2:85004/230‑1 (MQ=255)
                                                                                                                                        gTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGt                                                                                                                                   >  1:85004/1‑230 (MQ=255)
                                                                                                                                             aTGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATGCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGAt                                                                                                         <  1:231228/251‑1 (MQ=255)
                                                                                                                                               gATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGGTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAATCCTGTCCTCCGCACAAAATTTCACGCAGGATCg                                                                                                       >  2:193189/1‑251 (MQ=255)
                                                                                                                                                           ggCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCCTCTGCCCGTGACCGATGGCGATCAGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGCATGTTTTCCACATCATTGTCTATATACTAAACATGTCATCCGCGCAAAATTTCACGCCGGACCGCCTcacgaacc                                                                                            >  2:369602/1‑243 (MQ=255)
                                                                                                                                                                        tCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCCCCATGCTGTCAt                                                                               >  2:400561/1‑250 (MQ=255)
                                                                                                                                                                            tCCCGGTTTCGCTAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTcc                                                                          >  2:329298/1‑250 (MQ=255)
                                                                                                                                                                                     cGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGt                                                                                                                                                    >  2:397214/1‑168 (MQ=255)
                                                                                                                                                                                     cGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGt                                                                                                                                                    <  1:397214/168‑1 (MQ=255)
                                                                                                                                                                                       aTAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGCTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGt                                                              <  1:214013/251‑1 (MQ=255)
                                                                                                                                                                                       aTAATGGTAGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCCATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCAATGTAGAGATAATAACCCTGTCATCCGCGCAAAATTTCACGCTGGATCGCCTCCCGGACCACACTGCTGTCTTACTCACGGACAAAggg                                                              >  2:306973/1‑250 (MQ=255)
                                                                                                                                                                                                   cacaacAGCACGTTGAAACGTTGAGGATGGAAAGCATTCATCTCCCTTTCCATTCCGCACTCGGGTGTCTGCCCGCGTCCGATGGCTATCCTTGTTTGTGGCACCAGGTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATTAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCGCCATGCTGTCATACTCACTGACATAGGTTTTAACTGTCAg                                                  <  2:489112/246‑1 (MQ=255)
                                                                                                                                                                                                     aaCCAGCACGGTGAAACGGTGATGATGGAAATCATTCATCACCCGTTCCAAGTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCTCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAgac                                                <  2:37009/251‑1 (MQ=255)
                                                                                                                                                                                                           cccGTTGAAACGTTTATGATGGACATACTTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATCGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGACTGTATTCCACATCATTGTAGCGATAATAACCCTGTCCTCCGCGCCCAATTTCACGCAGGATCGCCTCCCGGACCACCATTCTATCAGAATCACGACAAAAGGTTTTAACTGCCAGCCGAcggg                                          >  2:371564/3‑251 (MQ=255)
                                                                                                                                                                                                             cGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGACAGACGAcgggcg                                        >  1:504709/1‑251 (MQ=255)
                                                                                                                                                                                                                 gaaacgtggattaggtaaaTCAGTCATCACCCGTTCCAATGCGCGCGCGCGCCTCTGCCCGTGTCCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTGCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAACTGTTTTAACTGCCAGACGACGggcgggcg                                    <  2:357771/235‑1 (MQ=255)
                                                                                                                                                                                                                 gAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACgggcgggc                                     <  1:11592/250‑1 (MQ=255)
                                                                                                                                                                                                                     cGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGt                                >  1:109744/1‑251 (MQ=255)
                                                                                                                                                                                                                             atgGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATc                                                                                                        <  1:87884/172‑1 (MQ=255)
                                                                                                                                                                                                                             atgGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATc                                                                                                        >  2:87884/1‑172 (MQ=255)
                                                                                                                                                                                                                                      attcatCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCCCGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGAAAggg               >  2:484545/1‑250 (MQ=255)
                                                                                                                                                                                                                                            caCCCGTTCCACTTCGCGCTTGCGAATCTGCCCGTGACCGTTGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGATCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACATGTCACGCa         <  1:199725/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cccGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGACCCGCAtt       >  1:299811/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   tCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCGcc  >  1:155911/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GCTTCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCGCT  >  NC_000913/1171170‑1171663

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