Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,016,797 A→G Y119C (TAT→TGT)  ubiE → UbiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,016,7970AG99.7% 1211.7 / ‑5.6 379Y119C (TAT→TGT) ubiEUbiE
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (182/196);  total (182/197)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 393 total aligned reads displayed.

GCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGAC  >  NC_000913/4016707‑4016885
                                                                                          |                                                                                        
   aaGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTc                                                                                     >  1:6142692/1‑93 (MQ=255)
   aaGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTc                                                                                     <  1:6629949/93‑1 (MQ=255)
    aGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCa                                                                                    <  1:585009/93‑1 (MQ=255)
    aGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCa                                                                                    <  2:3347466/93‑1 (MQ=255)
    aGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCa                                                                                    <  1:2645410/93‑1 (MQ=255)
    aGTGGTCCTTGCTGATATCAATGAATCCATACTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCa                                                                                    <  1:1842461/93‑1 (MQ=255)
     gTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAg                                                                                   <  1:6697979/93‑1 (MQ=255)
      tGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAgg                                                                                  <  2:1328592/93‑1 (MQ=255)
      tGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAgg                                                                                  >  1:6854305/1‑93 (MQ=255)
       ggTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGc                                                                                 <  1:1764494/93‑1 (MQ=255)
        gtgcTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCg                                                                                <  1:1010444/90‑1 (MQ=255)
        gTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCg                                                                                <  2:4583789/93‑1 (MQ=255)
          ccTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGaa                                                                              <  2:5652774/93‑1 (MQ=255)
           cTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAAc                                                                             >  1:7536457/1‑93 (MQ=255)
             tGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGc                                                                           <  1:406618/93‑1 (MQ=255)
             tGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGc                                                                           >  1:3517276/1‑93 (MQ=255)
             tGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGc                                                                           <  2:3157252/93‑1 (MQ=255)
              gCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCt                                                                          >  2:5000917/1‑93 (MQ=255)
              gCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCt                                                                          <  2:7850469/93‑1 (MQ=255)
               cTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTg                                                                         >  2:2488750/1‑93 (MQ=255)
                tGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGa                                                                        <  1:1047761/93‑1 (MQ=255)
                tGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGa                                                                        <  1:1664494/93‑1 (MQ=255)
                tGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGa                                                                        <  1:363392/93‑1 (MQ=255)
                tGATATCAATGAATCCATGCTCAAAAGGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGa                                                                        <  1:2119722/93‑1 (MQ=255)
                 gATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGAGAACGCTGAg                                                                       >  2:4447388/1‑93 (MQ=255)
                  atatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAgg                                                                      <  1:943957/93‑1 (MQ=255)
                   tatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGgc                                                                     >  2:3881879/1‑93 (MQ=255)
                   tatCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGgc                                                                     >  2:7630060/1‑93 (MQ=255)
                     tCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGgcgc                                                                   >  1:4673214/1‑93 (MQ=255)
                      cAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCt                                                                  >  2:7262270/1‑93 (MQ=255)
                       aaTGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCta                                                                 >  2:3527758/1‑92 (MQ=255)
                        aTGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCtac                                                                >  1:2497163/1‑91 (MQ=255)
                            aTCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGtt                                                            <  2:1855135/93‑1 (MQ=255)
                             tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                           >  2:806459/1‑93 (MQ=255)
                             tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                           >  1:7461036/1‑93 (MQ=255)
                             tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                           >  1:6661635/1‑93 (MQ=255)
                             tCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTc                                                           >  1:530012/1‑93 (MQ=255)
                               cATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTccc                                                         <  2:3671358/93‑1 (MQ=255)
                                 tGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCgg                                                       >  1:3479989/1‑93 (MQ=255)
                                 tGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCgg                                                       >  1:7188068/1‑93 (MQ=255)
                                  gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGa                                                      <  1:1301058/90‑1 (MQ=255)
                                  gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGa                                                      <  2:5568942/90‑1 (MQ=255)
                                  gctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGCGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGa                                                      >  1:7378423/4‑93 (MQ=255)
                                   ctcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGAt                                                     >  2:1237945/3‑93 (MQ=255)
                                    tcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATa                                                    >  2:4286017/2‑93 (MQ=255)
                                    tcAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATa                                                    >  2:4868416/2‑93 (MQ=255)
                                     cAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATaa                                                   >  2:832124/1‑93 (MQ=255)
                                      aaaaTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAac                                                  >  2:6494115/1‑93 (MQ=255)
                                         aggggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACAcc                                               <  2:109622/91‑1 (MQ=255)
                                         aTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACAcc                                               >  2:6286588/1‑93 (MQ=255)
                                         aTGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACAcc                                               >  1:553528/1‑93 (MQ=255)
                                          tGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCt                                              >  2:6315312/1‑93 (MQ=255)
                                          tGGGACGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCt                                              <  1:3134606/93‑1 (MQ=255)
                                            ggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCttt                                            <  2:1700925/93‑1 (MQ=255)
                                            ggCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCttt                                            <  2:4186778/93‑1 (MQ=255)
                                             gCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTg                                           <  2:3382096/93‑1 (MQ=255)
                                              ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                          <  2:4665569/93‑1 (MQ=255)
                                              ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                          <  2:2119249/93‑1 (MQ=255)
                                              ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                          <  1:762314/93‑1 (MQ=255)
                                              ccGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGa                                          >  1:5189585/1‑93 (MQ=255)
                                                   agaAGCGGCGTAATATCGGTGTGATTGGCAACGTTGAGTTTGTTCAGGCGAACGCTCACGCGCTCCCGTTCCCGGAGAACGCCTTTGATGGCa                                     <  2:7127261/93‑1 (MQ=255)
                                                    gaAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCAt                                    >  1:1384580/1‑93 (MQ=255)
                                                         tGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCa                               <  1:2319523/93‑1 (MQ=255)
                                                         tGCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCAGGATAACACCTTTGATTGCATCACCa                               <  1:335031/93‑1 (MQ=255)
                                                          gCGTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTATCGTTCCCGGATAACACCTTTGATTGCATCACCAt                              <  2:6404473/93‑1 (MQ=255)
                                                           cGTAATACCGGTGTGATTGGCAACGTTGAGTGTGTTCAGCCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAtt                             <  1:5023261/93‑1 (MQ=255)
                                                            gTAATATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                            >  2:571029/1‑93 (MQ=255)
                                                              aaTATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCg                          >  2:375778/1‑93 (MQ=255)
                                                              aaTATCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCg                          >  1:3285293/1‑93 (MQ=255)
                                                               atatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGt                         <  1:4524548/93‑1 (MQ=255)
                                                                tatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGtt                        >  2:4071765/1‑93 (MQ=255)
                                                                tatCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGtt                        >  2:2241513/1‑93 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  2:6897898/93‑1 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  1:1218839/93‑1 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       >  2:6783133/1‑93 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  1:4570056/93‑1 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  1:4792253/93‑1 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  1:3794666/93‑1 (MQ=255)
                                                                 atCGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGttt                       <  2:2328636/93‑1 (MQ=255)
                                                                   cGGTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTgg                     <  2:6984801/93‑1 (MQ=255)
                                                                    ggTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTggc                    >  2:2641488/1‑92 (MQ=255)
                                                                    ggTGTGATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGt                    <  2:283846/93‑1 (MQ=255)
                                                                      tgtgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCt                  <  2:5518736/93‑1 (MQ=255)
                                                                      gttgATTGGCAAAGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCt                  <  2:1508813/91‑1 (MQ=255)
                                                                        tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                >  2:5464622/1‑93 (MQ=255)
                                                                        tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                >  2:4334304/1‑93 (MQ=255)
                                                                        tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                <  1:5743920/93‑1 (MQ=255)
                                                                        tgATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                >  1:6413572/1‑93 (MQ=255)
                                                                         gATTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCg               >  2:7501788/1‑93 (MQ=255)
                                                                          aTTGGCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt              <  1:7024425/93‑1 (MQ=255)
                                                                             ggCAACGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAAc           <  1:6245906/93‑1 (MQ=255)
                                                                                aaCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc        >  1:6518360/1‑93 (MQ=255)
                                                                                aaCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc        >  2:2698895/1‑93 (MQ=255)
                                                                                aaCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc        >  1:5643335/1‑93 (MQ=255)
                                                                                 aCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCa       >  1:2868128/1‑93 (MQ=255)
                                                                                 aCGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCAGTTCGTTTGGTCTGCGTAACGTCa       >  2:7301368/1‑93 (MQ=255)
                                                                                  cGTTGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCAc      >  1:7746760/1‑93 (MQ=255)
                                                                                    ttGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg    <  2:5325722/93‑1 (MQ=255)
                                                                                    ttGAGTGTGTTCAGGCGAACGCTGAGGCGCTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCg    >  2:5896729/1‑93 (MQ=255)
                                                                                     tGAGTGTCTTCAGGCGAACGCTGAGGCACTACCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGa   <  2:2531232/93‑1 (MQ=255)
                                                                                          |                                                                                        
GCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCCCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGAC  >  NC_000913/4016707‑4016885

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: