Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,421 2 bp→GC coding (381‑382/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4210CG100.0% 1091.3 / NA 346G128R (GGA→CGA) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (180/166);  total (180/166)
*NC_0009134,001,4220GC99.7% 1084.2 / NA 346D127E (GAC→GAGrarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (180/165);  total (180/166)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508
                                                                                       ||                                                                                      
taataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAAtcgctcgc                                                                                    >  2:4188866/1‑93 (MQ=255)
   taGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                 >  1:7515213/1‑93 (MQ=255)
   taGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                 >  2:5579401/1‑93 (MQ=255)
   agggTAGAAAACCAAAAGCCCAAAGCTGGATGAACACGCCACATATCGCCAGAACCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                 <  2:206149/91‑1 (MQ=255)
         gCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                                           <  1:5115620/93‑1 (MQ=255)
         gCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                                           >  1:2166002/1‑93 (MQ=255)
            aaCCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatca                                                                        <  1:56946/93‑1 (MQ=255)
               cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                                     >  2:1273706/1‑93 (MQ=255)
                 aaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                                                   <  2:3540680/93‑1 (MQ=255)
                 aaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                                                   >  1:7588485/1‑93 (MQ=255)
                 aaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                                                   >  2:4969744/1‑93 (MQ=255)
                 aaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                                                   <  2:538136/93‑1 (MQ=255)
                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  1:7179902/93‑1 (MQ=255)
                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  2:291792/93‑1 (MQ=255)
                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  2:7299948/93‑1 (MQ=255)
                    gTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGc                                                                >  1:1911532/1‑93 (MQ=255)
                    gTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGc                                                                >  2:673772/1‑93 (MQ=255)
                    gTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGc                                                                <  1:6535590/93‑1 (MQ=255)
                     tCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGca                                                               <  2:4411305/93‑1 (MQ=255)
                         caGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAAt                                                           >  1:2958660/1‑93 (MQ=255)
                          cgCTGGACTACACCGCCACATATCGCACGACTCACAGCCAGCCGTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcacagta                                                          <  2:7013304/92‑4 (MQ=255)
                           gCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                                         <  1:4599382/93‑1 (MQ=255)
                            cTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                                        <  2:3476843/93‑1 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                                      >  2:1611571/1‑93 (MQ=255)
                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCTCAGCACAATGTTCAc                                                     >  2:7834107/1‑93 (MQ=255)
                                aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                                    >  1:6623677/1‑93 (MQ=255)
                                aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                                    >  2:4318925/1‑93 (MQ=255)
                                aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                                    >  1:122/1‑93 (MQ=255)
                                 cTACCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   <  2:1500246/93‑1 (MQ=255)
                                 cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   >  2:6395319/1‑93 (MQ=255)
                                 cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   <  1:22907/93‑1 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  >  1:104405/1‑93 (MQ=255)
                                    acacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                                                <  1:1093744/93‑1 (MQ=255)
                                    acacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCg                                                <  1:1643468/93‑1 (MQ=255)
                                     cacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCCGCACAATGTTCACCAGCgg                                               >  1:6409595/1‑93 (MQ=255)
                                      ccGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                                              <  1:1706370/92‑1 (MQ=255)
                                       cGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGt                                             <  2:1114334/93‑1 (MQ=255)
                                        gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            >  1:5542588/1‑93 (MQ=255)
                                        gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            >  1:921381/1‑93 (MQ=255)
                                         ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                                           >  2:7641265/1‑93 (MQ=255)
                                         ccACATATCCCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                                           <  1:1302290/93‑1 (MQ=255)
                                          cacaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                                          >  2:1825574/1‑93 (MQ=255)
                                           acaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaat                                         <  2:3256232/93‑1 (MQ=255)
                                            caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaata                                        <  1:6058328/93‑1 (MQ=255)
                                             atatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaataa                                       >  1:1841427/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      >  1:5891429/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      >  2:6307091/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  2:1409421/93‑1 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      >  1:4971615/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  1:7077539/93‑1 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAACGTTCACCAGCGGGTTAATaaa                                      >  2:2413823/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCACCAGCACAATGTTCACCAGCGGGATAATaaa                                      >  1:3106472/1‑93 (MQ=255)
                                                tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                                    <  2:3003539/93‑1 (MQ=255)
                                                 cGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTa                                   >  1:2557019/1‑93 (MQ=255)
                                                  gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa                                  >  2:7651347/1‑93 (MQ=255)
                                                  gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa                                  >  2:3778193/1‑93 (MQ=255)
                                                   ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc                                 >  1:4652744/1‑93 (MQ=255)
                                                   ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc                                 <  2:1645769/93‑1 (MQ=255)
                                                    cAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc                                <  2:6949146/93‑1 (MQ=255)
                                                    cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc                                >  2:7104085/1‑93 (MQ=255)
                                                      gAATCACGCCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              <  1:68362/93‑1 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              >  1:5732863/1‑93 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              >  2:2634264/1‑93 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              <  2:2572073/93‑1 (MQ=255)
                                                       aaTCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             <  1:4418078/93‑1 (MQ=255)
                                                       aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             >  1:2476795/1‑93 (MQ=255)
                                                        aTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg                            >  2:5290981/1‑93 (MQ=255)
                                                        aTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg                            >  1:7811970/1‑93 (MQ=255)
                                                        aTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCACCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg                            >  1:3231803/1‑93 (MQ=255)
                                                         tCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGc                           <  2:3934769/93‑1 (MQ=255)
                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          >  1:4306067/1‑93 (MQ=255)
                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          >  1:27194/1‑93 (MQ=255)
                                                           cccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         <  2:6202622/92‑1 (MQ=255)
                                                           aCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         >  1:4184320/1‑93 (MQ=255)
                                                              gccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct                      <  1:6238893/93‑4 (MQ=255)
                                                              cccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct                      <  2:4756308/92‑4 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   <  1:3565516/93‑1 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  1:667110/1‑93 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  2:4826994/1‑93 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  1:6936808/1‑93 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  1:397915/1‑93 (MQ=255)
                                                                  gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                  <  2:3729204/93‑1 (MQ=255)
                                                                  gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                  >  1:2357066/1‑93 (MQ=255)
                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 <  1:988027/93‑1 (MQ=255)
                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 >  2:3348898/1‑93 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               >  2:2232555/1‑93 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               >  2:7156396/1‑93 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               <  2:5090724/93‑1 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               <  2:6666154/93‑1 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               >  1:5313628/1‑93 (MQ=255)
                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGcct              >  1:2835472/1‑91 (MQ=255)
                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  2:6883392/1‑93 (MQ=255)
                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  2:4603184/1‑93 (MQ=255)
                                                                        gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat            >  1:4557894/1‑93 (MQ=255)
                                                                         cATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg           <  2:5979766/93‑1 (MQ=255)
                                                                           tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATggg         >  1:5620821/1‑92 (MQ=255)
                                                                             cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa       >  2:4057772/1‑93 (MQ=255)
                                                                              ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt      >  1:4914198/1‑93 (MQ=255)
                                                                              ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt      >  2:2838808/1‑93 (MQ=255)
                                                                                  gAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGGGGGTTAATAAAGTAACGAAGGCTCGCTTCCAGCATATCGTGGTtttt  >  2:5707257/1‑93 (MQ=255)
                                                                                       ||                                                                                      
TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508

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: