Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,205,438 T→C intergenic (‑328/‑287) purH ← / → rrsE phosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,205,4380TC99.4% 555.1 / ‑5.5 177intergenic (‑328/‑287)purH/rrsEphosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (86/90);  total (87/90)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 184 total aligned reads displayed.

TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAG  >  NC_000913/4205348‑4205512
                                                                                          |                                                                             
taAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCcccc                                                                             >  1:2341053/1‑93 (MQ=255)
 aaaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCa                                                                            >  1:2242260/1‑93 (MQ=255)
  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCAt                                                                           <  1:2704921/93‑1 (MQ=255)
  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCAt                                                                           >  2:3155830/1‑93 (MQ=255)
      gggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGac                                                                       >  1:593724/1‑93 (MQ=255)
      gggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGac                                                                       >  1:682065/1‑93 (MQ=255)
       ggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGTACTCCCTATAATGCGCCCCCATCGaca                                                                      >  2:2422944/1‑93 (MQ=255)
       ggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGaca                                                                      <  2:2050241/93‑1 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                    <  1:3754681/93‑1 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                    <  2:4418071/93‑1 (MQ=255)
        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                    <  1:1838248/93‑1 (MQ=255)
         cTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcg                                                                   <  1:206355/93‑1 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                  <  1:3967961/93‑1 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                  <  2:1128253/93‑1 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                  >  2:2964818/1‑93 (MQ=255)
          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                  <  2:4430943/93‑1 (MQ=255)
           tttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggc                                                                 <  1:1716596/93‑1 (MQ=255)
            ttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcg                                                                <  1:2978433/93‑1 (MQ=255)
             tGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcgg                                                               >  1:4267073/1‑93 (MQ=255)
              tcTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGa                                                              <  1:746582/92‑1 (MQ=255)
              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGa                                                              <  1:2531846/93‑1 (MQ=255)
                tgaggAATGCAGCATTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCTCCCCCATCGACACGGCGGAtg                                                            <  1:3658108/90‑1 (MQ=38)
                tGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGGGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtg                                                            >  1:4622820/1‑93 (MQ=255)
                tGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtg                                                            <  1:3115507/93‑1 (MQ=255)
                tGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtg                                                            <  2:959687/93‑1 (MQ=255)
                 ggggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgt                                                           >  2:2006605/1‑93 (MQ=255)
                 ggggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgt                                                           <  2:174282/93‑1 (MQ=255)
                  gggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgtg                                                          >  2:77344/1‑93 (MQ=255)
                  gggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgtg                                                          >  1:532449/1‑93 (MQ=255)
                     aaTGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAAt                                                       <  2:2177156/93‑1 (MQ=255)
                     aaTGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAAt                                                       <  2:4742058/93‑1 (MQ=255)
                         cagcagTTAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACt                                                  <  1:1377643/93‑1 (MQ=255)
                          agcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACtt                                                 <  2:2734526/93‑1 (MQ=255)
                           gcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                                <  1:783316/93‑1 (MQ=255)
                           gcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                                <  2:4315677/93‑1 (MQ=255)
                           gcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                                <  2:3991063/93‑1 (MQ=255)
                            gagtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTca                                               <  2:293184/92‑1 (MQ=255)
                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                              <  1:4229649/93‑1 (MQ=255)
                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                              >  1:3596276/1‑93 (MQ=255)
                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                              <  2:2342527/93‑1 (MQ=255)
                              gtcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcaca                                             <  2:2904816/93‑1 (MQ=255)
                              atcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcaca                                             <  1:2713446/92‑1 (MQ=255)
                               tcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCTCTTcacac                                            <  1:3614225/93‑1 (MQ=255)
                                cagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacaca                                           >  2:3135668/1‑93 (MQ=255)
                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGTGAACCCATCGACACGGCGGGAGTGAATCACTTCACGCaa                                          >  1:4433958/1‑93 (MQ=25)
                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                          >  1:2999262/1‑93 (MQ=255)
                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                          >  2:1788870/1‑93 (MQ=255)
                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                          >  2:4196934/1‑93 (MQ=255)
                                  gtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGGACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaaa                                         >  2:198748/1‑93 (MQ=255)
                                  gtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaaa                                         <  2:456441/93‑1 (MQ=255)
                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAc                                        <  1:4295420/93‑1 (MQ=255)
                                    caTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACa                                      <  2:619758/93‑1 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCTGCCTGCGGAGAACTCCCTATAATGCGCCCCCATGGACACGGCGGATGTGAATCACTTCACACAAACAGc                                    <  1:4705008/93‑1 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                    <  2:4124948/93‑1 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                    <  2:3102047/93‑1 (MQ=255)
                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                    <  1:3556423/93‑1 (MQ=255)
                                        tttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcg                                  <  1:2496206/93‑1 (MQ=255)
                                         ttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcgg                                 >  1:95995/1‑93 (MQ=255)
                                          tCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggt                                <  1:6275/93‑1 (MQ=255)
                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                               <  1:2184863/93‑1 (MQ=255)
                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                               <  2:194847/93‑1 (MQ=255)
                                            tGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttc                              >  1:3416703/1‑93 (MQ=255)
                                             gCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcg                             >  1:2201250/1‑93 (MQ=255)
                                              cAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcgg                            >  2:4458594/1‑93 (MQ=255)
                                               aaTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTTACACAAACAGCcgggtcggt                           >  2:4645962/1‑93 (MQ=255)
                                               aaTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTGACACAAACAACcggttcggt                           >  1:518408/1‑93 (MQ=255)
                                               aaTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggt                           >  2:2262341/1‑93 (MQ=255)
                                                  ttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGa                        <  1:4117941/93‑1 (MQ=255)
                                                  ttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGa                        <  2:1371214/93‑1 (MQ=255)
                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                       <  1:118462/93‑1 (MQ=255)
                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                       <  2:2835885/93‑1 (MQ=255)
                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                       <  1:334185/93‑1 (MQ=255)
                                                    ttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAag                      >  2:2916092/1‑93 (MQ=255)
                                                      cTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAagag                    >  2:4117828/1‑93 (MQ=255)
                                                      cTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCACGACACCGGCGGTGTGAATCACCTTCACCAAACAACCCGGTCGGGTTGaaaga                    >  1:186719/1‑90 (MQ=12)
                                                       tATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAgaga                   >  2:3702360/1‑93 (MQ=255)
                                                       tATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAgaga                   >  2:316484/1‑93 (MQ=255)
                                                       tATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAgaga                   >  1:692163/1‑93 (MQ=255)
                                                       tATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAgaga                   >  1:3236130/1‑93 (MQ=255)
                                                       tATTGCGGCCTGCGGAGAACTCCCTATAATGCCCCCCCATCGACACGGCGGATGTGAATCACTTCACACCAACAACCGGTTCGGTTGAAgaga                   >  1:2811946/1‑93 (MQ=17)
                                                        aTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                  >  1:2933760/1‑93 (MQ=255)
                                                            cGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAAt              >  2:587347/1‑93 (MQ=255)
                                                            cGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAAt              >  2:421733/1‑93 (MQ=255)
                                                            cGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAAt              <  1:3570102/93‑1 (MQ=255)
                                                             ggCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGCTCGGTTGAAGAGAAAAATc             >  1:1430394/1‑93 (MQ=255)
                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc            >  2:1712754/1‑93 (MQ=255)
                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc            <  1:1058970/93‑1 (MQ=255)
                                                               ccTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCt           >  1:1171748/1‑93 (MQ=255)
                                                               ccTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCt           >  2:4270333/1‑93 (MQ=255)
                                                                 tGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGa         <  2:665807/93‑1 (MQ=255)
                                                                   cGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGCATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGaaa       >  1:828397/1‑93 (MQ=255)
                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAtt     <  2:2450640/93‑1 (MQ=21)
                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAtt     >  2:2745505/1‑93 (MQ=21)
                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAtt     >  1:193011/1‑93 (MQ=21)
                                                                      agaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc    <  2:692858/93‑1 (MQ=18)
                                                                      agaACTCCCTATAATGCGCCCCCATCGACACGGAGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc    >  1:1693688/1‑93 (MQ=16)
                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa   >  2:2571304/1‑93 (MQ=17)
                                                                        aaCTCCCTCTAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAg  >  1:2619031/1‑93 (MQ=9)
                                                                        aaCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAg  <  1:1655304/93‑1 (MQ=14)
                                                                        aaCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAg  >  2:91471/1‑93 (MQ=14)
                                                                                          |                                                                             
TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAG  >  NC_000913/4205348‑4205512

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