Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 780,918 T→C intergenic (+43/‑390) lysQ → / → nadA tRNA‑Lys/quinolinate synthetase A protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913780,9180TC99.7% 749.5 / ‑5.9 303intergenic (+43/‑390)lysQ/nadAtRNA‑Lys/quinolinate synthetase A protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (177/125);  total (177/126)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 313 total aligned reads displayed.

TGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGC  >  NC_000913/780828‑780999
                                                                                          |                                                                                 
tGATTTTTAATCAATGGTTCGCAGGTCCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGacac                                                                                 <  1:3720542/93‑1 (MQ=11)
  aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTc                                                                               <  1:9059305/93‑1 (MQ=21)
  aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTc                                                                               >  1:132855/1‑93 (MQ=18)
  aCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTc                                                                               <  1:6197052/93‑1 (MQ=18)
     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            <  2:9598530/93‑1 (MQ=21)
     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            <  1:7666655/93‑1 (MQ=21)
     tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            >  2:6552132/1‑93 (MQ=18)
     ttaaaTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCaaa                                                                            <  2:360151/90‑1 (MQ=255)
       tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGa                                                                          >  1:1507085/1‑93 (MQ=18)
         aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGAtt                                                                        <  1:331651/93‑1 (MQ=21)
           cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCg                                                                      >  1:9950126/1‑93 (MQ=18)
            aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGa                                                                     >  2:6098484/1‑93 (MQ=18)
            aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGa                                                                     >  1:9177038/1‑93 (MQ=18)
             aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGac                                                                    >  2:1937203/1‑92 (MQ=18)
             aTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGaa                                                                    >  2:4623360/1‑93 (MQ=18)
              ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAc                                                                   <  1:6274197/93‑1 (MQ=18)
               tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAcc                                                                  <  1:7027451/93‑1 (MQ=18)
               tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAcc                                                                  <  2:4969545/93‑1 (MQ=21)
               tGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAAcc                                                                  <  2:8413080/93‑1 (MQ=21)
                ggTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCt                                                                 <  1:3126321/93‑1 (MQ=21)
                  gcGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgc                                                               <  1:4773227/92‑1 (MQ=21)
                   cGCAGTTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgca                                                              <  1:1842734/93‑1 (MQ=16)
                    gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGACCCTgcag                                                             <  1:4957152/93‑1 (MQ=16)
                     cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcagc                                                            <  1:2027820/93‑1 (MQ=21)
                     cAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTgcagc                                                            <  1:2378692/93‑1 (MQ=21)
                       ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGcagcag                                                          >  2:3499782/1‑93 (MQ=18)
                        gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         >  2:4915456/1‑93 (MQ=18)
                        gTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAgg                                                         >  1:8396799/1‑93 (MQ=18)
                         ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCCGGt                                                        >  2:8884915/1‑93 (MQ=16)
                         ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGt                                                        >  1:527388/1‑93 (MQ=18)
                          tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGtt                                                       <  2:8510940/93‑1 (MQ=18)
                          tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGtt                                                       <  1:611312/93‑1 (MQ=18)
                           cGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGttt                                                      >  1:6692460/1‑93 (MQ=18)
                            gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttg                                                     <  2:2946902/93‑1 (MQ=21)
                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttga                                                    <  1:7348349/93‑1 (MQ=21)
                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttga                                                    <  1:4051716/93‑1 (MQ=18)
                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttga                                                    >  1:4370776/1‑93 (MQ=18)
                             aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTCCAGCAGGTttga                                                    >  2:7738741/1‑93 (MQ=16)
                               tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                  <  1:914533/93‑1 (MQ=21)
                                ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagtt                                                 <  2:8665222/93‑1 (MQ=18)
                                   gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                              <  1:5508167/93‑1 (MQ=21)
                                   gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTttgagttgag                                              >  2:4853028/1‑93 (MQ=18)
                                    cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgc                                             >  2:1426317/1‑93 (MQ=18)
                                     aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcg                                            <  1:2722253/93‑1 (MQ=18)
                                      cGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                           >  1:6450169/1‑93 (MQ=18)
                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                          >  2:9368450/1‑93 (MQ=18)
                                       gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                          <  2:596833/93‑1 (MQ=18)
                                         cccACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                        >  1:8011432/1‑93 (MQ=18)
                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                       >  1:9564285/1‑93 (MQ=18)
                                          ccaccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                       >  1:5916301/1‑93 (MQ=18)
                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                      >  2:8604845/1‑93 (MQ=18)
                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGGTTGAGTTGAGCGCAGCGa                                      >  1:4012476/1‑93 (MQ=16)
                                           caccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGGTTGAGTTGAGCGCAGCGa                                      >  2:2755921/1‑93 (MQ=16)
                                             ccaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                    >  1:8752183/1‑93 (MQ=21)
                                             ccaATGTAAAAAAGCGCCCTAAAGCCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaaa                                    <  2:5512211/93‑1 (MQ=16)
                                                aTGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaa                                 <  2:9805440/93‑1 (MQ=18)
                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                >  2:2259360/1‑93 (MQ=18)
                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                >  1:10028703/1‑93 (MQ=18)
                                                 tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaacaac                                >  2:5780876/1‑93 (MQ=21)
                                                  gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                               >  2:9274439/1‑93 (MQ=21)
                                                  gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACg                               <  2:9513210/93‑1 (MQ=18)
                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              <  2:8046728/93‑1 (MQ=18)
                                                   tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                              >  2:10020888/1‑93 (MQ=21)
                                                    aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                             >  1:1023332/1‑93 (MQ=18)
                                                    aaaaaaGCGCCCTAAAGGCGCGTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGa                             <  1:7815235/93‑1 (MQ=16)
                                                     aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAg                            <  2:9135199/93‑1 (MQ=18)
                                                      aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                           <  2:1652760/93‑1 (MQ=18)
                                                       aaaGCGCCCTAAAGTCTCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCATCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                          <  1:862129/93‑1 (MQ=11)
                                                       aaaGCGCCATAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                          <  2:3761899/93‑1 (MQ=16)
                                                        aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCg                         >  1:8086726/1‑93 (MQ=21)
                                                         aGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGc                        <  1:9296311/93‑1 (MQ=18)
                                                           ctcccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCCGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTc                      <  2:3648106/91‑1 (MQ=255)
                                                             cccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgc                    <  1:5134430/93‑1 (MQ=18)
                                                              ccTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                   >  2:428919/1‑93 (MQ=18)
                                                               cTAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAAAAACGGAGCCGCTCgcgg                  >  2:4055619/1‑93 (MQ=255)
                                                               cTAAAGGCGCTTTTTTGCTATCTACGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcgg                  >  2:1117119/1‑93 (MQ=16)
                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCCGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 >  1:9718366/1‑93 (MQ=16)
                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 <  1:4639017/93‑1 (MQ=18)
                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 >  2:2404606/1‑93 (MQ=18)
                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 >  1:9269687/1‑93 (MQ=21)
                                                                tAAAGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggc                 >  2:412559/1‑93 (MQ=21)
                                                                 aaaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggcg                >  1:2538895/1‑93 (MQ=21)
                                                                 aaaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggcg                >  1:2942074/1‑93 (MQ=18)
                                                                 aaaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCgcggcg                >  2:2189580/1‑93 (MQ=18)
                                                                 aaaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCACCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCccggcc                >  2:7227549/1‑92 (MQ=11)
                                                                  aaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa               >  2:7007434/1‑93 (MQ=18)
                                                                  aaGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGa               >  2:3868328/1‑93 (MQ=18)
                                                                   aGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc              >  2:6190196/1‑93 (MQ=21)
                                                                   aGGCGCTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc              >  1:7176863/1‑93 (MQ=18)
                                                                       gcTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGcc          <  1:9859673/93‑1 (MQ=18)
                                                                        cTTTTTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAAGAccc         >  2:6482514/1‑93 (MQ=255)
                                                                        cTTTGTTGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGccc         <  2:4736480/93‑1 (MQ=16)
                                                                          tttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa       >  1:5166216/1‑93 (MQ=21)
                                                                          tttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGa       >  2:6752024/1‑93 (MQ=18)
                                                                           ttttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa      <  2:2361498/93‑1 (MQ=18)
                                                                           tgttGCTATCTGCGACACTCAAAGATTCGAATCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGaa      <  1:5711716/91‑1 (MQ=255)
                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     >  2:7477403/1‑93 (MQ=18)
                                                                            tttGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg     >  2:5277231/1‑93 (MQ=18)
                                                                              tGCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAggg   <  1:4264956/93‑1 (MQ=17)
                                                                               gCTATCTGCGACACTCAAAGATTCGAACCTGCAGCAGGCTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgc  >  1:4996890/1‑93 (MQ=9)
                                                                                          |                                                                                 
TGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGC  >  NC_000913/780828‑780999

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