Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,003 T→C D179G (GAC→GGC)  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,0030TC99.0% 329.2 / NA 98D179G (GAC→GGC) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (43/53);  minor base A (1/0);  total (45/53)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGTCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCT  >  NC_000913/3628857‑3629150
                                                                                                                                                  |                                                                                                                                                   
ccAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCaa                                                                                                                                                  >  4:170353/1‑150 (MQ=255)
ccAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCaa                                                                                                                                                  >  2:765087/1‑150 (MQ=255)
    cACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAAccccc                                                                                                                                             >  4:364338/1‑151 (MQ=255)
       ttCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCggt                                                                                                                                          >  3:239895/1‑149 (MQ=255)
        tctcCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCgtt                                                                                                                                          <  3:406297/150‑1 (MQ=255)
          tcCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCgttgtt                                                                                                                                       >  2:94183/1‑151 (MQ=255)
          tcCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAAACCCCgttgtt                                                                                                                                       >  3:96866/1‑151 (MQ=255)
           ccGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCgttgtt                                                                                                                                       >  2:1289352/1‑150 (MQ=255)
           ccGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCgttgt                                                                                                                                        <  3:275147/149‑1 (MQ=255)
           ccGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCgttg                                                                                                                                         <  2:732805/148‑1 (MQ=255)
              gCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCt                                                                                                                                   <  3:1250737/151‑1 (MQ=255)
              gCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCt                                                                                                                                   <  2:1080788/151‑1 (MQ=255)
                  tCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCtcatc                                                                                                                               >  2:1039057/1‑151 (MQ=255)
                  tCATCGCTAACAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGGGACCAGCACGCTCATATTAATCTGCCGCTGGCGAATAATGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGAGGGCCAACCCCCGGTGTTGGCtc                                                                                                                                  >  3:909548/1‑148 (MQ=255)
                   cATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCtcatct                                                                                                                              <  3:502838/151‑2 (MQ=255)
                     tCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCtcatca                                                                                                                              >  1:326026/1‑149 (MQ=255)
                      cGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAg                                                                                                                            <  1:169338/150‑1 (MQ=255)
                         tACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTTGCTCATCAAGGa                                                                                                                          >  3:1356555/1‑149 (MQ=255)
                            cagccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAAc                                                                                                                     <  3:1103263/151‑1 (MQ=255)
                             agccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAAc                                                                                                                     >  4:281933/1‑150 (MQ=255)
                             agccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACa                                                                                                                    <  2:1299360/151‑1 (MQ=255)
                              gccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAg                                                                                                                   <  1:475517/151‑1 (MQ=255)
                              gccagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCCCATACTGCTCTGCCGCTGGCGAATACTGTCGATTAAATCCCAGAACTGGGAGCGGGGGAGCGGGGCAACCCCCGTTGTTGGCTCCTCCACGGTCAACa                                                                                                                    >  2:366686/1‑150 (MQ=255)
                                cagTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGtt                                                                                                                 <  2:421849/151‑1 (MQ=255)
                                 agTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTc                                                                                                                <  4:142620/151‑1 (MQ=255)
                                  gTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTcc                                                                                                               <  3:827054/151‑1 (MQ=255)
                                    cGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCgg                                                                                                             >  2:9695/1‑151 (MQ=255)
                                    cGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCg                                                                                                              <  2:125454/150‑1 (MQ=255)
                                     gAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCggg                                                                                                            >  3:792864/1‑151 (MQ=255)
                                      aaGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGt                                                                                                           <  2:1046177/151‑1 (MQ=255)
                                       aGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTc                                                                                                          <  2:455626/151‑1 (MQ=255)
                                       aGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTc                                                                                                          <  3:1085378/151‑1 (MQ=255)
                                       aGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTc                                                                                                          <  2:266621/151‑1 (MQ=255)
                                       aGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTc                                                                                                          <  2:947922/151‑1 (MQ=255)
                                        gCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCg                                                                                                         <  3:153430/151‑1 (MQ=255)
                                        gCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCCGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGGGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGATTCGGGTc                                                                                                          >  3:996065/1‑150 (MQ=255)
                                          gTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGGCGTg                                                                                                       >  4:332660/1‑151 (MQ=255)
                                            tCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTgg                                                                                                      <  3:379494/150‑1 (MQ=255)
                                               gCCTCTTCCATATAGGCGGTGGGGACAAGCACACACATATTGTTCTGCCGCTTGCGAATACTGTAGATCAGAACCCAAAAATGGGGGGGGGAGAGCGGGCCAACCCCCGTTGGTGGGTCATCAAGGGACAAAAGTTTCGGGGCGGGGAtt                                                                                                   >  3:1106336/1‑150 (MQ=255)
                                               gCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTa                                                                                                  <  2:1274804/151‑1 (MQ=255)
                                                ccTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTaa                                                                                                 <  3:603308/151‑1 (MQ=255)
                                                 ctctTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGGTGTTTGCTCATCAAGGATCAACAGTTCCGGGTCGTGGAtt                                                                                                   >  2:1100400/1‑148 (MQ=255)
                                                     tCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAAcgcg                                                                                             <  2:539271/150‑1 (MQ=255)
                                                     tCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAAcgc                                                                                              <  2:417915/149‑1 (MQ=255)
                                                        atatAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCgcagc                                                                                         <  1:170355/151‑1 (MQ=255)
                                                          ataGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCgcagca                                                                                        <  2:1299487/150‑1 (MQ=255)
                                                            aGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGcaca                                                                                      <  3:1066665/150‑1 (MQ=255)
                                                            aGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGcaca                                                                                      <  1:108856/150‑1 (MQ=255)
                                                              gCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGTGGGGGAGAGCGGGACAACCCCCGTTGGCGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACTCGCAGCACAGcc                                                                                   >  3:1143569/1‑151 (MQ=255)
                                                              gCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGcc                                                                                   >  3:1373399/1‑151 (MQ=255)
                                                                 gTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCaa                                                                                <  4:98320/151‑1 (MQ=255)
                                                                 gTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCa                                                                                 >  1:180434/1‑150 (MQ=255)
                                                                   ggCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAg                                                                               <  4:171690/150‑1 (MQ=255)
                                                                    gCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGtt                                                                             <  2:719496/151‑1 (MQ=255)
                                                                     cGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGttt                                                                            >  3:1027966/1‑151 (MQ=255)
                                                                      gACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGtttt                                                                           >  4:399264/1‑151 (MQ=255)
                                                                         cAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGc                                                                        >  3:1181105/1‑151 (MQ=255)
                                                                             aCGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTc                                                                     >  2:147536/1‑150 (MQ=255)
                                                                             aCGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTc                                                                     >  2:596485/1‑150 (MQ=255)
                                                                             aCGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTc                                                                     >  3:166803/1‑150 (MQ=255)
                                                                             aCGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCa                                                                    <  2:96851/151‑1 (MQ=255)
                                                                               gCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATc                                                                  >  2:1139674/1‑151 (MQ=255)
                                                                                 tCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATccc                                                                >  4:126211/1‑151 (MQ=255)
                                                                                   atatTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCccgc                                                              <  1:383689/151‑1 (MQ=255)
                                                                                   atatTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCccgc                                                              >  3:914887/1‑151 (MQ=255)
                                                                                       tGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCgg                                                           <  2:290644/150‑1 (MQ=255)
                                                                                          tctGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGgag                                                         >  2:467265/1‑149 (MQ=255)
                                                                                           ctGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGt                                                      >  3:345131/1‑151 (MQ=255)
                                                                                             gCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGttt                                                    <  1:281936/151‑1 (MQ=255)
                                                                                              ccGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTc                                                   >  4:139281/1‑151 (MQ=255)
                                                                                               cGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTcc                                                  >  4:248651/1‑151 (MQ=255)
                                                                                                   ggCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTg                                               <  3:155580/150‑1 (MQ=255)
                                                                                                   ggCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGc                                              >  3:197480/1‑151 (MQ=255)
                                                                                                      gAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCgg                                           <  3:871358/151‑1 (MQ=255)
                                                                                                       aaTACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCgg                                           >  2:880781/1‑150 (MQ=255)
                                                                                                         tACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACg                                        >  1:95101/1‑151 (MQ=255)
                                                                                                          aCTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGCGTTTCCCTgccggccg                                        >  3:1313071/1‑150 (MQ=255)
                                                                                                          aCTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACg                                        <  3:1139729/150‑1 (MQ=255)
                                                                                                            tGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATc                                     >  3:693046/1‑151 (MQ=255)
                                                                                                             gTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGAt                                      <  3:94197/149‑1 (MQ=255)
                                                                                                               cGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATcgcg                                  <  3:834426/151‑1 (MQ=255)
                                                                                                               cGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATcgc                                   <  4:283319/150‑1 (MQ=255)
                                                                                                                   cAGATCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAAc                              >  3:622040/1‑151 (MQ=255)
                                                                                                                      aTCCCAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGt                           >  4:13961/1‑151 (MQ=255)
                                                                                                                        cccAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGa                                          >  1:76496/1‑134 (MQ=255)
                                                                                                                        cccAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGa                                          <  4:76496/134‑1 (MQ=255)
                                                                                                                         ccAGAACTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCt                        <  3:1039104/151‑1 (MQ=255)
                                                                                                                              aCTGGGAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAAccc                   >  2:715004/1‑151 (MQ=255)
                                                                                                                                 gggAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGt                >  1:468620/1‑151 (MQ=255)
                                                                                                                                 gggAGCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGt                <  3:834646/151‑1 (MQ=255)
                                                                                                                                     gCGGGAGAGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCt             <  3:1289420/150‑1 (MQ=255)
                                                                                                                                         gagaGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTca        <  2:693021/151‑1 (MQ=255)
                                                                                                                                         gagaGCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTca        <  2:723143/151‑1 (MQ=255)
                                                                                                                                          agagCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTcag       <  2:166778/151‑1 (MQ=255)
                                                                                                                                           gagCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTcagc      <  3:1197525/151‑1 (MQ=255)
                                                                                                                                           gagCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTcag       <  3:765122/150‑1 (MQ=255)
                                                                                                                                            agCGGGCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTcagca     <  2:783910/151‑1 (MQ=255)
                                                                                                                                               gggCCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCt  <  1:18221/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCAGCACTTCTCCGGCATTCATCGCTACCAGCCAGTCGAAGCGTTCGGCCTCTTCCATATAGGCGGTGGCGACCAGCACGCTCATATTGCTCTGCCGCTGGCGAATACTGTCGATCAGATCCCAGAACTGGGAGCGGGAGAGCGGGTCAACCCCCGTTGTTGGCTCATCAAGGATCAACAGTTCCGGGTCGTGGATTAACGCGCAGCACAGCCCAAGTTTTTGCTTCATCCCGCCGGAGAGTTTCCCTGCCGGACGATCGCGAAACGGTGCTAACCCGGTGCTGGTCAGCAGCT  >  NC_000913/3628857‑3629150

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑