Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,564,703 A→G R61G (AGA→GGA)  yffS → CPZ‑55 prophage; uncharacterized protein YffS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,7030AG99.6% 1069.0 / NA 259R61G (AGA→GGA) yffSCPZ‑55 prophage; uncharacterized protein YffS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (108/150);  minor base T (0/1);  total (108/151)
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 259 total aligned reads displayed.

AATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATGGATTCGAGTTCA  >  NC_000913/2564605‑2564800
                                                                                                  |                                                                                                 
 aTGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGAc                                                                                                <  2:3890935/101‑1 (MQ=255)
  tGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACg                                                                                               >  2:3891499/1‑101 (MQ=255)
    cGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACggg                                                                                             <  1:3890937/101‑1 (MQ=255)
     gaggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGa                                                                                            <  1:3890904/101‑1 (MQ=255)
     gaggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGa                                                                                            >  2:3891504/1‑101 (MQ=255)
      aggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAg                                                                                           <  2:3891325/101‑1 (MQ=255)
       ggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGc                                                                                          <  1:3890870/101‑1 (MQ=255)
       ggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGc                                                                                          <  2:3890827/101‑1 (MQ=255)
       ggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGc                                                                                          <  1:3891004/101‑1 (MQ=255)
         agaTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCa                                                                                        >  1:3891509/1‑101 (MQ=255)
          gaTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCaa                                                                                       >  2:3891510/1‑101 (MQ=255)
             cGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAAtt                                                                                    <  2:3891383/101‑1 (MQ=255)
             cGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAAtt                                                                                    <  1:3891281/101‑1 (MQ=255)
               cAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTc                                                                                  <  2:3891355/101‑1 (MQ=255)
                aaaCGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCg                                                                                 <  1:3891167/101‑1 (MQ=255)
                aaaCGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCg                                                                                 >  1:3891513/1‑101 (MQ=255)
                aaaCGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCg                                                                                 <  2:3891106/101‑1 (MQ=255)
                  aCGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGac                                                                               >  1:3891516/1‑101 (MQ=255)
                   cGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGaca                                                                              <  1:3891431/101‑1 (MQ=255)
                   cGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGaca                                                                              <  2:3891100/101‑1 (MQ=255)
                   cGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGaca                                                                              >  2:3891518/1‑101 (MQ=255)
                   cGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGaca                                                                              >  2:3891520/1‑101 (MQ=255)
                     aGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACAct                                                                            <  1:3891179/101‑1 (MQ=255)
                     aGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACAct                                                                            <  2:3891291/101‑1 (MQ=255)
                      gTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACActc                                                                           >  1:3891524/1‑101 (MQ=255)
                       ttAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACActct                                                                          <  1:3891150/101‑1 (MQ=255)
                        tAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTg                                                                         <  2:3891131/101‑1 (MQ=255)
                        tAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTg                                                                         <  1:3891062/101‑1 (MQ=255)
                          aCTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGc                                                                       <  1:3891248/101‑1 (MQ=255)
                          aCTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGc                                                                       >  1:3891527/1‑101 (MQ=255)
                           cTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCt                                                                      <  2:3891253/101‑1 (MQ=255)
                             cGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGa                                                                    >  2:3891530/1‑101 (MQ=255)
                              gCCGGGGAGTTAAAAACAACAAGGGGGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGaa                                                                   >  1:3891532/1‑101 (MQ=255)
                                cGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACt                                                                 <  1:3891140/101‑1 (MQ=255)
                                cGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTCAACt                                                                 <  1:3890896/101‑1 (MQ=255)
                                 ggggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACtg                                                                <  1:3891246/101‑1 (MQ=255)
                                  gggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACActc                                                                           <  1:3890723/89‑1 (MQ=255)
                                  gggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACtgt                                                               <  2:3891123/101‑1 (MQ=255)
                                  gggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACtgt                                                               <  2:3891200/101‑1 (MQ=255)
                                   ggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACtgtg                                                              <  1:3891156/101‑1 (MQ=255)
                                    gAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGc                                                             <  2:3891231/101‑1 (MQ=255)
                                     aGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCg                                                            <  1:3891369/101‑1 (MQ=255)
                                     aGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCg                                                            >  1:3891541/1‑101 (MQ=255)
                                     aGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCg                                                            >  2:3891538/1‑101 (MQ=255)
                                     aGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCg                                                            <  1:3891220/101‑1 (MQ=255)
                                     aGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCg                                                            <  1:3891118/101‑1 (MQ=255)
                                       ttAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGAc                                                          <  2:3891216/101‑1 (MQ=255)
                                        tAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACa                                                         <  2:3891153/101‑1 (MQ=255)
                                        tAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACa                                                         >  2:3891543/1‑101 (MQ=255)
                                         aaaaaCAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACAt                                                        <  2:3891171/101‑1 (MQ=255)
                                         aaaaaCAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACAt                                                        >  2:3891546/1‑101 (MQ=255)
                                          aaaaCAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATg                                                       <  1:3891360/101‑1 (MQ=255)
                                           aaaCAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGt                                                      >  2:3891548/1‑101 (MQ=255)
                                            aacaacAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACAt                                                        <  1:3890916/98‑1 (MQ=255)
                                            aacaacAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGtt                                                     >  1:3891551/1‑101 (MQ=255)
                                             acaacaAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTc                                                    >  2:3891552/1‑101 (MQ=255)
                                              caacaaGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCg                                                   >  2:3891554/1‑101 (MQ=255)
                                                acaaGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGc                                                 <  1:3891027/101‑1 (MQ=255)
                                                  aaGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGt                                               <  2:3891022/101‑1 (MQ=255)
                                                   aGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGtt                                              <  1:3891340/101‑1 (MQ=255)
                                                    gggAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGttt                                             >  2:3891558/1‑101 (MQ=255)
                                                     ggAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTa                                            <  2:3891081/101‑1 (MQ=255)
                                                        gagGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGc                                         >  1:3891562/1‑101 (MQ=255)
                                                          ggTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATg                                          >  2:3891566/1‑98 (MQ=255)
                                                          ggTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGccc                                       >  2:3891564/1‑101 (MQ=255)
                                                             tattaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGa                                    <  1:3891128/101‑1 (MQ=255)
                                                             tattaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGa                                    <  1:3891393/101‑1 (MQ=255)
                                                              attaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                                   >  1:3891569/1‑101 (MQ=255)
                                                              attaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                                   <  1:3891190/101‑1 (MQ=255)
                                                               ttaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACtgtg                                                              >  2:3891575/1‑73 (MQ=255)
                                                               ttaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTAt                                           >  1:3891574/1‑92 (MQ=255)
                                                               ttaCTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATg                                  <  2:3891054/101‑1 (MQ=255)
                                                                  cTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCg                               >  2:3891579/1‑101 (MQ=255)
                                                                  cTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCg                               >  1:3891577/1‑101 (MQ=255)
                                                                     aGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATg                            >  1:3891583/1‑101 (MQ=255)
                                                                     aGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATg                            <  1:3891306/101‑1 (MQ=255)
                                                                     aGGGCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATg                            >  2:3891586/1‑101 (MQ=255)
                                                                      gggCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGc                           >  1:3891588/1‑101 (MQ=255)
                                                                      gggCGGCGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGc                           <  1:3891230/101‑1 (MQ=255)
                                                                       ggcggcGTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCg                          <  2:3891245/101‑1 (MQ=255)
                                                                         cggcgTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTg                        <  2:3891303/101‑1 (MQ=255)
                                                                         cggcgTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTg                        <  2:3891168/101‑1 (MQ=255)
                                                                           gcgTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGc                      <  2:3890990/101‑1 (MQ=255)
                                                                            cgTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCg                     >  1:3891592/1‑101 (MQ=255)
                                                                             gTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCgg                    >  2:3891596/1‑101 (MQ=255)
                                                                             gTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCgg                    >  1:3891595/1‑101 (MQ=255)
                                                                             gTTAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCgg                    <  1:3891453/101‑1 (MQ=255)
                                                                               tAGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTa                  >  2:3891598/1‑101 (MQ=255)
                                                                                aGCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTaa                 <  2:3891269/101‑1 (MQ=255)
                                                                                 gCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATTCCCTGATGCCGATGCGTGGCGGTaaa                <  2:3891092/101‑1 (MQ=255)
                                                                                 gCCATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTaaa                <  1:3891262/101‑1 (MQ=255)
                                                                                  ccATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTg                        >  2:3891603/1‑92 (MQ=255)
                                                                                  ccATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAAt               >  1:3891602/1‑101 (MQ=255)
                                                                                   cATTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATg              <  1:3891111/101‑1 (MQ=255)
                                                                                    aTTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATgg             <  1:3891201/101‑1 (MQ=255)
                                                                                    aTTGGCTGGAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGGAAATgg             >  2:3891605/1‑101 (MQ=255)
                                                                                            gAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATGGATTCGAGt     >  2:3891610/1‑101 (MQ=255)
                                                                                            gAAGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATGGATTCGAGt     <  1:3891318/101‑1 (MQ=255)
                                                                                             aaGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATGGATTCGAGtt    >  1:3891611/1‑101 (MQ=255)
                                                                                              aGGTGGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCAGATGCCGATGCGTGGCGGTAAATGGATTCGAGTTc   <  2:3891193/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATGGATTCGAGTTCA  >  NC_000913/2564605‑2564800

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑