Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,284,989 (C)5→4 coding (1932/2511 nt) yejO ← putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,284,9850C.99.2% 586.6 / ‑3.1 130coding (1936/2511 nt)yejOputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (63/66);  total (64/66)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 146 total aligned reads displayed.

ACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATT  >  NC_000913/2284893‑2285073
                                                                                            |                                                                                        
aCTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGc                                                                                          >  1:4293424/1‑93 (MQ=255)
 cTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGcc                                                                                         >  1:719726/1‑93 (MQ=255)
  tCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACCTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccc                                                                                        <  1:523845/93‑1 (MQ=255)
    aTCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCACCCGGCATAAGcccca                                                                                      >  1:3386344/1‑92 (MQ=255)
    aTCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGcccca                                                                                      <  1:325560/93‑2 (MQ=255)
     tCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAGCCGGCATAAGccccag                                                                                     <  2:748594/93‑3 (MQ=255)
     tCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCAACCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATACGcccccg                                                                                     >  2:3862408/1‑92 (MQ=255)
      cTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccagg                                                                                    <  1:2041467/93‑4 (MQ=255)
      cTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccagg                                                                                    <  1:2103754/93‑4 (MQ=255)
      cTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccagg                                                                                    <  2:180951/93‑4 (MQ=255)
       ttGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccaggg                                                                                   >  1:698479/1‑89 (MQ=255)
        tGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGt                                                                                 <  2:3904253/93‑1 (MQ=255)
        tGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGt                                                                                 >  1:2938168/1‑93 (MQ=255)
        tGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGt                                                                                 <  1:987309/93‑1 (MQ=255)
          ccATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGtat                                                                               <  1:1613623/93‑1 (MQ=255)
           cATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGtat                                                                               >  2:4531524/1‑93 (MQ=255)
           cATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGtata                                                                              <  1:3613752/93‑1 (MQ=255)
           cATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGtata                                                                              >  2:2172470/1‑93 (MQ=255)
            atggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAg                                                                             >  1:857339/1‑93 (MQ=255)
            atggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAg                                                                             >  2:2615584/1‑93 (MQ=255)
            atggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAg                                                                             <  1:665186/93‑1 (MQ=255)
             tggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGc                                                                            <  2:3937727/93‑1 (MQ=255)
             tggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCAGCATAAGCCCCAGGGTATAGc                                                                            >  1:4265739/1‑93 (MQ=255)
              ggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCt                                                                           >  1:2759526/1‑93 (MQ=255)
              ggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGAATAGca                                                                           >  2:3349150/1‑92 (MQ=255)
                atggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTAt                                                                         >  1:3495770/1‑93 (MQ=255)
                  ggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGCATAGCTATCg                                                                       >  1:2409357/1‑93 (MQ=255)
                   gtgtTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGa                                                                      <  1:33312/93‑1 (MQ=255)
                    tgtTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGat                                                                     <  1:914719/93‑1 (MQ=255)
                    tgtTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGat                                                                     >  2:2287/1‑93 (MQ=255)
                     gtTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGata                                                                    <  1:1497488/93‑1 (MQ=255)
                       tGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATAtt                                                                  <  2:1487135/93‑1 (MQ=255)
                       tGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATAtt                                                                  >  2:4255307/1‑93 (MQ=255)
                        ggCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTa                                                                 <  1:620161/93‑1 (MQ=255)
                          cAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATTCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATGGATATTAcc                                                               <  1:3302556/93‑1 (MQ=255)
                           aaaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCt                                                              <  2:174502/93‑1 (MQ=255)
                           aaaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCt                                                              >  2:2293301/1‑93 (MQ=255)
                            aaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCtt                                                             <  2:1845186/93‑1 (MQ=255)
                            aaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCtt                                                             >  1:2764880/1‑93 (MQ=255)
                             aaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCttt                                                            >  1:4026995/1‑93 (MQ=255)
                              aCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTg                                                           >  1:2690533/1‑93 (MQ=255)
                              aCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTg                                                           >  1:4748211/1‑93 (MQ=255)
                               cGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTgc                                                          <  1:4506668/93‑1 (MQ=255)
                                cgTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTgcc                                                         <  1:2148988/92‑1 (MQ=255)
                                 gTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTgccg                                                        >  1:470925/1‑93 (MQ=255)
                                    aaCTTTCACCACCCCATCAACGTAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGccgccg                                                     >  2:2020220/1‑93 (MQ=255)
                                     aCTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATGGATATTACCTTTCcccccgg                                                    >  1:1348987/1‑92 (MQ=255)
                                     aCTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCcgccgc                                                    >  2:4435446/1‑93 (MQ=255)
                                         tCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATc                                                <  1:1560379/93‑1 (MQ=255)
                                           accaccCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCaa                                              >  1:651590/1‑93 (MQ=255)
                                           accaccCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCaa                                              <  1:2272337/93‑1 (MQ=255)
                                           accaccCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCaa                                              >  2:3754813/1‑93 (MQ=255)
                                              cccccATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAc                                           <  2:3520897/92‑1 (MQ=255)
                                              accCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAc                                           >  1:3858409/1‑93 (MQ=255)
                                              accCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAc                                           >  1:2371032/1‑93 (MQ=255)
                                                cccATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCa                                         >  2:2335294/1‑93 (MQ=255)
                                                cccATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCa                                         <  1:4149339/93‑1 (MQ=255)
                                                 ccATCAACACAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCaa                                        <  1:267614/93‑1 (MQ=255)
                                                  cATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAat                                       <  2:3234863/93‑1 (MQ=255)
                                                  cATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAat                                       >  2:2826362/1‑93 (MQ=255)
                                                    tCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAatat                                     >  2:1065687/1‑93 (MQ=255)
                                                    tCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAatat                                     >  2:2041436/1‑93 (MQ=255)
                                                     cAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATc                                    <  1:3081612/93‑1 (MQ=255)
                                                      aaCATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCa                                   <  1:757555/93‑1 (MQ=255)
                                                      aaCATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCa                                   <  2:4681328/93‑1 (MQ=255)
                                                        cATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGa                                 <  1:3703425/93‑1 (MQ=255)
                                                         aTAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGaa                                >  2:615644/1‑93 (MQ=255)
                                                         aTAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGaa                                <  1:1643425/93‑1 (MQ=255)
                                                          tAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAAt                               >  2:231708/1‑93 (MQ=255)
                                                            ggCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATga                             <  2:2213826/93‑1 (MQ=255)
                                                            ggCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATga                             <  1:2049896/93‑1 (MQ=255)
                                                            ggCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATga                             >  2:170180/1‑93 (MQ=255)
                                                              cACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATgaga                           <  1:110530/93‑1 (MQ=255)
                                                              cACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATgaga                           <  1:1057428/93‑1 (MQ=255)
                                                               aCCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagag                          <  2:3852323/93‑1 (MQ=255)
                                                               aCCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagag                          >  2:3246282/1‑93 (MQ=255)
                                                                 cGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTa                        <  1:3694008/93‑1 (MQ=255)
                                                                 cGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTa                        <  2:4624731/93‑1 (MQ=255)
                                                                   ttCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAAc                      >  1:4227508/1‑93 (MQ=255)
                                                                   ttCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAAc                      <  1:3912152/93‑1 (MQ=255)
                                                                   ttCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAAc                      >  2:613393/1‑93 (MQ=255)
                                                                     cTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCa                    >  1:1339625/1‑93 (MQ=255)
                                                                      tGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCaa                   >  1:45819/1‑93 (MQ=255)
                                                                       gATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCaaa                  <  1:4073240/93‑1 (MQ=255)
                                                                         tGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGa                <  2:350366/93‑1 (MQ=255)
                                                                          gCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGAt               <  1:1217860/93‑1 (MQ=255)
                                                                          gCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGAt               <  2:3605444/93‑1 (MQ=255)
                                                                          gCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGAt               <  1:4130350/93‑1 (MQ=255)
                                                                           cTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATc              <  1:1186864/93‑1 (MQ=255)
                                                                             cccAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCaa            <  1:9529/93‑1 (MQ=255)
                                                                                aaCCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGcc         <  2:956804/93‑1 (MQ=255)
                                                                                aaCCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGcc         >  2:3824711/1‑93 (MQ=255)
                                                                                aaCCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGcc         >  2:1124334/1‑93 (MQ=255)
                                                                                aaCCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGcc         <  2:1506515/93‑1 (MQ=255)
                                                                                  ccGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCcgc       <  2:4527960/93‑1 (MQ=255)
                                                                                    ggCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGACGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGACGCGa     >  2:1753042/1‑93 (MQ=255)
                                                                                     gCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGaa    >  2:3379640/1‑93 (MQ=255)
                                                                                       aTAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAAtt  >  1:138392/1‑93 (MQ=255)
                                                                                       aTAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAAtt  >  1:2885908/1‑93 (MQ=255)
                                                                                       aTAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAAtt  <  2:2573214/93‑1 (MQ=255)
                                                                                            |                                                                                        
ACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATT  >  NC_000913/2284893‑2285073

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