Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,860,740 C→G 11.1% M165I (ATG→ATC ygbI ← DeoR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,860,7400CG11.1% 207.7 / 16.4 63M165I (ATG→ATCygbIDeoR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base G (4/3);  total (34/29)
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.67e-01

ATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAG  >  NC_000913/2860602‑2860886
                                                                                                                                          |                                                                                                                                                  
atatTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTg                                                                                                                                           >  3:902637/1‑148 (MQ=255)
    ttCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATtgatca                                                                                                                                                >  4:352713/1‑139 (MQ=255)
         cgCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAatcat                                                                                                                                                  <  2:745906/132‑1 (MQ=255)
         cgCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAatcat                                                                                                                                                  >  1:745906/1‑132 (MQ=255)
                           tCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCg                                                                                                                <  4:902637/148‑1 (MQ=255)
                                  cggcTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAgg                                                                                                         >  4:1095826/1‑148 (MQ=255)
                                        gcactggcaATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGt                                                                                                   <  3:1095826/148‑1 (MQ=255)
                                                 aTCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATTGATCAGGCt                                                                                                                                            <  3:1040345/98‑1 (MQ=255)
                                                 aTCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATTGATCAGGCt                                                                                                                                            >  4:1040345/1‑98 (MQ=255)
                                                          ttCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTgcg                                                                                                                                         <  3:1135341/92‑1 (MQ=255)
                                                          ttCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTgcg                                                                                                                                         >  4:1135341/1‑92 (MQ=255)
                                                                gTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGc                                                                              >  3:285594/1‑145 (MQ=255)
                                                                gTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGc                                                                              <  4:285594/145‑1 (MQ=255)
                                                                      ttATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATTGATCAGGc                                                                                                                                             <  4:986004/76‑1 (MQ=255)
                                                                      ttATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATGATTGATCAGGc                                                                                                                                             >  3:986004/1‑76 (MQ=255)
                                                                      ttATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGc                                                                                                                          >  1:278232/1‑95 (MQ=255)
                                                                      ttATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGc                                                                                                                          <  2:278232/95‑1 (MQ=255)
                                                                         tCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGatcaatc                                                                                                                                                    <  1:74596/66‑1 (MQ=255)
                                                                         tCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGatcaatc                                                                                                                                                    >  2:74596/1‑66 (MQ=255)
                                                                            tCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGacaca                                                                                       >  4:655090/1‑124 (MQ=255)
                                                                            tCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGacaca                                                                                       <  3:655090/124‑1 (MQ=255)
                                                                                      tAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCGTGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATtgt                                                              <  4:570311/139‑1 (MQ=255)
                                                                                      tAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATtgt                                                              >  3:570311/1‑139 (MQ=255)
                                                                                          aaTCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCaa                                                          <  1:532149/139‑1 (MQ=255)
                                                                                          aaTCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCaa                                                          >  2:532149/1‑139 (MQ=255)
                                                                                               cccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTc                                                                                                                    <  2:258624/76‑1 (MQ=255)
                                                                                               cccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTc                                                                                                                    >  1:258624/1‑76 (MQ=255)
                                                                                               cccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACgg                                                                                                    <  1:405630/92‑1 (MQ=255)
                                                                                               cccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACgg                                                                                                    >  2:405630/1‑92 (MQ=255)
                                                                                               cccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAAtt                                                        >  3:562489/1‑136 (MQ=255)
                                                                                               cccGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAAtt                                                        <  4:562489/136‑1 (MQ=255)
                                                                                                 cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaataat                                                                 <  1:201959/125‑1 (MQ=255)
                                                                                                 cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaataat                                                                 >  4:322084/1‑125 (MQ=255)
                                                                                                 cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaataat                                                                 <  3:322084/125‑1 (MQ=255)
                                                                                                 cGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaataat                                                                 >  2:201959/1‑125 (MQ=255)
                                                                                                   cacaCTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAag                                        >  2:989874/1‑148 (MQ=255)
                                                                                                       cTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAagt                                       >  2:994803/1‑145 (MQ=255)
                                                                                                       cTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAagt                                       <  1:994803/145‑1 (MQ=255)
                                                                                                         ccATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTc                                                                                               >  1:795088/1‑87 (MQ=255)
                                                                                                         ccATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTc                                                                                               <  2:795088/87‑1 (MQ=255)
                                                                                                         ccATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTcc                                  >  4:804417/1‑148 (MQ=255)
                                                                                                               cGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATa                            <  3:804417/148‑1 (MQ=255)
                                                                                                                 aTGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACtgttg                                                 <  2:1024208/125‑1 (MQ=255)
                                                                                                                 aTGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACtgttg                                                 >  1:1024208/1‑125 (MQ=255)
                                                                                                                      gAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATa                            <  2:638325/141‑1 (MQ=255)
                                                                                                                      gAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATa                            >  1:638325/1‑141 (MQ=255)
                                                                                                                      gAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAAt                     >  4:402666/1‑148 (MQ=255)
                                                                                                                      gAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAAt                     >  4:124756/1‑148 (MQ=255)
                                                                                                                      gAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAAt                     >  1:929784/1‑148 (MQ=255)
                                                                                                                       aaaTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAgtg                                                                         >  1:371578/1‑95 (MQ=255)
                                                                                                                       aaaTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAgtg                                                                         <  2:371578/95‑1 (MQ=255)
                                                                                                                       aaaTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATtg                                                               <  1:485744/105‑1 (MQ=255)
                                                                                                                       aaaTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATtg                                                               >  2:485744/1‑105 (MQ=255)
                                                                                                                         aTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGt                  <  3:124756/148‑1 (MQ=255)
                                                                                                                         aTAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGt                  <  3:402666/148‑1 (MQ=255)
                                                                                                                           aaaaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaataat                                                                 >  2:642278/1‑99 (MQ=255)
                                                                                                                           aaaaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGaataat                                                                 <  1:642278/99‑1 (MQ=255)
                                                                                                                           aaaaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGttt                >  2:1056149/1‑148 (MQ=255)
                                                                                                                             aaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtg               >  2:860644/1‑147 (MQ=255)
                                                                                                                             aaGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtg               <  1:860644/147‑1 (MQ=255)
                                                                                                                              aGCCTGATGATTGATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgtg             >  4:1050459/1‑148 (MQ=255)
                                                                                                                              aGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTtgtg             <  1:1056149/148‑1 (MQ=255)
                                                                                                                               gCCTGATGATTGATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGa                        >  2:20655/1‑136 (MQ=255)
                                                                                                                               gCCTGATGATTGATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGa                        <  1:20655/136‑1 (MQ=255)
                                                                                                                               gCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAAt                     <  4:394325/139‑1 (MQ=255)
                                                                                                                               gCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAAt                     >  3:394325/1‑139 (MQ=255)
                                                                                                                                   gATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcac        >  4:953408/1‑148 (MQ=255)
                                                                                                                                   gATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACcac        >  2:118307/1‑148 (MQ=255)
                                                                                                                                    atcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAgtgt                                                                        >  4:229287/1‑83 (MQ=255)
                                                                                                                                    atcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAgtgt                                                                        <  3:229287/83‑1 (MQ=255)
                                                                                                                                     tcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCg                   <  3:458068/135‑1 (MQ=255)
                                                                                                                                     tcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCg                   >  4:458068/1‑135 (MQ=255)
                                                                                                                                     tcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGt                  >  1:715207/1‑136 (MQ=255)
                                                                                                                                     tcaatcaTCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGt                  <  2:715207/136‑1 (MQ=255)
                                                                                                                                        atcatcAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtca   >  4:658446/1‑148 (MQ=255)
                                                                                                                                         tcatcaGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCTGCCACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcag  >  4:127525/1‑148 (MQ=255)
                                                                                                                                         tcatcaGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcag  >  3:996282/1‑148 (MQ=255)
                                                                                                                                         tcatcaGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcag  <  1:118307/148‑1 (MQ=255)
                                                                                                                                         tcatcaGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAcagtcag  <  3:953408/148‑1 (MQ=255)
                                                                                                                                         tcatcaGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCAaagtcag  <  2:929784/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
ATATTTCGTCGCATCACAGACCAAAACTCGCTGGCGGCTGGCACTGGCAATCGCCCGTTTCACCGTGACTTTATCTTCCGCTGGCGTAGAAATCCCCCGCACACTCCATGACGATGCAGAAATAAAAGCCTGATCAATCATCAGGCTGCGCAGCATGGTCGCAGCGGCTTCCCCGACACAGGAACGGTTTTCCCGACACACTGCACCGCCAGTGTGAATAATTGTGCAATTACTGTTGTCGAGCAAGTAGTCCGCAATAACGAAATCGTTTGTGACCACAGTCAG  >  NC_000913/2860602‑2860886

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