Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,002,214 A→G L32P (CTG→CCG)  yigI ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,2140AG99.6% 727.4 / ‑4.8 232L32P (CTG→CCG) yigIhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (111/120);  total (111/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

Only 100 of 244 total aligned reads displayed.

CCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCA  >  NC_000913/4002124‑4002303
                                                                                          |                                                                                         
ccgccgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGt                                                                                         >  1:719912/1‑93 (MQ=255)
ccgccgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGt                                                                                         >  1:4657901/1‑93 (MQ=255)
 cgccgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGtt                                                                                        >  1:216569/1‑93 (MQ=255)
  gccgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTc                                                                                       <  2:4125226/93‑1 (MQ=255)
    cgTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCa                                                                                     <  2:2535553/93‑1 (MQ=255)
      tGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACGCTTTTTCGTAACGCTCCGGTTCCATc                                                                                   >  2:3793204/1‑93 (MQ=255)
       gCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATcc                                                                                  >  2:3010639/1‑93 (MQ=255)
         aaaaTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATcccc                                                                                >  2:3287450/1‑93 (MQ=255)
          aaaTGCTTTGCGCGCAGTTGCCCACCATCATTGGCGGATTTTTAAAGGCCAGCTGTGCCATCTCTTTTTCGTGACGCTCCGGTTCCATCCCCa                                                                               <  2:1586956/93‑1 (MQ=255)
          aaaTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCa                                                                               >  2:598628/1‑93 (MQ=255)
           aaTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCaa                                                                              <  1:532690/93‑1 (MQ=255)
           aaTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCaa                                                                              <  1:3621945/93‑1 (MQ=255)
           aaTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCaa                                                                              <  1:1533646/93‑1 (MQ=255)
           aaTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCaa                                                                              <  1:3029155/93‑1 (MQ=255)
            aTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAAt                                                                             >  1:2985987/1‑93 (MQ=255)
            aTGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAAt                                                                             <  1:2026423/93‑1 (MQ=255)
             tGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATg                                                                            <  1:585283/93‑1 (MQ=255)
              gCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATgc                                                                           <  2:556191/93‑1 (MQ=255)
               cTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCGGCTGTGCGAACTCTTTTTCCTAACGCTCCGGTTCCATCCCCAATgcg                                                                          >  2:704188/1‑93 (MQ=255)
                 ttGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGcgcg                                                                        >  2:4150490/1‑93 (MQ=255)
                   gcgcCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGt                                                                      <  1:2097372/93‑1 (MQ=255)
                     gcCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCACTGCGCGGTTa                                                                    >  2:1472947/1‑93 (MQ=255)
                      cccAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTaa                                                                   >  2:1108848/1‑93 (MQ=255)
                       ccAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTaaa                                                                  >  1:398301/1‑93 (MQ=255)
                       ccAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTaaa                                                                  >  2:432016/1‑93 (MQ=255)
                        cAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAAt                                                                 >  2:2604980/1‑93 (MQ=255)
                         aGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATg                                                                <  2:3389389/93‑1 (MQ=255)
                         aGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATg                                                                <  2:1952098/93‑1 (MQ=255)
                          gTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATgg                                                               >  2:3238085/1‑93 (MQ=255)
                            tGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCa                                                             <  2:1027602/93‑1 (MQ=255)
                            tGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCa                                                             >  2:1878369/1‑93 (MQ=255)
                             gCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCAt                                                            >  2:1735639/1‑93 (MQ=255)
                               ccaccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCAtgt                                                          >  1:4279244/1‑93 (MQ=255)
                               ccaccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCAtgt                                                          <  2:4265720/93‑1 (MQ=255)
                               ccaccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCAtgt                                                          <  2:3103094/93‑1 (MQ=255)
                                caccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCAtgtg                                                         >  2:3186429/1‑93 (MQ=255)
                                 cccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGa                                                        <  2:2488969/92‑1 (MQ=255)
                                 accaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGa                                                        <  2:3550231/93‑1 (MQ=255)
                                 accaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGa                                                        >  1:3579130/1‑93 (MQ=255)
                                 accaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGa                                                        <  1:3389907/93‑1 (MQ=255)
                                  ccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGAt                                                       <  2:3015230/93‑1 (MQ=255)
                                  ccaTCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGAt                                                       >  1:1033388/1‑93 (MQ=255)
                                   catcatTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATa                                                      <  1:1506678/93‑1 (MQ=255)
                                     tcatTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaa                                                    >  2:1442576/1‑93 (MQ=255)
                                     tcatTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaa                                                    >  1:3544036/1‑93 (MQ=255)
                                       atTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaaca                                                  >  2:2791556/1‑93 (MQ=255)
                                        ttGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTTACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaacaa                                                 >  1:2517177/1‑93 (MQ=255)
                                        ttGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaacaa                                                 >  2:1787512/1‑93 (MQ=255)
                                        ttGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaacaa                                                 >  2:1520760/1‑93 (MQ=255)
                                        ttGGCTGATTTTAAAAGCCAGCTGTGCGGACTTTTTTTTGGAACGCTTCGGGTTCATTCCCCATGCGCGGTTTAATGGCATGTGATaaacaaa                                                >  2:4609223/1‑93 (MQ=255)
                                         tGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaacaaa                                                <  2:4422219/93‑1 (MQ=255)
                                         tGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaacaaa                                                <  1:1112202/93‑1 (MQ=255)
                                         tGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATaaacaaa                                                <  1:514998/93‑1 (MQ=255)
                                            cTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACAt                                             <  2:3648441/93‑1 (MQ=255)
                                             tGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACAtc                                            <  2:289119/93‑1 (MQ=255)
                                             tGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACAtc                                            <  1:3253133/93‑1 (MQ=255)
                                               aTTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACAtctc                                          >  2:1135115/1‑93 (MQ=255)
                                                tttttAAGGGCTGGCTGTGAGAACTTTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCa                                         <  2:84986/93‑1 (MQ=255)
                                                tttttAAATGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCa                                         >  2:4581981/1‑93 (MQ=255)
                                                tttttAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCa                                         <  1:102890/93‑1 (MQ=255)
                                                tttttAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCa                                         >  2:50653/1‑93 (MQ=255)
                                                 ttttAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCAc                                        <  2:3761678/93‑1 (MQ=255)
                                                  tttAAAGGCGATCTGTGCGAACTCTTTTTCCTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCAcc                                       <  2:3947498/93‑1 (MQ=255)
                                                  tttAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCAcc                                       <  1:2344131/93‑1 (MQ=255)
                                                    tAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGCTTAATGGCATGTGATAAACAAACATCTCACCCa                                     >  1:4133751/1‑93 (MQ=255)
                                                     aaaGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCAc                                    >  2:3135467/1‑93 (MQ=255)
                                                        ggCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTaa                                 <  2:3596123/93‑1 (MQ=255)
                                                        ggCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTaa                                 <  1:958192/93‑1 (MQ=255)
                                                         gCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAAt                                <  2:4001892/93‑1 (MQ=255)
                                                         gCCAGCTGTGCGAACTCTTTGTCGTATCTCTCAGGTTCAATCCCCAATGCGCGGTTAAATTGCATGTGATAACCAAACATCTCACCCACCAAt                                <  2:2747828/93‑1 (MQ=255)
                                                           cAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAAttt                              <  1:333608/93‑1 (MQ=255)
                                                           cAGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAAttt                              <  2:3737780/93‑1 (MQ=255)
                                                            aGCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTc                             >  2:2362238/1‑93 (MQ=255)
                                                             gCTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCa                            <  1:1081442/93‑1 (MQ=255)
                                                              cTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAg                           <  1:1676073/93‑1 (MQ=255)
                                                              cTGTGCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAg                           >  1:1280680/1‑93 (MQ=255)
                                                               tgtgCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGGTAAATGGCATGTGATAAACAAACATCTCACCCACCAATTTCAgg                          >  2:615932/1‑93 (MQ=255)
                                                                 tgCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGc                        <  1:4531098/93‑1 (MQ=255)
                                                                 tgCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGc                        <  1:1460113/93‑1 (MQ=255)
                                                                 tgCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGc                        >  1:1479763/1‑93 (MQ=255)
                                                                 tgCGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGc                        <  2:1971333/93‑1 (MQ=255)
                                                                   cGAACTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCtt                      <  1:2711085/93‑1 (MQ=255)
                                                                     aaCTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGt                    >  2:2439185/1‑93 (MQ=255)
                                                                     aaCTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGt                    >  1:4222410/1‑93 (MQ=255)
                                                                      aCTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGtt                   >  1:114302/1‑93 (MQ=255)
                                                                      aCTCTTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGtt                   <  1:2400912/93‑1 (MQ=255)
                                                                       ctctTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTc                  >  2:226081/1‑93 (MQ=255)
                                                                       ctctTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTc                  >  2:1640043/1‑93 (MQ=255)
                                                                         ctTTTTCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAg                >  1:603648/1‑93 (MQ=255)
                                                                          tttttCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGc               >  1:2634700/1‑93 (MQ=255)
                                                                              tCGTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTc           <  1:89311/93‑1 (MQ=255)
                                                                               cGTAACGGTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCa          <  1:2222326/93‑1 (MQ=255)
                                                                                gTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAg         >  1:1297813/1‑93 (MQ=255)
                                                                                gTAACGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAg         >  2:1192769/1‑93 (MQ=255)
                                                                                   aCGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTAc      <  1:3747600/93‑1 (MQ=255)
                                                                                   aCGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTAc      <  1:4355689/93‑1 (MQ=255)
                                                                                   aCGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTAc      <  1:2800947/93‑1 (MQ=255)
                                                                                    cGCTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCTCCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACg     <  1:3736599/93‑1 (MQ=255)
                                                                                      cTCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGc   >  2:3458900/1‑93 (MQ=255)
                                                                                       tCCGGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCa  >  2:3068652/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATCCCCAATGCGCGGTTAAATGGCATGTGATAAACAAACATCTCACCCACTAATTTCAGGGCTTGTTCAGCGGTCAGTACGGCA  >  NC_000913/4002124‑4002303

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