Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,629,498 T→C 9.9% Q14R (CAA→CGA)  rbbA ← ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,4980TC 9.9% 297.0 / 19.3 102Q14R (CAA→CGA) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base T (50/41);  new base C (6/4);  total (57/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 106 total aligned reads displayed.

GACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGT  >  NC_000913/3629350‑3629646
                                                                                                                                                     |                                                                                                                                                    
gACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGtt                                                                                                                                                      <  1:127658/149‑1 (MQ=255)
         aCCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAg                                                                                                                                             >  1:643782/1‑149 (MQ=255)
           cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGc                                                                                                                                           >  4:508833/1‑149 (MQ=255)
           cGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGc                                                                                                                                           >  4:784829/1‑149 (MQ=255)
             gAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgg                                                                                                                                         <  1:557755/149‑1 (MQ=255)
             aaaaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCTCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgg                                                                                                                                         >  2:540425/2‑149 (MQ=255)
               aaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTCGCGCAACTGGcgggc                                                                                                                                       <  2:292749/149‑2 (MQ=255)
               aaTCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCTGTCGCGCAACTGGcgggc                                                                                                                                       <  2:810088/149‑2 (MQ=255)
                 tCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacg                                                                                                                                     >  3:480551/1‑149 (MQ=255)
                  cAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                    <  1:917229/149‑1 (MQ=255)
                   aaCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacggg                                                                                                                                   >  1:297081/1‑149 (MQ=255)
                    aCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                                                  <  2:511348/149‑1 (MQ=255)
                    aCGACAACAAGCTCGACTTCCCGACACCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGTTCCCGTCCCGCCTGTCGCGCAACTGGcgggcgtga                                                                                                                                   <  4:3776/149‑6 (MQ=37)
                      gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                                                                 >  2:845262/1‑148 (MQ=255)
                      gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                                                                 <  1:845262/148‑1 (MQ=255)
                      gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAAc                                                                                                                                >  1:400596/1‑149 (MQ=255)
                         aacaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAg                                                                                                                             >  3:489617/1‑149 (MQ=255)
                          acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGt                                                                                                                            >  3:393027/1‑149 (MQ=255)
                          acaaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGt                                                                                                                            >  1:886300/1‑149 (MQ=255)
                           caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTgc                                                                                                                                                    >  2:321459/1‑124 (MQ=255)
                           caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTgc                                                                                                                                                    <  1:321459/124‑1 (MQ=255)
                           caaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGtt                                                                                                                           >  4:333095/1‑149 (MQ=255)
                             aGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGAGCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                                                                         >  4:103891/1‑149 (MQ=255)
                              gCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCa                                                                                                                        >  4:106071/1‑149 (MQ=255)
                               cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTACATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAg                                                                                                                       >  1:97489/1‑149 (MQ=255)
                                tCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                                                                      >  3:451295/1‑149 (MQ=255)
                                 cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCTGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                                                                     <  2:809499/149‑1 (MQ=255)
                                   aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGc                                                                                                                   <  4:480551/149‑1 (MQ=255)
                                   aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGc                                                                                                                   <  4:668824/149‑1 (MQ=255)
                                     ttCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                                 <  2:252279/149‑1 (MQ=255)
                                     ttCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGt                                                                                                                 <  4:385413/149‑1 (MQ=255)
                                       cccGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                                                                               >  2:380265/1‑149 (MQ=255)
                                        ccGACGCCGTCCGGGCCAATCAGCCCGACAATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCAt                                                                                                              >  4:346146/1‑149 (MQ=255)
                                          gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTg                                                                                                            >  4:127056/1‑149 (MQ=255)
                                             gCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                    >  3:660292/1‑122 (MQ=255)
                                             gCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                                                    <  4:660292/122‑1 (MQ=255)
                                                gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                                                                      >  2:508849/1‑149 (MQ=255)
                                                gTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCaa                                                                                                      <  1:384242/149‑1 (MQ=255)
                                                       ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAt                                                                                                                     >  3:393082/1‑127 (MQ=255)
                                                       ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacc                                                                                               >  4:541952/1‑149 (MQ=255)
                                                        cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                                                                              <  4:802605/149‑1 (MQ=255)
                                                        cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTcacca                                                                                              >  1:697866/1‑149 (MQ=255)
                                                          aTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGtt                                                                                                                                                      >  4:831519/1‑91 (MQ=255)
                                                          aTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGtt                                                                                                                                                      <  3:831519/91‑1 (MQ=255)
                                                          aTCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTGGCGCGACAGGCGGCACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACg                                                                                            >  1:721131/1‑149 (MQ=255)
                                                             aGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                                                                         <  3:508833/149‑1 (MQ=255)
                                                              gCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgt                                                                                        >  4:407533/1‑149 (MQ=255)
                                                               cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                       <  2:676127/149‑1 (MQ=255)
                                                               cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                       >  3:575629/1‑149 (MQ=255)
                                                               cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtc                                                                                       <  2:524881/149‑1 (MQ=255)
                                                                ccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcg                                                                                      >  3:925642/1‑149 (MQ=255)
                                                                 cGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                                                     <  3:905433/149‑1 (MQ=255)
                                                                  gACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgtc                                                                                    <  1:521496/149‑1 (MQ=255)
                                                                  gACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGAGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgtc                                                                                    >  3:762384/1‑149 (MQ=255)
                                                                   aCCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCa                                                                                   <  2:646893/149‑1 (MQ=255)
                                                                             gggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                         <  3:117398/149‑1 (MQ=255)
                                                                               gCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGt                                                                       >  2:767175/1‑149 (MQ=255)
                                                                                   gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                                                                   <  4:575629/149‑1 (MQ=255)
                                                                                   gAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                                                                   <  1:707289/149‑1 (MQ=255)
                                                                                    aaTATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttcttc                                                                  <  3:333095/149‑1 (MQ=255)
                                                                                        tCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcattc                                                              >  3:657972/1‑149 (MQ=255)
                                                                                          gagaGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAc                                                            <  3:794591/149‑1 (MQ=255)
                                                                                            gagTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAccc                                                          >  2:66643/1‑149 (MQ=255)
                                                                                                 atatTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc                                                     <  3:503738/149‑1 (MQ=255)
                                                                                                   atTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCa                                                   >  3:97234/1‑149 (MQ=255)
                                                                                                             gcAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGccc                                         <  3:574292/149‑1 (MQ=255)
                                                                                                               aaCGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCgg                                       >  2:355546/1‑149 (MQ=255)
                                                                                                                aCGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGc                                      >  2:42291/1‑149 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                     >  1:327372/1‑108 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                                                                     <  2:327372/108‑1 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                            <  1:380265/149‑1 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                            >  1:354072/1‑149 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                            <  4:393027/149‑1 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                            <  2:237613/149‑1 (MQ=255)
                                                                                                                          cATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGttt                            >  1:171961/1‑149 (MQ=255)
                                                                                                                           aTAATGCTGGCTCACG‑CCCGCCTGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCCTTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCTGCCAAACCGGtttt                           >  2:759018/1‑149 (MQ=255)
                                                                                                                           aTAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                         <  3:826582/103‑1 (MQ=255)
                                                                                                                           aTAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                                                                         >  4:826582/1‑103 (MQ=255)
                                                                                                                           aTAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                           >  2:557575/1‑149 (MQ=255)
                                                                                                                           aTAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                           >  4:479788/1‑149 (MQ=255)
                                                                                                                           aTAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                           <  2:293019/149‑1 (MQ=255)
                                                                                                                               tGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACa                       >  1:29776/1‑149 (MQ=255)
                                                                                                                                     ttcccgtcccgccTGTCGCGCAACTGGCTGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTc                  >  2:831589/8‑149 (MQ=255)
                                                                                                                                      tCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                >  2:190894/1‑149 (MQ=255)
                                                                                                                                      tCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                >  4:123404/1‑149 (MQ=255)
                                                                                                                                      tCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                <  2:643782/149‑1 (MQ=255)
                                                                                                                                      tCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCa                >  3:633603/1‑149 (MQ=255)
                                                                                                                                      tCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGAGGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGACACGGAAGGTCTTCATTAACCCGCACCCACGCTACGCCCGGCAAACCAGTTTTGACATATTCCa                >  4:752738/1‑149 (MQ=255)
                                                                                                                                        aCG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                           <  4:285043/136‑1 (MQ=255)
                                                                                                                                        aCG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtttt                           >  3:285043/1‑136 (MQ=255)
                                                                                                                                          gtcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAt             >  2:388003/3‑149 (MQ=255)
                                                                                                                                          gtcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAt             >  2:629231/3‑149 (MQ=255)
                                                                                                                                           tcccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtg            >  3:11869/2‑149 (MQ=255)
                                                                                                                                            cccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtt                             >  4:558789/1‑131 (MQ=255)
                                                                                                                                            cccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGtt                             <  3:558789/131‑1 (MQ=255)
                                                                                                                                             ccgccTGTCGCGCAACTGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGAtgct          >  2:799214/1‑149 (MQ=255)
                                                                                                                                                   gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGa    <  4:451295/149‑1 (MQ=255)
                                                                                                                                                   gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGa    <  2:697866/149‑1 (MQ=255)
                                                                                                                                                   gTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGa    <  1:190894/149‑1 (MQ=255)
                                                                                                                                                     tGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGt  <  4:395324/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACG‑CCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAACCGGTTTTGACATATTCCAGATGCTGCTGGAGT  >  NC_000913/3629350‑3629646

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