Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,051 A→C 5.3% intergenic (+19/+77) trpT → / ← hdfR tRNA‑Trp/flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,0510AC 5.3% 391.6 / 10.4 133intergenic (+19/+77)trpT/hdfRtRNA‑Trp/flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (57/69);  new base C (4/3);  total (61/72)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

Only 100 of 134 total aligned reads displayed.

CCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAG  >  NC_000913/3946903‑3947199
                                                                                                                                                    |                                                                                                                                                    
cctgTCACGCAGGGGGTCGCGGGGTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAg                                                                                                                                                     <  2:84059/150‑1 (MQ=255)
      aCGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAgg                                                                                                                                               <  1:878946/150‑1 (MQ=255)
        gCAGGGGGTCGCGGTGTCGAGTCCCGTCCGTTCCGCCACCCTACTTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTGCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGAt                                                                                                                                             <  1:272053/150‑1 (MQ=255)
        gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGAt                                                                                                                                             <  1:728949/150‑1 (MQ=255)
         cAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCACCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGAtg                                                                                                                                            <  2:1664149/150‑2 (MQ=39)
            ggggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGATATCATCCTTAGCGAAAGCTAAGGAttttt                                                                                                                                         >  1:108021/1‑150 (MQ=255)
             ggggCGCGGGCTCGAGTCCCGTCCGTACCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTTGGAGTTCGAGTCTCTCCTCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                        <  1:64140/150‑1 (MQ=255)
             gggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTCGAGCACCGGGCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                        <  1:925512/150‑1 (MQ=255)
             gggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                        <  2:1559855/150‑1 (MQ=255)
             gggTCGCGGGTTCGAGTCCCGGCCGTGCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTGCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                        <  2:1271335/150‑1 (MQ=255)
             gggTCGCGGGTGCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGGCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                                                                                                        <  2:2041360/150‑1 (MQ=255)
                gcgcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                                                                                                     <  1:1714622/149‑1 (MQ=255)
                  gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATc                                                                                                                                   >  1:826439/1‑150 (MQ=255)
                  gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATc                                                                                                                                   <  1:1907823/150‑1 (MQ=255)
                  gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATc                                                                                                                                   <  1:961047/150‑1 (MQ=255)
                   cgGGTTCGAGTACCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTATCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCt                                                                                                                                  >  1:1306769/1‑150 (MQ=255)
                     ggTTCGAGTCCCGTCCGAACCGCCACCCTAATTCGGGGCGTAGTTCATTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGCGAGTTCGAGGCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGa                                                                                                                                <  1:861781/150‑1 (MQ=255)
                         cGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATa                                                                                                                            >  1:778966/1‑150 (MQ=255)
                           aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAAc                                                                                                                          >  2:2016652/1‑150 (MQ=255)
                              cccGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCct                                                                                                                       >  1:621604/1‑150 (MQ=255)
                              cccGTCCGTTCCGCCACCCTAATTAGGGGAGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCAGCCCCTGCCAGATATCATCCTTAGCGAAAGCTAAGGATATTTTTTATCTGAAATGACCct                                                                                                                       >  2:387781/1‑150 (MQ=255)
                               ccgtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCctc                                                                                                                      >  1:557686/1‑150 (MQ=255)
                                    cgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGa                                                                                                                 <  2:243029/150‑1 (MQ=255)
                                      ttCCGCCACCCTAATTAGGTGGGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCAGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAg                                                                                                               <  2:1197391/150‑1 (MQ=255)
                                      ttCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAg                                                                                                               <  2:374640/150‑1 (MQ=255)
                                      ttCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAg                                                                                                               <  2:1848844/150‑1 (MQ=255)
                                        ccGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTa                                                                                                             <  2:431940/150‑1 (MQ=255)
                                         cGCCACACTAATTAGTGGCGTAGTTCATTTGGTAGAGCACCGGGCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTaa                                                                                                            <  2:622911/150‑1 (MQ=255)
                                           ccAGCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAAt                                                                                                          <  1:581448/150‑1 (MQ=255)
                                           ccACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAAt                                                                                                          <  1:479461/150‑1 (MQ=255)
                                              cccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                       <  1:1004143/150‑1 (MQ=255)
                                              cccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                       <  2:495399/150‑1 (MQ=255)
                                                    ttAGGGGGGTAGTTTAATTGGGAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTTGCGAAATCTAAGGATTTTTTTTATCTGAATTAACCCTCTCCGAAGTAAATCCTTCTAcc                                                                                                 <  1:1271590/150‑1 (MQ=255)
                                                    ttAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAcc                                                                                                 <  2:2022050/150‑1 (MQ=255)
                                                     tAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCg                                                                                                >  1:1859660/1‑150 (MQ=255)
                                                      aGGTGCGTAGTTCAATTGGTAGAGCACCGGGCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                                                                                               <  2:1462262/150‑1 (MQ=255)
                                                      aGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                                                                                               <  1:1955635/150‑1 (MQ=255)
                                                        gggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGTATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCa                                                                                             <  2:679643/150‑1 (MQ=255)
                                                        gggCGGAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGTGAGTTCGAGTCTCTCCGACCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTATTAATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCa                                                                                             <  1:287758/150‑1 (MQ=255)
                                                         ggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCGTTCTACCGGCAt                                                                                            >  2:109570/1‑150 (MQ=255)
                                                          gCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATc                                                                                           <  2:615492/150‑1 (MQ=255)
                                                           cGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATcc                                                                                          <  2:529189/150‑1 (MQ=255)
                                                           cGTAGTTCAATTGGTAGAGCACCGGTATCCAAAACCGGGTGGTGGTAGTACGGGTCTCTCCGCCCCAGCAAGAAATCATCCTTAGCGAAAGCTAAGGTTTTTTTTTATCTGAAATACCCCTCTCCGAAGTAAATCCTTCTACCGGCATcc                                                                                          >  1:1014564/1‑150 (MQ=255)
                                                               gTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgc                                                                                      >  1:1190647/1‑150 (MQ=255)
                                                                ttCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcc                                                                                     <  2:1049899/150‑1 (MQ=255)
                                                                 tCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcca                                                                                    >  1:1143494/1‑150 (MQ=255)
                                                                  cAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccag                                                                                   >  1:1637172/1‑150 (MQ=255)
                                                                   aaTTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc                                                                                  <  1:662314/150‑1 (MQ=255)
                                                                        gTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTc                                                                             <  1:634485/150‑1 (MQ=255)
                                                                           gagCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCata                                                                          <  2:1973524/150‑1 (MQ=255)
                                                                            agCACCGGTCTCCAAAACCGGGTGTTGTTAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                         <  1:1247444/150‑1 (MQ=255)
                                                                            agCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                         <  1:760850/150‑1 (MQ=255)
                                                                             gCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGCCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATAtt                                                                        <  2:1884833/150‑1 (MQ=255)
                                                                               aCCGGTCTCCAAAACCGGGTTTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTaa                                                                      <  1:1882125/150‑1 (MQ=255)
                                                                                ccGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                     >  2:1305007/1‑150 (MQ=255)
                                                                                 cGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCAGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGCAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATa                                                                    <  1:98190/150‑1 (MQ=255)
                                                                                  ggTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATac                                                                   <  1:379362/150‑1 (MQ=255)
                                                                                   gTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTAGCCAGCCATTCATATTAATaca                                                                  <  1:340880/150‑1 (MQ=255)
                                                                                      tcCAAAACCGGGTGTTGGGAGTTCGAGGCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACtt                                                               <  2:1743994/150‑1 (MQ=255)
                                                                                       ccAAAACCTGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                              <  1:2043419/150‑1 (MQ=255)
                                                                                       ccAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                              >  1:287436/1‑150 (MQ=255)
                                                                                        cAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCa                                                             >  2:34872/1‑150 (MQ=255)
                                                                                         aaaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCAt                                                            >  2:883813/1‑150 (MQ=255)
                                                                                          aaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATc                                                           >  2:1966511/1‑150 (MQ=255)
                                                                                            aCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCa                                                         <  1:539228/150‑1 (MQ=255)
                                                                                            aCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCa                                                         <  2:1306768/150‑1 (MQ=255)
                                                                                             ccGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAg                                                        <  1:1488999/150‑1 (MQ=255)
                                                                                              cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGc                                                       <  2:1875130/150‑1 (MQ=255)
                                                                                                   gtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGtt                                                  >  1:1122251/1‑150 (MQ=255)
                                                                                                   gtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGtt                                                  <  1:522063/150‑1 (MQ=255)
                                                                                                     tGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTTTCTGAAATAACCCTCTCCGAAGGACATCCGGCTACCCGCAGCCTTGCCAGCCATTCATATTAATACACTTACTACAGCACGTTaa                                                >  2:1237558/1‑150 (MQ=255)
                                                                                                           ttCGAGTCTCTCCGACCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCa                                          <  2:142759/150‑1 (MQ=255)
                                                                                                            tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaa                                         >  1:2016814/1‑150 (MQ=255)
                                                                                                              gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCATCTCCGAAGTAAATCCTTCTACCGGCCTCCTTGCCAGCCCTTCATATTAATACACTTCATACAGCACGTTAATTTTCaaaa                                       >  1:1569881/1‑150 (MQ=255)
                                                                                                               aGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAg                                      <  2:778965/150‑1 (MQ=255)
                                                                                                                 tctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                    <  2:25903/150‑1 (MQ=255)
                                                                                                                 tctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                    <  2:1956194/150‑1 (MQ=255)
                                                                                                                   tctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTACGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAATGATcg                                  >  2:1247774/1‑150 (MQ=255)
                                                                                                                    ctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACCCTTCATCCAGCACGTTAATTTTCAAAAGATcgc                                 >  2:1684850/1‑150 (MQ=255)
                                                                                                                    ctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcgc                                 >  2:273092/1‑150 (MQ=255)
                                                                                                                    ctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcgc                                 >  1:723051/1‑150 (MQ=255)
                                                                                                                     cccGCCCCTGACACAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCCCCATCCTTGCCAGCCATTCATATTAATACACATCATCCAGCACGTTAATTTTCAAAAGATcgcg                                <  2:168248/149‑1 (MQ=255)
                                                                                                                      ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                               >  2:443960/1‑150 (MQ=255)
                                                                                                                      ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATGCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                               >  2:1040358/1‑150 (MQ=255)
                                                                                                                       cGCACCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCAACCTTGCCAGCCATTCATATTAATACACTTCATCCAGCAAGTTAATTTTCAAAAGATCGCGaa                              >  2:1390789/1‑150 (MQ=255)
                                                                                                                         ccccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATc                            <  2:1825812/150‑1 (MQ=255)
                                                                                                                           ccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGAATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCaa                          >  2:613076/1‑150 (MQ=255)
                                                                                                                              gCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGc                       >  2:1195980/1‑150 (MQ=255)
                                                                                                                              gCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTAATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGc                       >  1:138739/1‑150 (MQ=255)
                                                                                                                                 aGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAtt                    >  1:435861/1‑150 (MQ=255)
                                                                                                                                  taaaTCATCCTTAGAGAAAGCTAAGGATTTTTTTTATATGAAATAACCCTCTACGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                   <  1:350241/149‑1 (MQ=255)
                                                                                                                                  gAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCGTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                   >  2:1816284/1‑150 (MQ=255)
                                                                                                                                  aaaaaaatcCTTAGCTTTCGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAttt                   >  1:900309/7‑150 (MQ=255)
                                                                                                                                      tcatcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTAt               >  1:1391390/1‑150 (MQ=255)
                                                                                                                                       catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc              <  2:1563227/150‑1 (MQ=255)
                                                                                                                                        atcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCg             <  2:1453880/150‑1 (MQ=255)
                                                                                                                                         tcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGc            <  2:1550719/150‑1 (MQ=255)
                                                                                                                                           cTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTa          <  2:431993/150‑1 (MQ=255)
                                                                                                                                                 gAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGc    <  2:833675/150‑1 (MQ=255)
                                                                                                                                                 gAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGc    >  2:88024/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGCAG  >  NC_000913/3946903‑3947199

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: