Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,033,600 A→G intergenic (+335/‑257) b1966 → / → yedU putative outer membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,033,6000AG83.3% 583.6 / 101.7 263intergenic (+335/‑257)b1966/yedUputative outer membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/14);  new base G (129/90);  total (159/104)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 264 total aligned reads displayed.

TAGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGC  >  NC_000913/2033509‑2033692
                                                                                           |                                                                                            
  gAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTa                                                                                           >  1:3916693/1‑93 (MQ=255)
   aaaTATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTaa                                                                                          <  2:2992243/93‑1 (MQ=255)
   aaaTATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTaa                                                                                          >  1:2668536/1‑93 (MQ=255)
   aaaTATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTaa                                                                                          >  1:1653755/1‑93 (MQ=255)
    aaTATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTaaa                                                                                         >  1:94564/1‑93 (MQ=255)
    aaTATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTaaa                                                                                         >  1:1788114/1‑93 (MQ=255)
      tatTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtc                                                                                       <  1:7104411/93‑1 (MQ=255)
      tatTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtc                                                                                       <  1:6507535/93‑1 (MQ=255)
      tatTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtc                                                                                       >  2:3243773/1‑93 (MQ=255)
      tatTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtc                                                                                       <  2:6602383/93‑1 (MQ=255)
      tatTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtc                                                                                       >  1:2863089/1‑93 (MQ=255)
       atTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtct                                                                                      <  1:913298/93‑1 (MQ=255)
        tttGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGGTAAAtctc                                                                                     <  1:6956543/93‑1 (MQ=255)
        tttGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtctc                                                                                     <  2:5531938/93‑1 (MQ=255)
        tttGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtctc                                                                                     >  1:2548000/1‑93 (MQ=255)
        tttGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtctc                                                                                     <  2:639293/93‑1 (MQ=255)
          tGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAAtctctc                                                                                   >  1:5401717/1‑93 (MQ=255)
           gCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTcc                                                                                  <  1:2660345/93‑1 (MQ=255)
           gCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTcc                                                                                  <  2:2754568/93‑1 (MQ=255)
           gCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTcc                                                                                  <  2:1674602/93‑1 (MQ=255)
            cAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTccc                                                                                 <  2:7453381/93‑1 (MQ=255)
            cAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTccc                                                                                 <  2:4014238/93‑1 (MQ=255)
             aaaaaTAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTcccc                                                                                >  2:3772163/1‑93 (MQ=255)
             aaaaaTAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTcccc                                                                                <  1:7464032/93‑1 (MQ=255)
               aaaTAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGc                                                                              >  2:7149570/1‑93 (MQ=255)
                aaTAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGcc                                                                             >  1:3636406/1‑93 (MQ=255)
                    gTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCAccc                                                                         >  1:1733768/1‑93 (MQ=255)
                      gCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCg                                                                       <  2:1748133/93‑1 (MQ=255)
                        aTGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCGGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTa                                                                     <  1:1936442/93‑1 (MQ=255)
                        aTGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTa                                                                     >  1:7414153/1‑93 (MQ=255)
                         tGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACTCCGTAc                                                                    >  1:5004574/1‑93 (MQ=255)
                         tGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTAc                                                                    <  1:7210468/93‑1 (MQ=255)
                         tGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTAc                                                                    >  1:1258299/1‑93 (MQ=255)
                          gAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTAcc                                                                   >  2:15140/1‑93 (MQ=255)
                            aaTAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACctc                                                                 >  2:7138307/1‑93 (MQ=255)
                              tAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg                                                               >  2:4236426/1‑93 (MQ=255)
                              tAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTg                                                               >  2:673364/1‑93 (MQ=255)
                               aaGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGa                                                              >  2:5385801/1‑93 (MQ=255)
                                 gCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATa                                                            >  2:3533842/1‑93 (MQ=255)
                                    ttttCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATg                                                         >  2:2051386/1‑93 (MQ=255)
                                     tttCTCTACCAATCCGGGATATAATTATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATgg                                                        >  2:62079/1‑93 (MQ=255)
                                       tctcTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTc                                                      >  2:6967094/1‑93 (MQ=255)
                                       tctcTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTc                                                      >  1:529412/1‑93 (MQ=255)
                                        ctctACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCt                                                     >  2:1797591/1‑93 (MQ=255)
                                          ctACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTaa                                                   >  1:5811912/1‑93 (MQ=255)
                                          ctACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTaa                                                   <  2:3608978/93‑1 (MQ=255)
                                            aCCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCGCCCCGTACCTCTGATAATGGTCTaaaa                                                 <  2:3485554/93‑1 (MQ=255)
                                            aCCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTaaaa                                                 >  2:896742/1‑93 (MQ=255)
                                               aaTCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCa                                              <  2:6845073/93‑1 (MQ=255)
                                                  ccGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTg                                           <  1:791543/93‑1 (MQ=255)
                                                  ccGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTg                                           <  1:4581882/93‑1 (MQ=255)
                                                   cGGGATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGa                                          >  1:5146985/1‑93 (MQ=255)
                                                    gggATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGaa                                         <  1:1810265/93‑1 (MQ=255)
                                                    gggATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGaa                                         >  1:4873895/1‑93 (MQ=255)
                                                     ggATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAg                                        <  2:247822/93‑1 (MQ=255)
                                                      gATATAAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGc                                       <  2:5720008/93‑1 (MQ=255)
                                                        tataAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCa                                     >  1:1143423/1‑93 (MQ=255)
                                                        tataAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCa                                     >  2:4780778/1‑93 (MQ=255)
                                                         ataAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCAc                                    >  2:2544722/1‑93 (MQ=255)
                                                         ataAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCAc                                    >  2:6514670/1‑93 (MQ=255)
                                                          taAATATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACt                                   >  1:2031682/1‑93 (MQ=255)
                                                          taAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACt                                   >  2:151790/1‑93 (MQ=255)
                                                           aaaTATTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACtt                                  <  2:5231014/93‑1 (MQ=255)
                                                              tatTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcg                               >  2:4702861/1‑93 (MQ=255)
                                                              tatTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcg                               >  1:593468/1‑93 (MQ=255)
                                                              tatTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcg                               >  2:1081011/1‑93 (MQ=255)
                                                              tatTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcg                               >  1:7184892/1‑93 (MQ=255)
                                                               atTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcga                              <  1:3466534/93‑1 (MQ=255)
                                                               atTCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcga                              <  2:5261022/93‑1 (MQ=255)
                                                                ttCATCAGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgac                             >  2:4454909/1‑93 (MQ=255)
                                                                ttCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgac                             <  2:5906440/93‑1 (MQ=255)
                                                                  catcaGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                           <  2:3094698/93‑1 (MQ=255)
                                                                  catcaGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                           <  1:1388928/93‑1 (MQ=255)
                                                                  catcaGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                           <  2:5048133/93‑1 (MQ=255)
                                                                  catcaGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGcgacga                           >  1:2157255/1‑93 (MQ=255)
                                                                    tcaGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGacgacg                         >  1:315408/1‑93 (MQ=255)
                                                                    tcaGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGacgacg                         >  1:2258342/1‑93 (MQ=255)
                                                                     caGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGt                        <  2:2096055/93‑1 (MQ=255)
                                                                      aGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGtt                       <  1:4625025/93‑1 (MQ=255)
                                                                      aGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGtt                       >  1:6567453/1‑93 (MQ=255)
                                                                      aGCACCACAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGtt                       <  2:7389027/93‑1 (MQ=255)
                                                                         accacAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCt                    >  1:7064414/1‑93 (MQ=255)
                                                                          ccacAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTc                   >  1:1425576/1‑93 (MQ=255)
                                                                          ccacAATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTc                   >  2:1281850/1‑93 (MQ=255)
                                                                          ccacAATATTCCACAGCACTAAATCTCCCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTc                   <  2:5668087/93‑1 (MQ=255)
                                                                            acaATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCaa                 >  1:1045872/1‑93 (MQ=255)
                                                                             caATATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAAt                >  2:7482980/1‑93 (MQ=255)
                                                                              aaTATTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATa               >  2:1588306/1‑93 (MQ=255)
                                                                               atatTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATaa              <  1:3989019/93‑1 (MQ=255)
                                                                               atatTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATaa              >  1:124704/1‑93 (MQ=255)
                                                                               atatTCCACAGCACTAAATCCCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATaa              <  1:4364171/93‑1 (MQ=255)
                                                                                tatTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAg             <  2:6292394/93‑1 (MQ=255)
                                                                                tatTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAg             >  1:1556538/1‑93 (MQ=255)
                                                                                 atTCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAgg            >  2:3871109/1‑93 (MQ=255)
                                                                                  ttCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGa           >  2:7028414/1‑93 (MQ=255)
                                                                                  ttCCACAGCACTAAATCCCCCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGa           <  1:6293701/93‑1 (MQ=255)
                                                                                   tCCACAGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGaa          <  1:6690730/93‑1 (MQ=255)
                                                                                   tCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGaa          <  2:6132831/93‑1 (MQ=255)
                                                                                      acaGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAAt       >  2:4389242/1‑93 (MQ=255)
                                                                                      acaGCGCTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAAt       >  2:3290844/1‑93 (MQ=255)
                                                                                           |                                                                                            
TAGAAATATTTGCAAAAATAGTGCATGAAATAAGCTTTTTCTCTACCAATCCGGGATATAAATATTCATCAGCACCACAATATTCCACAGCACTAAATCTCTCCCCGCCACCCCGTACCTCTGATAATGGTCTAAAATCATTGAAGCCACTTGCGACGACGTTGCTCAATAAGGAAAATTTGGC  >  NC_000913/2033509‑2033692

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