Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,926 T→C T3A (ACC→GCC)  b1297 ← putative glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,9260TC99.5% 609.6 / NA 184T3A (ACC→GCC) b1297putative glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (143/40);  minor base A (1/0);  total (144/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

GCGCTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGTTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTA  >  NC_000913/1358834‑1359014
                                                                                            |                                                                                        
gcgcTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGt                                                                                          >  2:1601154/1‑93 (MQ=255)
   cTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTc                                                                                       >  1:2807604/1‑93 (MQ=255)
    tCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTTCTTCAACGATATTGGCTTct                                                                                      >  1:3142030/1‑92 (MQ=255)
    tCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTcc                                                                                      >  2:3669388/1‑93 (MQ=255)
    cgcgggTCGCGTTTCGCTCCTTGCAAAAGGGCCCTCCCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTcc                                                                                      <  2:295846/89‑1 (MQ=255)
      cAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCAt                                                                                    >  2:4191328/1‑93 (MQ=255)
       aGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATg                                                                                   >  2:2926653/1‑93 (MQ=255)
        ggTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGa                                                                                  >  1:384671/1‑93 (MQ=255)
        ggTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGa                                                                                  >  1:208599/1‑93 (MQ=255)
          tAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGAtt                                                                                >  2:4053217/1‑93 (MQ=255)
           aGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGAttc                                                                               >  2:3841364/1‑93 (MQ=255)
            gCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGAttct                                                                              >  2:2427663/1‑93 (MQ=255)
                ttCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGc                                                                          >  2:2165177/1‑93 (MQ=255)
                 tCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGcc                                                                         >  1:1126834/1‑93 (MQ=255)
                  cACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCt                                                                        >  1:2559077/1‑93 (MQ=255)
                   aCTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCtt                                                                       >  2:3739379/1‑93 (MQ=255)
                   aCTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCtt                                                                       >  2:1687926/1‑93 (MQ=255)
                    cTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCttt                                                                      >  2:1055775/1‑93 (MQ=255)
                     ttCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTg                                                                     >  2:4557649/1‑93 (MQ=255)
                     ttCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTg                                                                     >  1:4670347/1‑93 (MQ=255)
                      tCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTgg                                                                    >  1:3124187/1‑93 (MQ=255)
                        cTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtt                                                                  >  2:2890125/1‑93 (MQ=255)
                         tGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgttt                                                                 >  1:3046715/1‑93 (MQ=255)
                         tGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgttt                                                                 >  1:485543/1‑93 (MQ=255)
                          gCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttg                                                                >  1:4614739/1‑93 (MQ=255)
                            gTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttgtt                                                              >  1:4122043/1‑93 (MQ=255)
                              aaaGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGtttt                                                            >  2:1188202/1‑93 (MQ=255)
                                    gcTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTc                                                      >  2:3042544/1‑93 (MQ=255)
                                    cctcccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTc                                                      <  1:2954080/89‑1 (MQ=255)
                                     cTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCg                                                     >  1:3858496/1‑93 (MQ=255)
                                     cTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCg                                                     >  1:1397715/1‑93 (MQ=255)
                                     cTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGCTTGTTTTGCGCttc                                                     >  2:1848727/1‑91 (MQ=255)
                                      tACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGt                                                    >  2:214489/1‑93 (MQ=255)
                                      ccccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGt                                                    <  2:3754840/91‑1 (MQ=255)
                                       cccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGtt                                                   <  2:4151375/92‑1 (MQ=255)
                                       aCCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGtt                                                   >  2:3486120/1‑93 (MQ=255)
                                        ccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGttt                                                  >  1:1624557/1‑92 (MQ=255)
                                        ccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTa                                                  >  1:4736356/1‑93 (MQ=255)
                                         cccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTAt                                                 <  2:1962123/93‑1 (MQ=255)
                                         cccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTAt                                                 <  2:4530616/93‑1 (MQ=255)
                                         cccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTAt                                                 <  2:162734/93‑1 (MQ=255)
                                         cccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTAt                                                 >  1:2714544/1‑93 (MQ=255)
                                          ccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATc                                                >  1:1820692/1‑93 (MQ=255)
                                           cGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCa                                               >  2:2924881/1‑93 (MQ=255)
                                            gCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCaa                                              >  2:2035192/1‑93 (MQ=255)
                                            gCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCaa                                              <  2:893143/93‑1 (MQ=255)
                                            gCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCaa                                              >  1:196353/1‑93 (MQ=255)
                                               tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGc                                           >  2:4482915/1‑93 (MQ=255)
                                               tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGc                                           <  1:2888294/93‑1 (MQ=255)
                                                 tcttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGt                                         >  2:4276536/1‑93 (MQ=255)
                                                 tcttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGt                                         >  2:3206996/1‑93 (MQ=255)
                                                   ttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTaa                                       <  1:3467543/93‑1 (MQ=255)
                                                   ttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTaa                                       >  1:3503398/1‑93 (MQ=255)
                                                    tctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTaaa                                      >  2:2736113/1‑93 (MQ=255)
                                                      tGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAat                                    >  1:3383456/1‑93 (MQ=255)
                                                       gACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAata                                   >  2:597207/1‑93 (MQ=255)
                                                       gACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAata                                   <  1:2034477/93‑1 (MQ=255)
                                                       gACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAata                                   >  1:339106/1‑93 (MQ=255)
                                                           gctgAACAAAGTTCTCTACTTCAACGATATTGGCTTTCATGATTATTCGCCTTTGGGTTGCTTTACCCTCGTTATCAAAGCGTAAACTATAAt                               >  2:786201/1‑93 (MQ=255)
                                                           gctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAAt                               <  1:940279/93‑1 (MQ=255)
                                                           gctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAAt                               >  1:1051020/1‑93 (MQ=255)
                                                           gctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAAt                               >  1:3757398/1‑93 (MQ=255)
                                                           gctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAAt                               <  1:4008697/93‑1 (MQ=255)
                                                             tgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGa                             >  2:524973/1‑93 (MQ=255)
                                                             tgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGa                             >  2:4491274/1‑93 (MQ=255)
                                                                aCAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGacca                          <  2:609980/93‑1 (MQ=255)
                                                                aCAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGacca                          >  2:3469949/1‑93 (MQ=255)
                                                                 cAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGaccac                         >  1:4248274/1‑93 (MQ=255)
                                                                 cAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGCAAAATATAATGaccac                         >  1:3523217/1‑93 (MQ=255)
                                                                  aaaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGaccacc                        >  1:3951309/1‑93 (MQ=255)
                                                                  aaaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGaccacc                        >  2:2665061/1‑93 (MQ=255)
                                                                  aaaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGaccacc                        >  2:1676758/1‑93 (MQ=255)
                                                                  aaaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGaccacc                        <  1:3177399/93‑1 (MQ=255)
                                                                   aaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGAccacca                       <  1:883257/93‑1 (MQ=255)
                                                                   aaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGAccacca                       >  1:1235016/1‑93 (MQ=255)
                                                                    aGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCAt                      <  2:4205629/93‑1 (MQ=255)
                                                                      ttCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTc                    >  1:2447749/1‑93 (MQ=255)
                                                                      ttCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTc                    >  2:1474742/1‑93 (MQ=255)
                                                                       tctcTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCg                   <  1:467452/93‑1 (MQ=255)
                                                                        ctctACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGa                  >  1:827377/1‑93 (MQ=255)
                                                                        ctctACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGa                  >  2:437154/1‑93 (MQ=255)
                                                                          ctACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAAt                >  1:432185/1‑93 (MQ=255)
                                                                              ttCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGt            >  1:1662299/1‑93 (MQ=255)
                                                                              ttCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGt            >  2:1290549/1‑93 (MQ=255)
                                                                              ttCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGt            >  2:3750938/1‑93 (MQ=255)
                                                                              ttCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGt            >  1:3328413/1‑93 (MQ=255)
                                                                                cAACGATATTGGCTTCCATGATGCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAAGATAATGACCACCATTCGAATCTGTAt          <  2:3755803/93‑1 (MQ=255)
                                                                                 aaCTATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg         >  2:1808328/1‑93 (MQ=255)
                                                                                 aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg         >  2:4010253/1‑93 (MQ=255)
                                                                                 aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg         >  1:4494547/1‑93 (MQ=255)
                                                                                 aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg         >  2:1548162/1‑93 (MQ=255)
                                                                                 aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg         <  2:1871561/93‑1 (MQ=255)
                                                                                  aCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGc        >  2:1399370/1‑93 (MQ=255)
                                                                                  aCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGc        >  1:2204263/1‑93 (MQ=255)
                                                                                  aCGATATTGGCTTCCATGATTCTTCGCCTTTAGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATGTAATGACCACCATTCGAATCTGTATGc        >  1:3710784/1‑93 (MQ=255)
                                                                                   cGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCa       >  1:506575/1‑93 (MQ=255)
                                                                                      tatTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAAc    <  1:1869312/93‑1 (MQ=255)
                                                                                      tatTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAAc    >  2:1310739/1‑93 (MQ=255)
                                                                                      tatTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAAc    >  1:1358991/1‑93 (MQ=255)
                                                                                        ttGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTa  >  2:335461/1‑93 (MQ=255)
                                                                                            |                                                                                        
GCGCTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGTTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTA  >  NC_000913/1358834‑1359014

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