Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,581 G→T 9.8% intergenic (‑424/+19) rrsG ← / ← clpB 16S ribosomal RNA/ClpB chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,5810GT 9.8% 263.4 / 12.2 81intergenic (‑424/+19)rrsG/clpB16S ribosomal RNA/ClpB chaperone
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base T (5/3);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGA  >  NC_000913/2731437‑2731729
                                                                                                                                                |                                                                                                                                                    
ttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAggg                                                                                                                                                 >  1:225728/1‑150 (MQ=255)
ttCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAggg                                                                                                                                                 <  1:1553935/150‑1 (MQ=255)
    gtgtCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCg                                                                                                                                             >  2:1932056/1‑150 (MQ=255)
      gtCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtt                                                                                                                                           <  1:1504641/150‑1 (MQ=255)
      gtCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtt                                                                                                                                           <  2:1800656/150‑1 (MQ=255)
        cAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCTTCGGGGCTCGtttt                                                                                                                                         <  2:1346838/150‑1 (MQ=255)
        cAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGtttt                                                                                                                                         >  2:637408/1‑150 (MQ=255)
          gtggtgGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATTTAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTAt                                                                                                                                       <  2:453925/150‑1 (MQ=255)
          gtggtgGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTAt                                                                                                                                       >  2:1105033/1‑150 (MQ=255)
            ggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCa                                                                                                                                     <  2:1363564/150‑1 (MQ=255)
            ggtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCa                                                                                                                                     <  1:446196/150‑1 (MQ=255)
             gtggCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGtttttgtct                                                                                                                                    <  2:225728/150‑6 (MQ=255)
                gcgcATTATAGGGAGTTATTCCGGCCTGACAAGAGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTa                                                                                                                                 >  1:279452/1‑150 (MQ=255)
                 cgcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTAc                                                                                                                                >  2:1663054/1‑150 (MQ=255)
                  gcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                               <  1:1567066/150‑1 (MQ=255)
                  gcATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACt                                                                                                                               <  1:1580048/150‑1 (MQ=255)
                   cATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGCCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTg                                                                                                                              <  1:148726/150‑1 (MQ=255)
                          gggAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCg                                                                                                                       >  1:1254937/1‑150 (MQ=255)
                              gTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACaa                                                                                                                   <  1:1105033/150‑1 (MQ=255)
                                 aTTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATcc                                                                                                                <  2:1016915/150‑1 (MQ=255)
                                     cGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTc                                                                                                            >  1:655914/1‑150 (MQ=255)
                                      ggCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCt                                                                                                           <  1:168644/150‑1 (MQ=255)
                                      ggCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCt                                                                                                           <  1:1871370/150‑1 (MQ=255)
                                      ggCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCt                                                                                                           <  2:1615254/150‑1 (MQ=255)
                                        ccTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTc                                                                                                         <  1:1600377/150‑1 (MQ=255)
                                         cTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCa                                                                                                        >  2:1478170/1‑150 (MQ=255)
                                           gACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCAtt                                                                                                      >  1:1274092/1‑150 (MQ=255)
                                               aGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACt                                                                                                  <  1:166802/150‑1 (MQ=255)
                                                gCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACtt                                                                                                 >  2:1756020/1‑150 (MQ=255)
                                                  gAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTcc                                                                                               <  2:279452/150‑1 (MQ=255)
                                                   aaaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCa                                                                                              >  1:844517/1‑150 (MQ=255)
                                                    aaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCtta                                                                                                         <  1:1064517/138‑2 (MQ=255)
                                                    aaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTc                                                                                                         >  2:1064517/1‑138 (MQ=255)
                                                    aaaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAg                                                                                             >  2:1362319/1‑150 (MQ=255)
                                                     aaaTATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAgg                                                                                            >  2:746509/1‑150 (MQ=255)
                                                       atatAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGGCCCGAAGGGCGCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTGATTAACTTCCAGGCg                                                                                          <  2:1415242/150‑1 (MQ=255)
                                                        tataAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGa                                                                                         >  2:451165/1‑150 (MQ=255)
                                                          taAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACGTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAAt                                                                                       >  2:44947/1‑150 (MQ=255)
                                                            aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGAc                                                                                                                           <  2:422121/112‑1 (MQ=255)
                                                            aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGAc                                                                                                                           >  1:422121/1‑112 (MQ=255)
                                                            aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTcc                                                                                               <  1:871549/140‑1 (MQ=255)
                                                            aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTcc                                                                                               >  2:871549/1‑140 (MQ=255)
                                                            aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCa                                                                                     <  2:1274092/150‑1 (MQ=255)
                                                            aaaaCTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCa                                                                                     >  2:661194/1‑150 (MQ=255)
                                                                     tCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCg                                                                            >  2:711844/1‑150 (MQ=255)
                                                                       gTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGa                                                                          <  1:711844/150‑1 (MQ=255)
                                                                         tCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAAc                                                                        <  2:1574105/150‑1 (MQ=255)
                                                                            cTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCaa                                                                     >  1:1273700/1‑150 (MQ=255)
                                                                            cTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCaa                                                                     >  1:1311817/1‑150 (MQ=255)
                                                                             tCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAt                                                                    <  2:885397/150‑1 (MQ=255)
                                                                              cACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAtt                                                                   >  2:601215/1‑150 (MQ=255)
                                                                              cACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAAtt                                                                   >  1:1735695/1‑150 (MQ=255)
                                                                                  tttCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCAcc                                                               <  1:1985979/150‑1 (MQ=255)
                                                                                   ttCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCa                                                              >  2:280877/1‑150 (MQ=255)
                                                                                      aGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGAc                                                           <  1:686430/150‑1 (MQ=255)
                                                                                         cAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGt                                                        >  1:1657198/1‑150 (MQ=255)
                                                                                         cAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGt                                                        >  1:862581/1‑150 (MQ=255)
                                                                                                           gTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCggg                                      >  1:1604686/1‑150 (MQ=255)
                                                                                                                ttCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTc                                 <  1:1276344/150‑1 (MQ=255)
                                                                                                                 tCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTAccc                                                                             <  1:907302/105‑1 (MQ=255)
                                                                                                                 tCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCAGGTCTTCATTAACTTCCAGGCGAATCACTTTAccc                                                                             >  2:907302/1‑105 (MQ=255)
                                                                                                                   cGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGAt                              >  1:371264/1‑150 (MQ=255)
                                                                                                                    gTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCTGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATc                             <  1:1182974/150‑1 (MQ=255)
                                                                                                                     tCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATct                            <  1:1478170/150‑1 (MQ=255)
                                                                                                                       tAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgc                          >  1:975251/1‑150 (MQ=255)
                                                                                                                          cTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgc                       >  2:1798439/1‑150 (MQ=255)
                                                                                                                            tataGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATctgctgctg                     >  2:404466/1‑150 (MQ=255)
                                                                                                                              taGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGaa                   >  1:1911799/1‑150 (MQ=255)
                                                                                                                              aatgataaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGaa                   >  2:1095623/7‑150 (MQ=255)
                                                                                                                               atgataaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAAt                  >  1:1384809/6‑150 (MQ=255)
                                                                                                                                 aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtttt                                  <  2:133070/132‑1 (MQ=255)
                                                                                                                                 aCAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGtttt                                  >  1:133070/1‑132 (MQ=255)
                                                                                                                                  cAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGc               >  1:442880/1‑150 (MQ=255)
                                                                                                                                   taaaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCa              >  1:14516/2‑150 (MQ=255)
                                                                                                                                    aaaaCGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCAc             <  2:1311817/150‑1 (MQ=255)
                                                                                                                                     aaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACg            >  2:654501/1‑150 (MQ=255)
                                                                                                                                     aaaCGAGCCCTTCGGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACg            >  1:1307831/1‑150 (MQ=255)
                                                                                                                                           gCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCag      <  2:1254937/150‑1 (MQ=255)
                                                                                                                                           gCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCag      >  2:192137/1‑150 (MQ=255)
                                                                                                                                              ccGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAgg   >  1:173665/1‑150 (MQ=255)
                                                                                                                                               cGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGa  <  2:1657198/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCCGTGTCAGTGGTGGCGCATTATAGGGAGTTATTCCGGCCTGACAAGCGAAAAATATAAAAACTTTATCGTTCGCTCACTTTTCAGGCAAAACATCTTAAATATAGTCTTTTCCGTCTAACTTATAGACAAAAACGAGCCCCGAAGGGCTCGTTTTATCATTTACTGGACGGCGACAATCCGGTCTTCATTAACTTCCAGGCGAATCACTTTACCCGGAACCAATTCACCAGACAGTATTTGCTGTGCCAGCGGGTTTTCGATCTGCTGCTGAATTGCACGTTTCAGAGGA  >  NC_000913/2731437‑2731729

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