Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,591,413 A→G 5.2% N56S (AAC→AGC)  yffB → putative ArsC family reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,591,4130AG 5.2% 239.5 / 2.8 77N56S (AAC→AGC) yffBputative ArsC family reductase
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (2/2);  total (39/38)
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.15e-01

AAAAATTGTGACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGG  >  NC_000913/2591268‑2591560
                                                                                                                                                 |                                                                                                                                                   
aaaaaTTGTGACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAacac                                                                                                                                                  >  7:134694/1‑149 (MQ=255)
aaaaaTTGTGACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAacac                                                                                                                                                  >  7:369910/1‑149 (MQ=255)
aaaaaTTGTGACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAAGTATTGAACGATTTTATCAACGAGTTAGGCTGGGAAGCGTTACTCAacac                                                                                                                                                  >  8:283141/1‑149 (MQ=255)
 aaaaTTGTGACGCCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACAcc                                                                                                                                                 >  7:259318/1‑149 (MQ=255)
 aaaaTTGTGACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACAcc                                                                                                                                                 >  3:358402/1‑149 (MQ=255)
        tgACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGt                                                                                                                                           <  3:331746/148‑1 (MQ=255)
        tgACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTa                                                                                                                                          <  2:333078/149‑1 (MQ=255)
           cacCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACaa                                                                                                                                       <  1:317429/149‑1 (MQ=255)
            acCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAAc                                                                                                                                      >  1:348008/1‑149 (MQ=255)
             ccATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAAcc                                                                                                                                     >  8:165004/1‑149 (MQ=255)
              cATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGACGCGTTACTCAACACCCGTGGTACAACCt                                                                                                                                    >  4:138110/1‑149 (MQ=255)
              cATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTAACACAACACCCGTGGTACACCCt                                                                                                                                    >  2:155413/1‑149 (MQ=255)
               aTTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAGCACCCGTGGTACAACCTg                                                                                                                                   >  7:287735/1‑149 (MQ=255)
               aTTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAGCACCCGTGGTACAACCTg                                                                                                                                   >  7:282321/1‑149 (MQ=255)
                ttAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTgg                                                                                                                                  >  6:66958/1‑149 (MQ=255)
                  aaaaaGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGc                                                                                                                                 <  1:336244/148‑1 (MQ=255)
                    aaaGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGt                                                                                                                                           <  3:120203/136‑1 (MQ=255)
                    aaaGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGt                                                                                                                                           >  4:120203/1‑136 (MQ=255)
                           cGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGa                                                                                                                       >  1:234785/1‑149 (MQ=255)
                            gCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGa                                                                                                                       <  5:66958/148‑1 (MQ=255)
                            gCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGAc                                                                                                                      >  3:163732/1‑149 (MQ=255)
                              cGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGa                                                                                                                       >  1:380492/1‑146 (MQ=255)
                              cGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGa                                                                                                                       <  2:380492/146‑1 (MQ=255)
                               gTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGaa                                                                                                                   <  7:196018/149‑1 (MQ=255)
                               gTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGaa                                                                                                                   <  4:163732/149‑1 (MQ=255)
                                         aGCCAATAACATCGACTATCGATTTAATGATTACCCCTTCGATGGTCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGTGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCa                                                                                                         <  1:310558/149‑1 (MQ=255)
                                           ccAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAAt                                                                                                       <  5:192439/149‑1 (MQ=255)
                                             aaTAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACTATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATa                                                                                                      >  1:157137/1‑148 (MQ=255)
                                               tAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGACACCACCCGCAATaaaa                                                                                                   >  3:185843/1‑149 (MQ=255)
                                                 aCATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAAt                                                                                                  <  2:348008/148‑1 (MQ=255)
                                                   aTCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCAc                                                                                               <  8:134694/149‑1 (MQ=255)
                                                      gACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGa                                                                                            <  8:259318/149‑1 (MQ=255)
                                                       aCTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGAt                                                                                           <  8:19309/149‑1 (MQ=255)
                                                             gTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGcc                                                                                     <  4:142811/149‑1 (MQ=255)
                                                             gTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGcc                                                                                     >  8:394369/1‑149 (MQ=255)
                                                                ttCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCctct                                                                                  >  7:176401/1‑149 (MQ=255)
                                                                  cATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATGAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCTATGCGGCCTCTgc                                                                                >  1:129866/1‑149 (MQ=255)
                                                                   aTGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTgcg                                                                               >  7:255982/1‑149 (MQ=255)
                                                                     gATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTgcggc                                                                             >  1:234473/1‑149 (MQ=255)
                                                                     gATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTgcggc                                                                             >  8:346661/1‑149 (MQ=255)
                                                                     gATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTgcggc                                                                             >  7:183908/1‑149 (MQ=255)
                                                                      aTTACCGCGTCGATGGGATGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATTCGGCCTCTgcggcg                                                                            >  8:273622/1‑149 (MQ=255)
                                                                         aCCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCggcggc                                                                          <  8:319292/148‑1 (MQ=255)
                                                                          ccGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCAt                                                                        <  2:300098/149‑1 (MQ=255)
                                                                           cgcgTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCAtt                                                                       <  1:200142/149‑1 (MQ=255)
                                                                             cgTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAGCACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTaa                                                                     <  8:287735/149‑1 (MQ=255)
                                                                             cgTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAGCACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTaa                                                                     <  8:282321/149‑1 (MQ=255)
                                                                              gTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAAt                                                                    <  1:181902/149‑1 (MQ=255)
                                                                              gTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCCGCCGCAGAGGCCGCAtcggt                                                                    <  2:234785/149‑5 (MQ=255)
                                                                                  aTGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACt                                                                <  8:255982/149‑1 (MQ=255)
                                                                                   tGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTg                                                               <  8:176401/149‑1 (MQ=255)
                                                                                      gCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGaaa                                                            <  8:55992/149‑1 (MQ=255)
                                                                                          gACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGc                                                         <  1:155413/148‑1 (MQ=255)
                                                                                               cGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCaa                                                   <  1:354712/149‑1 (MQ=255)
                                                                                                   ttattGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTAt                                               <  7:165004/149‑1 (MQ=255)
                                                                                                     attGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATc                                              <  2:53364/148‑1 (MQ=255)
                                                                                                     attGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATc                                              >  4:309739/1‑148 (MQ=255)
                                                                                                      ttGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAACCCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCaa                                            <  8:369910/149‑1 (MQ=255)
                                                                                                         aaCGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACg                                         >  6:385293/1‑149 (MQ=255)
                                                                                                          aCGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGt                                        >  3:248842/1‑149 (MQ=255)
                                                                                                           cGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTc                                       >  3:237031/1‑149 (MQ=255)
                                                                                                              ttttATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCAt                                    >  3:152670/1‑149 (MQ=255)
                                                                                                                     aaCGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTg                             >  7:293624/1‑149 (MQ=255)
                                                                                                                     aaCGAATTAGGCTGGGAAGCGTTACTCAACACCCGTCGTACAACCTGGCCTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCACTGCGCCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCAATGctctc                             >  2:343997/1‑148 (MQ=255)
                                                                                                                      aCGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGc                            <  3:138110/149‑1 (MQ=255)
                                                                                                                      aCGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGCAAACTGTACGAAACCACCCGCAATAAAATCACCGATGCGGCCGCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGc                            <  7:273622/149‑1 (MQ=255)
                                                                                                                             aGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCg                     >  5:177827/1‑149 (MQ=255)
                                                                                                                              ggCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCgg                    >  7:226047/1‑149 (MQ=255)
                                                                                                                               gCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGt                   <  4:358402/149‑1 (MQ=255)
                                                                                                                                 tGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGTCATCTGCGGCGGCATTAATGACTGAAATCACAGCAATTATCAAACGTCCATTGCTCTCCGTGCCAGGt                   >  6:26403/1‑147 (MQ=255)
                                                                                                                                 tGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCg                     <  1:275571/145‑1 (MQ=255)
                                                                                                                                 tGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCg                     >  2:275571/1‑145 (MQ=255)
                                                                                                                                 tGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGt                   >  1:138246/1‑147 (MQ=255)
                                                                                                                                   ggAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTaa                 <  3:309739/147‑1 (MQ=255)
                                                                                                                                   ggAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGc               <  7:394369/149‑1 (MQ=255)
                                                                                                                                    gAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGcc              >  5:381001/1‑149 (MQ=255)
                                                                                                                                         gTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTAtgc         >  1:1918/1‑149 (MQ=255)
                                                                                                                                          ttACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTAtgtt        >  3:210322/1‑147 (MQ=255)
                                                                                                                                              tCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGctgctg    >  4:237655/1‑149 (MQ=255)
                                                                                                                                               cAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTgg   <  5:385293/149‑1 (MQ=255)
                                                                                                                                                aaCACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTggg  >  3:18553/1‑149 (MQ=255)
                                                                                                                                                aaCACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTggg  >  8:11648/1‑149 (MQ=255)
                                                                                                                                                aaCACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTggg  >  1:72232/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAAAATTGTGACACCATTAAAAAGGCTCGCCGTTGGCTGGAAGCCAATAACATCGACTATCGTTTTCATGATTACCGCGTCGATGGGCTGGACAGCGAATTATTGAACGATTTTATCAACGAATTAGGCTGGGAAGCGTTACTCAACACCCGTGGTACAACCTGGCGTAAACTGGACGAAACCACCCGCAATAAAATCACCGATGCGGCCTCTGCGGCGGCATTAATGACTGAAATGCCTGCAATTATCAAACGTCCATTGCTCTGCGTGCCCGGTAAGCCTATGCTGCTGGG  >  NC_000913/2591268‑2591560

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