Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 522,463 T→C V683A (GTG→GCG)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913522,4630TC98.7% 271.3 / ‑5.1 76V683A (GTG→GCG) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (35/40);  minor base G (0/1);  total (35/41)
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

TTTCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGT  >  NC_000913/522329‑522597
                                                                                                                                      |                                                                                                                                       
tttCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtg                                                                                                                                     <  2:1139110/139‑1 (MQ=255)
 ttCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtgg                                                                                                                                    <  2:317903/139‑1 (MQ=255)
  tCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtggt                                                                                                                                   >  1:1001903/1‑139 (MQ=255)
            gaATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTc                                                                                                                         <  2:278050/139‑1 (MQ=255)
                      gaaaagTTGACACAACTCAACCGGCAGGTCCCGACTATTAGCCTGTTAGGTATTGGGACAATTTTAAAACAGGTGGGTCAGGTGCTGAAGCAGATAAGTCTTGCGCTAGAAGGGATGGTGGTACTGGTCACCGCCTTCg                                                                                                              <  1:1649996/134‑1 (MQ=255)
                                     cTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtgg                                                                                                                                    <  1:1284571/103‑1 (MQ=255)
                                     cTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtgg                                                                                                                                    >  2:1284571/1‑103 (MQ=255)
                                     cTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtgttgctgttgc                                                                                               >  1:32388/1‑139 (MQ=255)
                                     cTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTAtgttgctgttgc                                                                                               >  1:496737/1‑139 (MQ=255)
                                      tCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATgttgctgttgct                                                                                              <  2:1001903/139‑1 (MQ=255)
                                        aaCCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAt                                                                                                                                      <  2:261660/98‑1 (MQ=255)
                                        aaCCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAt                                                                                                                                      >  1:261660/1‑98 (MQ=255)
                                         aCCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGc                                                                                           <  2:1342378/139‑1 (MQ=255)
                                          ccGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                          >  1:105790/1‑139 (MQ=255)
                                            gCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGcaca                                                                                        <  1:286157/139‑1 (MQ=255)
                                                gTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                    <  2:32388/139‑1 (MQ=255)
                                                 ttCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgc                                                                                   <  1:761638/139‑1 (MQ=255)
                                                          aTTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGaa                                                                          <  2:496737/139‑1 (MQ=255)
                                                           ttAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                   >  1:474679/1‑97 (MQ=255)
                                                           ttAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                   <  2:474679/97‑1 (MQ=255)
                                                            tAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATg                                                                        >  1:96102/1‑139 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTgg                                                                                            >  2:1866909/1‑102 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTgg                                                                                            <  1:1866909/102‑1 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGcac                                                                                         <  1:1236006/105‑1 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGcac                                                                                         >  2:1236006/1‑105 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtc                                                                    >  2:683625/1‑126 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtc                                                                    <  1:683625/126‑1 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGa                                                       >  1:521108/1‑139 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGa                                                       >  2:772667/1‑139 (MQ=255)
                                                                             gcgcGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGa                                                       >  2:541784/1‑139 (MQ=255)
                                                                                cgATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgca                                                                                  <  2:1804452/109‑1 (MQ=255)
                                                                                cgATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtgca                                                                                  >  1:1804452/1‑109 (MQ=255)
                                                                                cgATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCt                                                    >  2:1452957/1‑139 (MQ=255)
                                                                                 gATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtgg                                                                                                                                    <  1:374489/59‑1 (MQ=255)
                                                                                 gATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGAtgg                                                                                                                                    >  2:374489/1‑59 (MQ=255)
                                                                                  aTTTTAAAACAGGTCGGTCAGGCGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCtg                                                   >  1:784459/1‑139 (MQ=255)
                                                                                     ttAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCtggtg                                                <  2:784459/139‑1 (MQ=255)
                                                                                      tAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTggtggt                                              <  2:1432260/139‑1 (MQ=255)
                                                                                       aaaaCAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                             <  1:772667/139‑1 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                   >  1:754850/1‑63 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                   <  2:754850/63‑1 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                          <  2:21234/88‑1 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                          >  1:21234/1‑88 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                    >  1:834593/1‑94 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggtg                                                                                    <  2:834593/94‑1 (MQ=255)
                                                                                             ggtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGgcgc                                       >  1:622971/1‑139 (MQ=255)
                                                                                              gtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGca                                      >  2:515392/1‑139 (MQ=255)
                                                                                              gtcggtcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGca                                      >  1:8970/1‑139 (MQ=255)
                                                                                                   tcAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTc                                 >  2:972353/1‑139 (MQ=255)
                                                                                                      ggTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGt                              <  2:105790/139‑1 (MQ=255)
                                                                                                        tGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                   <  2:96809/52‑1 (MQ=255)
                                                                                                        tGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGc                                                                                                                   >  1:96809/1‑52 (MQ=255)
                                                                                                        tGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCt                                                                                                                 <  2:673793/54‑1 (MQ=255)
                                                                                                        tGCTGGAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCt                                                                                                                 >  1:673793/1‑54 (MQ=255)
                                                                                                             gaggagGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggt                                                                                     >  1:1208604/1‑77 (MQ=255)
                                                                                                             gAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCAcaggt                                                                                     <  2:1208604/77‑1 (MQ=255)
                                                                                                             gAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                             <  2:1306336/117‑1 (MQ=255)
                                                                                                             gAGCAGGTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGgtggtg                                             >  1:1306336/1‑117 (MQ=255)
                                                                                                                  ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagcgtc                                                              >  2:1632013/1‑95 (MQ=255)
                                                                                                                  ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATgcgtcagcgtc                                                              <  1:1632013/95‑1 (MQ=255)
                                                                                                                  ggTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGaaaaaa                  <  2:1756962/139‑1 (MQ=255)
                                                                                                                   gTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAg                                                         <  1:1814775/99‑1 (MQ=255)
                                                                                                                   gTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAg                                                         >  2:1814775/1‑99 (MQ=255)
                                                                                                                   gTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAAc                 >  1:1286917/1‑139 (MQ=255)
                                                                                                                   gTAAGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAAc                 >  2:1622826/1‑139 (MQ=255)
                                                                                                                     aaGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAActg               <  2:96102/139‑1 (MQ=255)
                                                                                                                      aGTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAActgc              <  1:1622826/139‑1 (MQ=255)
                                                                                                                       gTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAActgct             <  1:541784/139‑1 (MQ=255)
                                                                                                                       gTCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAActgct             <  2:1286917/139‑1 (MQ=255)
                                                                                                                        tCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                          >  1:287726/1‑61 (MQ=255)
                                                                                                                        tCGGGCGCTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGca                                                                                          <  2:287726/61‑1 (MQ=255)
                                                                                                                            gcgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTa        <  1:373362/139‑1 (MQ=255)
                                                                                                                            gcgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTa        <  1:153089/139‑1 (MQ=255)
                                                                                                                             cgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTc                                 >  1:967885/1‑113 (MQ=255)
                                                                                                                             cgcTGGAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTc                                 <  2:967885/113‑1 (MQ=255)
                                                                                                                                  gAAGCGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGt  <  2:521108/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
TTTCCGCTGGGAGAATGGCAACGGCATGTTGACACAACTCAACCGCCAGTTCCCGACCATTAGCCTGTTAGATATTGGCGCGATTTTAAAACAGGTCGGTCAGGTGCTGGAGCAGGTAAGTCGGGCGCTGGAAGTGATGGTGGTACTGGTCACCGCCTGCGGTATGTTGCTGTTGCTGGCACAGGTGCAGGTGGGAATGCGTCAGCGTCATCAGGAGCTGGTGGTGTGGCGCACACTCGGTGCGGGGAAAAAACTGCTGCGTACCACGT  >  NC_000913/522329‑522597

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