Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,368,057 G→T 5.1% A7S (GCC→TCC)  arnA → fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,368,0570GT 5.1% 252.9 / 2.5 73A7S (GCC→TCC) arnAfused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base T (2/2);  total (40/33)
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.39e-01

TATTGGCGCTCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGC  >  NC_000913/2367910‑2368197
                                                                                                                                                   |                                                                                                                                            
tATTGGCGCTCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGc                                                                                                                                             <  3:600165/149‑1 (MQ=255)
    ggCGCTCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCGGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTAc                                                                                                                                         >  1:705556/1‑149 (MQ=255)
      cgcTCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTTCACCGATGTGCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGATATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCa                                                                                                                                       <  4:549655/149‑1 (MQ=255)
         tCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTTCCTACCACGa                                                                                                                                    >  3:292901/1‑149 (MQ=255)
         tCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTTCCTACCACGa                                                                                                                                    >  1:342006/1‑149 (MQ=255)
             tttATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATg                                                                                                                                >  4:631495/1‑149 (MQ=255)
              ttATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATgg                                                                                                                               <  2:788639/149‑1 (MQ=255)
              ttATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATgg                                                                                                                               <  2:656309/149‑1 (MQ=255)
              ttATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATgg                                                                                                                               <  4:667949/149‑1 (MQ=255)
                aTCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGa                                                                                                                             <  2:786072/149‑1 (MQ=255)
                    gCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGcc                                                                                                                         >  3:148490/1‑149 (MQ=255)
                          gATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTa                                                                                                                   >  3:662681/1‑149 (MQ=255)
                           aTTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTAt                                                                                                                  >  2:27621/1‑149 (MQ=255)
                                tCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAg                                                                                                             >  4:324969/1‑149 (MQ=255)
                                  ggTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGcc                                                                                                           <  3:393656/149‑1 (MQ=255)
                                      aaTATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCctgc                                                                                                       <  4:266043/149‑1 (MQ=255)
                                        tatatCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCctgctg                                                                                                     >  4:756693/1‑149 (MQ=255)
                                          tatCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTgctggc                                                                                                   <  4:662681/149‑1 (MQ=255)
                                               gCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCg                                                                                              >  4:177933/1‑149 (MQ=255)
                                                cAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCgg                                                                                             <  1:567852/149‑1 (MQ=255)
                                                     tCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGc                                                                                                            >  4:781143/1‑129 (MQ=255)
                                                     tCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGc                                                                                                            <  3:781143/129‑1 (MQ=255)
                                                        acacCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGaaa                                                                                     >  2:720642/1‑149 (MQ=255)
                                                          acCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAAtt                                                                                   >  3:430314/1‑149 (MQ=255)
                                                           ccGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTa                                                                                  <  2:256068/149‑1 (MQ=255)
                                                             gATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAgc                                                                                <  4:565916/149‑1 (MQ=255)
                                                             gATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAgc                                                                                <  1:489808/149‑1 (MQ=255)
                                                                        ccgcccCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTAc                                                                     >  4:662991/1‑149 (MQ=255)
                                                                          gcccCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATgg                                                                                                                               >  2:499126/1‑89 (MQ=255)
                                                                          gcccCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATgg                                                                                                                               <  1:499126/89‑1 (MQ=255)
                                                                          gcccCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTAccc                                                                   <  3:662991/149‑1 (MQ=255)
                                                                            ccccGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCAt                                                                 >  1:88804/1‑149 (MQ=255)
                                                                              ccGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATAc                                                               >  2:401842/1‑149 (MQ=255)
                                                                              ccGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATAc                                                               >  1:704603/1‑149 (MQ=255)
                                                                                gCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTg                                                             >  1:73553/1‑149 (MQ=255)
                                                                                  tATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGAt                                                           >  1:633354/1‑149 (MQ=255)
                                                                                  tATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGAt                                                           >  1:515175/1‑149 (MQ=255)
                                                                                            aGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTg                                                 >  1:290698/1‑149 (MQ=255)
                                                                                             gCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAgc                                                                                <  4:269918/117‑1 (MQ=255)
                                                                                             gCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAgc                                                                                >  3:269918/1‑117 (MQ=255)
                                                                                                 gTTATCCGTCCATCCAGCAAGGAAAATGCATAATGAAAACCGTCGTTTTTTCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaaaa                                            <  4:292901/149‑1 (MQ=255)
                                                                                                 gTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTTCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaaaa                                            <  2:342006/149‑1 (MQ=255)
                                                                                                 gTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaaaa                                            <  1:27621/149‑1 (MQ=255)
                                                                                                    aTCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGcc                                         >  3:774731/1‑149 (MQ=255)
                                                                                                    aTCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGcc                                         >  2:801772/1‑149 (MQ=255)
                                                                                                        gTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCtttt                                     >  1:440146/1‑149 (MQ=255)
                                                                                                         tccatccaGCAAGTAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATCCTGATAATCCCGGTGAAAAAGCCTTTTa                                    <  4:69086/149‑1 (MQ=255)
                                                                                                         tccatccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTa                                    >  1:243531/1‑149 (MQ=255)
                                                                                                           catccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATg                                  <  2:73553/149‑1 (MQ=255)
                                                                                                            atccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATgg                                 >  1:789445/1‑149 (MQ=255)
                                                                                                               caGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGATATTAGCGCCATTTTTTCCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTc                              <  2:243531/149‑1 (MQ=255)
                                                                                                               caGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTc                              >  4:390611/1‑149 (MQ=255)
                                                                                                                aGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCg                             <  3:378582/149‑1 (MQ=255)
                                                                                                                   aaGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATc                                                       <  3:354653/120‑1 (MQ=255)
                                                                                                                   aaGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATc                                                       >  4:354653/1‑120 (MQ=255)
                                                                                                                      gAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCt                       <  1:801772/149‑1 (MQ=255)
                                                                                                                            gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGa                                                <  3:221181/118‑1 (MQ=255)
                                                                                                                            gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGa                                                >  4:221181/1‑118 (MQ=255)
                                                                                                                            gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTg                 >  1:261476/1‑149 (MQ=255)
                                                                                                                            gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCGGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTCATAATCCCGGTGAAAAATCCTTTTATGGTTCGGGGGCTCGtgtg                 <  2:265300/149‑1 (MQ=255)
                                                                                                                               taatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTAc                                                                     >  4:92998/1‑94 (MQ=255)
                                                                                                                               taatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTAc                                                                     <  3:92998/94‑1 (MQ=255)
                                                                                                                               taatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcg              >  2:80990/1‑149 (MQ=255)
                                                                                                                               taatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcg              >  1:93436/1‑149 (MQ=255)
                                                                                                                                 ctGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcggc            <  1:720642/148‑1 (MQ=255)
                                                                                                                                  tGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGt                                                  >  4:433318/1‑110 (MQ=255)
                                                                                                                                  tGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGt                                                  <  3:433318/110‑1 (MQ=255)
                                                                                                                                  tGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGgcggcg           >  1:77708/1‑149 (MQ=255)
                                                                                                                                  tGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGgcggcg           >  1:317889/1‑149 (MQ=255)
                                                                                                                                   gAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGcggcgg          >  1:259556/1‑149 (MQ=255)
                                                                                                                                      aaCCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGaaa       >  2:714495/1‑149 (MQ=255)
                                                                                                                                        ccGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAaga     <  2:88804/149‑1 (MQ=255)
                                                                                                                                           tcgtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGc  <  2:178226/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TATTGGCGCTCAGTTTATCGGCATGGGATTACTCGGTGAATATATCGGCAGGATCTACACCGATGTCCGCGCCCGCCCCCGCTATTTTGTTCAGCAAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGAGGC  >  NC_000913/2367910‑2368197

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