Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,203 T→C 5.9% intergenic (+46/+25) smrB → / ← yfcO putative DNA endonuclease/DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,2030TC 5.9% 253.1 / 5.1 87intergenic (+46/+25)smrB/yfcOputative DNA endonuclease/DUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base C (3/2);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTG  >  NC_000913/2449055‑2449351
                                                                                                                                                    |                                                                                                                                                    
gCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTTCGGGTAATGGCAGGGGCGAATg                                                                                                                                                     >  2:1102411/1‑150 (MQ=255)
    gtgATGGCCTTTCATCAGTCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGcccc                                                                                                                                                 <  1:380169/150‑1 (MQ=255)
     tgATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCcct                                                                                                                                                >  2:651684/1‑150 (MQ=255)
         ggCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGCGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgcca                                                                                                                                            <  2:1706288/150‑2 (MQ=255)
          gCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCATTCGCCcctgccat                                                                                                                                           <  2:1439526/150‑3 (MQ=255)
             tttCATCAGGCACCAAATGAGTATGGCGGTGATGCGGCGTTGGTGGTGTTGATTGAAGTCGAAGAGTGGCCGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  2:878820/150‑1 (MQ=255)
             tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        >  2:1373609/1‑150 (MQ=255)
             tttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGGTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtg                                                                                                                                        <  1:1391844/150‑1 (MQ=255)
              ttCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGAGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGgtgt                                                                                                                                       >  1:1206605/1‑150 (MQ=255)
                 cacaGGCACCAAACGAGTATGGCGGTGATGCGGCGCGGGTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGGTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATc                                                                                                                                    <  1:1210783/148‑1 (MQ=255)
                 atcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGGTGTTGGTGTTGATTGAAGTCGAAGTGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATc                                                                                                                                    <  2:660005/150‑1 (MQ=255)
                  tcaGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCt                                                                                                                                   >  1:323266/1‑150 (MQ=255)
                   caGGCACCAAAAGAGGATGGCGGTGATGCGGCGTTGTTGGTGTTGAGTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCAATTAGGAGTTGAGCGCGTTGCGCCAGGTACGGGGAATGGCAGGGGCGAATGACCCTGCCTTgtgtttggt                                                                                                                                  >  2:1044729/1‑145 (MQ=255)
                        aCCAAAAGAGTATCGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGc                                                                                                                             <  1:1157092/150‑1 (MQ=255)
                          cAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTa                                                                                                                           >  2:726110/1‑150 (MQ=255)
                           aaaaGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTAc                                                                                                                          <  2:766628/150‑1 (MQ=255)
                               gagTAGGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGACGAGGGGCTGCCGCCGGAGTGGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACt                                                                                                                      <  1:1696189/150‑1 (MQ=255)
                                 gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTTCGGGTAATGGCAGGGGCGAATGCCACTGCCTGGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  1:339040/1‑150 (MQ=255)
                                 gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTTGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  2:302842/1‑150 (MQ=255)
                                 gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCa                                                                                                                    >  2:1474133/1‑150 (MQ=255)
                                 gTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAACAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGAGAGTTGCGACAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGGGTATCTTAGGGCTACAACTCa                                                                                                                    >  1:257488/1‑150 (MQ=255)
                                   aTGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCTGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATTGCCCTGCCTGGTGTTTATAAGGGGTTCaaaagctt                                                                                                                  >  1:1530922/1‑144 (MQ=255)
                                   aTGGCGGTGATGCGGCGTGGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAtt                                                                                                                  <  1:1437306/150‑1 (MQ=255)
                                    tGGCGGTGATGAGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCAtta                                                                                                                 >  2:1695773/1‑149 (MQ=255)
                                     ggCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCa                                                                                                                >  1:796237/1‑150 (MQ=255)
                                      gCGGTTAGGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGATGCCCTGAGGAGTTGAGCGCGGTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAg                                                                                                               <  1:1137056/150‑1 (MQ=255)
                                                     ggttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCAGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCg                                                                                                <  1:1082637/149‑1 (MQ=255)
                                                       ttgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGc                                                                                              >  1:1364096/1‑150 (MQ=255)
                                                        tgGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCa                                                                                             <  1:594951/150‑1 (MQ=255)
                                                             ttgattgaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAg                                                                                        <  2:568009/150‑1 (MQ=255)
                                                                 ttgaAGTCGCAGAGTGGCTGCCGCCGGAGTGGCCCTGAGGAGGTTAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATc                                                                                    <  2:578527/150‑1 (MQ=255)
                                                                 tggaAGTCGAAGTGGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGAAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATc                                                                                    <  2:1136530/148‑1 (MQ=255)
                                                                   gaAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATAGa                                                                                  >  1:347253/1‑150 (MQ=255)
                                                                     aggcGAAGAGGGGCTGCAGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGAtt                                                                                <  2:1240172/147‑1 (MQ=255)
                                                                      gTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTg                                                                               >  1:923087/1‑150 (MQ=255)
                                                                        cGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGc                                                                             >  2:410485/1‑150 (MQ=255)
                                                                         gAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAAAGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCCAGTATAATCGATTGGcg                                                                            >  2:1558119/1‑149 (MQ=255)
                                                                         gAAGAGGGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCt                                                                            <  1:1172580/150‑1 (MQ=255)
                                                                           agagTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCCGGGGCGAATGCCCCTGCCTGGTGTATCGAAGGGCTTCAACTCATTCAGTTTGTACTGGAACCGGCAAGTAGCATCGATTGGCTGc                                                                          >  1:555156/1‑150 (MQ=255)
                                                                            gtgTGGCTGCCTGCGGAGGTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCa                                                                         <  1:323536/148‑1 (MQ=255)
                                                                                  cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                   <  1:211182/150‑1 (MQ=255)
                                                                                  cTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgc                                                                   >  1:1161476/1‑150 (MQ=255)
                                                                                   tcccgccgGAGTTTCCATGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATTGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGTGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATActgct                                                                  <  2:1096552/148‑1 (MQ=255)
                                                                                     ccgccgGTGGTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTgg                                                                <  1:222831/150‑1 (MQ=255)
                                                                                       gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTa                                                              >  2:16148/1‑150 (MQ=255)
                                                                                       gccgGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTa                                                              <  1:450653/150‑1 (MQ=255)
                                                                                           gAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCaa                                                          <  2:323266/150‑1 (MQ=255)
                                                                                            aGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAac                                                         >  2:282498/1‑150 (MQ=255)
                                                                                                gCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgc                                                     >  1:1189891/1‑150 (MQ=255)
                                                                                                  ccTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcac                                                   >  2:809512/1‑150 (MQ=255)
                                                                                                   cggaggagGGGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGGATCTAACGGCGACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacg                                                  <  1:839750/148‑1 (MQ=255)
                                                                                                   cTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacg                                                  >  2:1195592/1‑150 (MQ=255)
                                                                                                    tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                    <  2:428443/147‑1 (MQ=255)
                                                                                                    tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgca                                                    >  1:428443/1‑147 (MQ=255)
                                                                                                    tGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACAGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAAcacgcacgc                                                 >  1:1234450/1‑149 (MQ=255)
                                                                                                            tGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATCCCCCTGCCTGGTGTATCTAAGGGCTACATCTCATTCAGTATGTACTGGATCCGGCAAGTACAAACGATTGGCTGCATACTGCAGGTATCAACACGCACTACCATATCg                                         >  1:270423/1‑150 (MQ=255)
                                                                                                             gAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGCCCATATCgt                                        >  1:969397/1‑150 (MQ=255)
                                                                                                               gcgcGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTCCAACTCAGTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCgtgt                                      >  2:1045968/1‑150 (MQ=255)
                                                                                                                 gcgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCAAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGCTTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTaa                                    >  2:391433/1‑150 (MQ=255)
                                                                                                                  cgTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAAc                                   >  2:35563/1‑150 (MQ=255)
                                                                                                                   gTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAAGTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACt                                  >  2:826347/1‑150 (MQ=255)
                                                                                                                   gTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACt                                  <  1:15359/150‑1 (MQ=255)
                                                                                                                    ttGCGCCAGGTACGGGTAATGGCAGGGGCGAAGCCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATACAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTg                                 <  1:882262/150‑1 (MQ=255)
                                                                                                                      gcgcCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTTGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGtt                               <  1:302842/150‑1 (MQ=255)
                                                                                                                          cAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTc                           >  2:1595944/1‑150 (MQ=255)
                                                                                                                             gTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGt                        <  1:1606988/150‑1 (MQ=255)
                                                                                                                             gTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGt                        >  2:682145/1‑150 (MQ=255)
                                                                                                                              tACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGtt                       <  1:1383849/150‑1 (MQ=255)
                                                                                                                                cGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTAATGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATAAACACGCACGACCATATCGTGTAACTGTTTGTCCGTtat                     >  2:30803/1‑150 (MQ=255)
                                                                                                                                cGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCGAAGGGCGACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGCTTGGGTGCATACTGCTGGTTTAAACACGCCCGACCATATCGTGGTACTGTTTGTCCGTgat                     >  2:721788/1‑150 (MQ=255)
                                                                                                                                 gggTAATTGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTtata                    <  2:62021/150‑1 (MQ=255)
                                                                                                                                  ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGCCAAGTAGAAACGATTGGCTGCATACTGCTGGTATCAACACGCACTACCATATCGTGTAACTGTTTGTCCGTTATAg                   >  2:58626/1‑150 (MQ=255)
                                                                                                                                  ggTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAAACGGCAAGTAGAATCGAGTGGCTACATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTTTAg                   >  2:953456/1‑150 (MQ=255)
                                                                                                                                   gTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                  <  2:1364096/150‑1 (MQ=255)
                                                                                                                                   caccaggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAgg                  >  2:960104/6‑150 (MQ=255)
                                                                                                                                    tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTTAGGGCTACAACTCATTCAGTATGTACTGGACCCGGCAAGTAGAATCGATTGGCTGCTTACTGCGGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg                 <  2:260468/150‑1 (MQ=255)
                                                                                                                                    tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg                 <  2:298738/150‑1 (MQ=255)
                                                                                                                                    tAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAggg                 <  1:815747/150‑1 (MQ=255)
                                                                                                                                     aaTGGCAGGGGCGAATTCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATAATGCTGGTATCAACAAGCACGACCATATCGTGTTACTTTTTGTCCGTTATAGGgt                >  2:1146675/1‑150 (MQ=255)
                                                                                                                                      aTGGCAGGGGTGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGgtg               <  1:472023/150‑1 (MQ=255)
                                                                                                                                        ggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGgtgta             >  1:228057/1‑149 (MQ=255)
                                                                                                                                        ggCAGGGGCATTCGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTc             >  1:590357/1‑150 (MQ=255)
                                                                                                                                         gCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCg            <  2:1134323/150‑1 (MQ=255)
                                                                                                                                                cGAATGCCCCTTCCTGGTGTATCTAAGGGCTACAACTCATTCTGTTTGTACTGGAACCTGCATGTTGATTGGATTGGCTGCATAGTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGc     <  2:1417120/150‑1 (MQ=255)
                                                                                                                                                 gAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCg    <  1:649237/150‑1 (MQ=255)
                                                                                                                                                  aaTGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGt   >  1:239050/1‑150 (MQ=255)
                                                                                                                                                   aTGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTg  <  1:1573004/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCATGTGATGGCCTTTCATCAGGCACCAAAAGAGTATGGCGGTGATGCGGCGTTGTTGGTGTTGATTGAAGTCGAAGAGTGGCTGCCGCCGGAGTTGCCCTGAGGAGTTGAGCGCGTTGCGCCAGGTACGGGTAATGGCAGGGGCGAATGCCCCTGCCTGGTGTATCTAAGGGCTACAACTCATTCAGTATGTACTGGAACCGGCAAGTAGAATCGATTGGCTGCATACTGCTGGTATCAACACGCACGACCATATCGTGTAACTGTTTGTCCGTTATAGGGTGTCGATTTGGCGTG  >  NC_000913/2449055‑2449351

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