Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,325,677 A→G 5.3% I420T (ATT→ACT)  ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,325,6770AG 5.3% 229.4 / 2.5 75I420T (ATT→ACT) ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

GCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAA  >  NC_000913/3325544‑3325825
                                                                                                                                     |                                                                                                                                                    
gCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCtt                                                                                                                                       >  8:281350/1‑149 (MQ=255)
 cTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTc                                                                                                                                      <  8:205968/149‑1 (MQ=255)
  tttGCAGTTTCTGTCGGCTGGGGCTGTTTGCGTCGCCCTCAGTCTAGAAGAAAGTCACACCCAGCTCACGTCGGCGTTGGATAATCGTCAGTTTGTGCACCGGTTCGTGTTCAGTCACCAGGCGACCTATAATCGCATGACCCGCTTCg                                                                                                                                     <  8:106164/149‑1 (MQ=255)
  tttCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCg                                                                                                                                     >  8:31291/1‑149 (MQ=255)
   ttCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGt                                                                                                                                    >  6:136025/1‑149 (MQ=255)
             tGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGt                                                                                                                          <  8:121093/149‑1 (MQ=255)
               aCGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGt                                                                                                                          <  8:329086/147‑1 (MQ=255)
                cGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGa                                                                                                                       <  8:311846/149‑1 (MQ=255)
                    tGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGCCCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGGGGTAAGCCGTCGAttct                                                                                                                   >  5:126640/1‑149 (MQ=255)
                      gcgcTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGAttctt                                                                                                                  <  1:331591/148‑1 (MQ=255)
                        gcTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGAttctttct                                                                                                               <  1:284783/149‑1 (MQ=255)
                        gcTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGAttctttct                                                                                                               <  1:87712/149‑1 (MQ=255)
                          tGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTgc                                                                                                             >  4:314608/1‑149 (MQ=255)
                            aTTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTgcgc                                                                                                           >  7:39620/1‑149 (MQ=255)
                                     cccTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGtcat                                                                                                  <  7:281350/149‑1 (MQ=255)
                                      ccTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCCCCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGtcat                                                                                                  >  8:334462/1‑148 (MQ=255)
                                            gCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATg                                                                                           <  7:31291/149‑1 (MQ=255)
                                            gCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATg                                                                                           >  8:81474/1‑149 (MQ=255)
                                            gCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATg                                                                                           >  3:124803/1‑149 (MQ=255)
                                              aagaagAAAGTCACACCCAGCGCACGTCCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAg                                                                                                                              <  6:250577/112‑1 (MQ=255)
                                              aagaagAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAg                                                                                                                              >  5:250577/1‑112 (MQ=255)
                                              aagaagAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCCCCATg                                                                                           >  3:163487/1‑147 (MQ=255)
                                              aagaagAAAGTCACACCCAGCGAACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCTATTCTTTCTGCGCTTCAGTCATCACCATg                                                                                           >  2:123924/1‑147 (MQ=255)
                                               agaagaAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCCCCGGATCGGGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCGTCGTGGTAAGCCGTCGATTCGTTCTGCGCTTCCGTCATCACCATggcg                                                                                        >  4:98614/1‑147 (MQ=255)
                                                 aagaaAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGc                                                                                       >  7:25827/1‑148 (MQ=255)
                                                   gaaAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGAGAATCGCATGACCCGCTTCGTGGTAAGCCGGCGATTCTTTCTGCGCTTCCGTCATCACCCTGGAGCGAc                                                                                    >  4:164806/1‑149 (MQ=255)
                                                    aaaGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCTGTCATCACCATGGAGCGACg                                                                                   >  2:155437/1‑149 (MQ=255)
                                                    aaaGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACg                                                                                   <  8:39620/149‑1 (MQ=255)
                                                          acacCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGc                                                                             <  7:250570/149‑1 (MQ=255)
                                                            acCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGTCGACCGATAATCGCATGAACCGCTTCGTAGTAAGACGTCGATTCTTTCTGCGCTTCCGTCATCACCATGTATCGACGTTCCTCAc                                                                           >  2:100739/1‑149 (MQ=255)
                                                              ccAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCAcc                                                                          >  6:337177/1‑148 (MQ=255)
                                                              ccAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCAcc                                                                          >  2:167976/1‑148 (MQ=255)
                                                              ccAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCAcc                                                                          >  1:95240/1‑148 (MQ=255)
                                                                      cGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGAt                                                                 >  7:161261/1‑149 (MQ=255)
                                                                       gACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTcc                                                                               <  1:150135/134‑1 (MQ=255)
                                                                       gACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTcc                                                                               >  2:150135/1‑134 (MQ=255)
                                                                         ccGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGAtttt                                                              <  6:171050/149‑1 (MQ=255)
                                                                         ccGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGGGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGAtttt                                                              <  6:126640/149‑1 (MQ=255)
                                                                             gTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCt                                                          >  1:345896/1‑149 (MQ=255)
                                                                                ggATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTc                                                       <  4:113937/149‑1 (MQ=255)
                                                                                ggATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTc                                                       <  5:28682/149‑1 (MQ=255)
                                                                                ggATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTc                                                       <  3:314608/149‑1 (MQ=255)
                                                                                      tCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTtc                                                 >  2:12722/1‑149 (MQ=255)
                                                                                         tCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCAc                                                                           <  4:196325/120‑1 (MQ=255)
                                                                                         tCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCAc                                                                           >  3:196325/1‑120 (MQ=255)
                                                                                           aCTTTGTGCACCGGATCGTGTTCCGGCACCAGTCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGaa                                            <  2:345896/149‑1 (MQ=255)
                                                                                             tttGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACt                                          >  3:261374/1‑149 (MQ=255)
                                                                                             tttGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACt                                          >  4:227830/1‑149 (MQ=255)
                                                                                               tgtgCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTc                                         <  7:81474/148‑1 (MQ=255)
                                                                                               tgtgCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTc                                         >  6:29804/1‑148 (MQ=255)
                                                                                                 tgCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCaa                                       <  8:161261/148‑1 (MQ=255)
                                                                                                 tgCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCaa                                       <  1:123924/148‑1 (MQ=255)
                                                                                                 tgCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAAc                                      >  8:20740/1‑149 (MQ=255)
                                                                                                   cACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAAcc                                     <  8:25827/148‑1 (MQ=255)
                                                                                                   cACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCa                                    >  3:20641/1‑149 (MQ=255)
                                                                                                    aCCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCAt                                   >  7:74507/1‑149 (MQ=255)
                                                                                                     ccGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTTGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATc                                  <  1:28576/149‑1 (MQ=255)
                                                                                                     ccGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGGCTTTCGCTTTCTCGAACTCAACCATc                                  <  2:198146/149‑1 (MQ=255)
                                                                                                        gATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGa                                >  1:121283/1‑148 (MQ=255)
                                                                                                         aTCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCCCCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGaca                              >  1:32588/1‑149 (MQ=255)
                                                                                                                tCCGGCACCAGGCGACCGATAGTCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTATGCGCTTCCGTCATCACCATGTAGCGACGATCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACAATCGACAAAAagcg                       >  6:195148/1‑149 (MQ=255)
                                                                                                                tCCGGCACCAGGCGACCGATAGTCGCATGACCCGCTTCGTGGAAAGCCGTCGAGTCTTTCTGCGCTTACCTCATCACCATGGAGCGACGTTCCGCACCCCTCATGATTTTGTCTTTCGCGTTCTCGAACTCAACCATCTCCACAAcgcg                       >  6:200023/1‑149 (MQ=255)
                                                                                                                    gCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATTGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttg                   <  7:334462/149‑1 (MQ=255)
                                                                                                                      aCCCGGCGACCGATAGTCGCATGACCCGCTTCGTGGTAAGCCGTCTATTCTTTATGAGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCTAACTCAACCATCGACACAACGCGTttgtt                 <  5:195148/149‑1 (MQ=255)
                                                                                                                      aCCAGGCGACCGGTAGTCGCATGTCCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGGGACGTTCCGCACCCATGATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttgtt                 <  5:200023/149‑1 (MQ=255)
                                                                                                                      aCCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCTGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttgtt                 <  1:155437/149‑1 (MQ=255)
                                                                                                                      aCCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttgtt                 >  8:169523/1‑149 (MQ=255)
                                                                                                                      aCCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGGCGATTCTTTCTGCGCGTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttgtt                 <  8:168707/149‑1 (MQ=255)
                                                                                                                        cAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttgtt                 <  4:81382/147‑1 (MQ=255)
                                                                                                                        cAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTttgtt                 >  3:81382/1‑147 (MQ=255)
                                                                                                                        cAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGc               >  3:260328/1‑149 (MQ=255)
                                                                                                                        cAGGCGACCGATAATCGCATGACCCGCTTCGTGCTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGc               >  4:63430/1‑149 (MQ=255)
                                                                                                                           gCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCAc            <  5:337177/149‑1 (MQ=255)
                                                                                                                             gACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGa          <  1:167976/149‑1 (MQ=255)
                                                                                                                                cGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGag         >  2:297358/1‑147 (MQ=255)
                                                                                                                                    aaTCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCATTTGTTGCCACGAGCAGCg    >  8:228690/1‑148 (MQ=255)
                                                                                                                                     aTCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGaa  <  4:124803/149‑1 (MQ=255)
                                                                                                                                     aTCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGaa  <  3:227830/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
GCTTTCCAGTTTCTGACGGCTGGCGCTGATTGCGTCGCCCTCAGGCAAGAAGAAAGTCACACCCAGCGCACGACCGCGTGGGATAATCGTCACTTTGTGCACCGGATCGTGTTCCGGCACCAGGCGACCGATAATCGCATGACCCGCTTCGTGGTAAGCCGTCGATTCTTTCTGCGCTTCCGTCATCACCATGGAGCGACGTTCCGCACCCATCATGATTTTGTCTTTCGCTTTCTCGAACTCAACCATCGACACAACGCGTTTGTTGCCACGAGCAGCGAA  >  NC_000913/3325544‑3325825

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