Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,641,005 G→A S48S (TCG→TCA dtpB → dipeptide and tripeptide permease B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,641,0050GA89.0% 434.6 / 2.7 139S48S (TCG→TCAdtpBdipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (42/80);  minor base C (15/0);  total (58/81)
Fisher's exact test for biased strand distribution p-value = 5.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 139 total aligned reads displayed.

AAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCAT  >  NC_000913/3640858‑3641151
                                                                                                                                                      |                                                                                                                                                    
aaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACaa                                                                                                                                                   <  2:553420/151‑1 (MQ=255)
  atatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAaga                                                                                                                                                 <  4:722126/151‑1 (MQ=255)
   tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTAAAAAGCTTGGGTTTTCACaa                                                                                                                                                   >  3:437416/1‑148 (MQ=255)
    atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAagag                                                                                                                                                <  2:1673402/150‑1 (MQ=255)
    atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAagag                                                                                                                                                <  4:480501/150‑1 (MQ=255)
      gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCa                                                                                                                                              <  4:1382050/150‑1 (MQ=255)
      gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAg                                                                                                                                             <  3:336541/151‑1 (MQ=255)
       aaTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTCTTTAAAAGCTTTGGTTTTCCCAAAAGCAgg                                                                                                                                            >  4:693968/1‑151 (MQ=255)
       aaTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAg                                                                                                                                             <  2:223613/150‑1 (MQ=255)
          acaacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCtt                                                                                                                                        <  2:1285398/151‑1 (MQ=255)
          acaacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCtt                                                                                                                                        <  2:592045/151‑1 (MQ=255)
           caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTAAAAAAGTTGGGTTTTCACAAAAGCAggggtt                                                                                                                                        >  2:84824/1‑147 (MQ=255)
            aacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCTTTTAAAAACTTGGATTTTCCCAAAAGCAGGCtttt                                                                                                                                      >  2:610272/1‑151 (MQ=255)
            aacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCtttt                                                                                                                                      <  4:721089/151‑1 (MQ=255)
              caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAGCTTGGATTCTCACAAGAGCAGGCtttt                                                                                                                                      >  3:126401/1‑149 (MQ=255)
               aacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAACAGCTTGGGTTTTCACAAAAGCAGGCgttttt                                                                                                                                    >  1:14674/1‑148 (MQ=255)
                acacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAACTTGGGTTCTCACAAGAGCAGGCtttttt                                                                                                                                    >  3:196733/1‑147 (MQ=255)
                    ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACttt                                                                                                                              <  4:13558/151‑1 (MQ=255)
                    ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACttt                                                                                                                              <  1:19399/151‑1 (MQ=255)
                     cATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTTAACAACTTGGGTTCTCCCAAAAGCAGGGCTTTTTCTCtttt                                                                                                                             >  2:927652/1‑151 (MQ=255)
                      aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTCTTCGTTTAACAGCCTTGATTCTCCCAAGAGGAGGGTTTTTTTACctttt                                                                                                                            >  3:342012/1‑150 (MQ=255)
                       tGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCCTTTAAAACCTTGGGTTTTCCCAAAAGCAGGCTTTTTTTCCtttt                                                                                                                             >  4:278148/1‑149 (MQ=255)
                        ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTgg                                                                                                                           <  2:1212720/150‑1 (MQ=255)
                         gggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTg                                                                                                                         <  2:725684/151‑1 (MQ=255)
                           gATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAAGTTTGGTTTTCCCAAGAGGAGGCTTTTTTTACtttttggggt                                                                                                                       >  4:365572/1‑145 (MQ=255)
                            aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCtt                                                                                                                      <  2:497397/151‑1 (MQ=255)
                                tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAACAGCTTGGGTTCTCACAAGAGCAGGGGTTTTTCACTTTTTGTGCtttt                                                                                                                    >  4:1004526/1‑149 (MQ=255)
                                 gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTTAAAAACTTTGAGTTTCACAAGAGAAGGGTTTTTTTTCTTTTTTTGTttttttt                                                                                                                 >  1:94604/1‑148 (MQ=255)
                                 gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtg                                                                                                                   <  4:65890/149‑1 (MQ=255)
                                  cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTCTTCTTAAACAAATTGGGGTTTCTCAAGAGGAGGCGTTTGTTACTTTTTGGGCtttttgt                                                                                                                 >  4:812790/1‑147 (MQ=255)
                                  cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgctg                                                                                                                <  4:731513/151‑1 (MQ=255)
                                   agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgctgc                                                                                                               <  4:150227/151‑1 (MQ=255)
                                   agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgctgc                                                                                                               <  2:1092838/151‑1 (MQ=255)
                                    gcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcg                                                                                                              <  2:1282859/151‑1 (MQ=255)
                                     caACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcgc                                                                                                             <  4:963668/151‑1 (MQ=255)
                                      aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTCTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCt                                                                                                            <  4:492260/151‑1 (MQ=255)
                                       aCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGGTGAAAAACTTTGGTTTTTACAAAAAGAGGCCTTTTTTACTTTTTGGGCTTTtgct                                                                                                                 >  3:440887/1‑145 (MQ=255)
                                          tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAAATTTGGTTTTCACAAAAGGAGGGTTTTGTTTCTTTTTGGGCGTTTTTTGCGCTgg                                                                                                          >  2:1365244/1‑149 (MQ=255)
                                           cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGt                                                                                                         <  4:710381/149‑1 (MQ=255)
                                           cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGt                                                                                                         <  2:340166/149‑1 (MQ=255)
                                           cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                                        <  2:622617/150‑1 (MQ=255)
                                             cccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTAt                                                                                                     <  2:101183/151‑1 (MQ=255)
                                              ccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAACATGGGTTCTCACAAAAAGAAGGTTTTTTTTCTTTTTGTGTTTTTTCTGCGCTgggttat                                                                                                     >  2:1366624/1‑145 (MQ=255)
                                               cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTTTTTAACAGCTTGGGTTTTCCCAAGAGCAGGCTTTTTTCACTTTTTGTGTTTTTTCTGCGCTGGGCTTTgg                                                                                                   >  1:345555/1‑151 (MQ=255)
                                                    ttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTc                                                                                               <  4:610842/150‑1 (MQ=255)
                                                     tCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCa                                                                                              <  4:1409138/150‑1 (MQ=255)
                                                      cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCAt                                                                                             <  4:1558361/150‑1 (MQ=255)
                                                      cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAACAGCTTTGGTTCTCACAAAAACAGGGTTTTTTGTATTTTTGTGCTTTTTTTGCGCTgggttttgggctcctt                                                                                            >  2:1369610/1‑137 (MQ=255)
                                                              tttttGTCGAGTTATGGGAGCTATTCGGCTACTACGGCGTTCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGTATTCTCTCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                                      <  4:1207139/149‑1 (MQ=255)
                                                              tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                                      <  4:779562/149‑1 (MQ=255)
                                                              tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAt                                                                                       <  3:381585/148‑1 (MQ=255)
                                                                  tGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTTGTTTAAAAGCTTTGGTTCTCCCAAAAGCAGGCGTTTTTTCCCTTTTGGGCTTTTTCTGCGCTGGTCTATTGGCTCCTTTCCATTtgcgg                                                                                 >  4:1017990/1‑150 (MQ=255)
                                                                   gTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTggcggc                                                                                <  2:1527098/150‑1 (MQ=255)
                                                                      gAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGCTTTTTTCGTTTAACAACTTGGGTTTTCACAAGAACCGGGGTTTTTTACTTTTTGTGATTTTTGTGCCCTGGGCTATGGCCCCCttttcctttgggg                                                                                 >  4:1105490/1‑136 (MQ=255)
                                                                       aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGt                                                                           <  4:1075229/151‑1 (MQ=255)
                                                                        gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTc                                                                          <  2:1072901/151‑1 (MQ=255)
                                                                        gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTc                                                                          <  3:55755/151‑1 (MQ=255)
                                                                          tATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCg                                                                         <  4:1139775/150‑1 (MQ=255)
                                                                           aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAACTTGGGTTTTCCCAAAAGCAGGGTTTTTTTACTTTTTGTGCTTTTTTTTCGCTGGTCTATTGCCTCCTTTTCATTtgcgggctttttgg                                                                        >  4:776516/1‑141 (MQ=255)
                                                                            tGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTGGAATCTCACAAAAGCAGGGTTTTTTCACTTTTTGTGGTTTTTCTGCGCTGGGCTTTGGCCTCCTTTTCATTggggggctttttggcg                                                                      >  3:453612/1‑140 (MQ=255)
                                                                             gggAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCg                                                                     <  2:1337490/150‑1 (MQ=255)
                                                                                            actacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTgggg                                                       <  2:1353211/149‑1 (MQ=255)
                                                                                            actacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGAc                                                     <  3:314543/151‑1 (MQ=255)
                                                                                               acGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCa                                                   <  2:496521/150‑1 (MQ=255)
                                                                                                 ggCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTAAACAGCCTTGATTCTCACAAGAGCAGGCGTTTGTCACTTTTTGTGCTTTTTCTGCGCTGGGCTATGGCCTCATTTCCATTGGGGGGTATGTCGGCGACCACCCCCTGGGGGCCcaa                                                 >  3:467741/1‑150 (MQ=255)
                                                                                                       cAGGGCGTACTGGCGGTTTTTTTCGTTAAAAAGCTTGGATTCTCCCAAAAGCCAGCTTTTGTCACTTTTTGTGGTTTTTCTTCGCTGGTCTATGGCCCCCTTTTCCTTGGCGGGCTTTTCGGCGGCCCACCTCCGGGGGCCCaacccccc                                            >  2:1015515/1‑144 (MQ=255)
                                                                                                       cAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCGATGTCGGCGACCACCTGCTGGGGACCAAACGCACCa                                          <  1:213286/151‑1 (MQ=255)
                                                                                                           gCGTACTGGCGGTTTTCTTTGTTTAAAAGCATTGGTTCTCCCAAGAGCAGGCGTTTGTTACTTTTTGGGGTTTTGTTGCGCCGGGCTATGGCCCCATTTTCATTTGCGGGGATTTTGGCGGCCACCCGCTGGGGGGCAcacccacc                                           >  2:629689/1‑146 (MQ=255)
                                                                                                            cGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGtt                                     <  4:1560922/151‑1 (MQ=255)
                                                                                                            cGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGtt                                     <  2:1606566/151‑1 (MQ=255)
                                                                                                               aCTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCtt                                  <  4:1000527/151‑1 (MQ=255)
                                                                                                                cTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCtt                                  >  4:571993/1‑150 (MQ=255)
                                                                                                                 tGGCGGTTTTCTTCGTTAAACCGCTTGGATTCTCACCAGAGCAGGCTTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGGTCTGTCGGCGGCCCCCTTCTGGGGGCCCAACACACCaattgtctttg                                >  2:795929/1‑141 (MQ=255)
                                                                                                                  ggCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAATGCACCATTGTTCTTGGa                               <  4:864072/151‑1 (MQ=255)
                                                                                                                  ggCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTg                                 <  3:41941/149‑1 (MQ=255)
                                                                                                                  ggCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGa                               <  4:1186106/151‑1 (MQ=255)
                                                                                                                   gCGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAg                              <  3:423440/151‑1 (MQ=255)
                                                                                                                    cGGTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGc                             >  4:1528412/1‑151 (MQ=255)
                                                                                                                       gTTTTCTTCGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCAc                           <  4:460424/151‑1 (MQ=255)
                                                                                                                           ttcttaGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtg                        >  3:19291/1‑151 (MQ=255)
                                                                                                                                cGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTg                   <  4:330366/151‑1 (MQ=255)
                                                                                                                                cGTTAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTg                   <  2:212477/151‑1 (MQ=255)
                                                                                                                                   ttAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCt                    >  2:1418145/1‑148 (MQ=255)
                                                                                                                                    tAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGc                 <  4:1676134/150‑1 (MQ=255)
                                                                                                                                    tAAACAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCg                >  2:944432/1‑151 (MQ=255)
                                                                                                                                     aaaCAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCg                <  4:921605/150‑1 (MQ=255)
                                                                                                                                     aaaCAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGa               <  3:321155/151‑1 (MQ=255)
                                                                                                                                     aaaCAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCCACCACCTGCTGGGGACCAAAAGCCCCACTGTTCTTGGAGCACTTGTGCTGGCg                >  2:709880/1‑150 (MQ=255)
                                                                                                                                       aCAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGAtt             <  1:126387/151‑1 (MQ=255)
                                                                                                                                       aCAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGAtt             >  4:1607554/1‑151 (MQ=255)
                                                                                                                                        cAGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTg            >  4:1299443/1‑151 (MQ=255)
                                                                                                                                         aGCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTg            <  2:811713/150‑1 (MQ=255)
                                                                                                                                          gCTTGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc          <  2:694749/151‑1 (MQ=255)
                                                                                                                                           cTTCGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTgg           >  1:445521/1‑149 (MQ=255)
                                                                                                                                            ttGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCt         >  1:228560/1‑150 (MQ=255)
                                                                                                                                            ttGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCt         >  2:15764/1‑150 (MQ=255)
                                                                                                                                             tGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTa        <  4:619359/150‑1 (MQ=255)
                                                                                                                                             tGGATTCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTa        <  4:84980/150‑1 (MQ=255)
                                                                                                                                                 ttCTCACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTg                                 <  1:296262/121‑1 (MQ=255)
                                                                                                                                                  tctcACAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCAt  >  4:591167/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCAT  >  NC_000913/3640858‑3641151

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