Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,561,503 A→G intergenic (+366/‑94) b2447 → / → b2448 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,5030AG96.8% 570.6 / 4.6 189intergenic (+366/‑94)b2447/b2448hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (118/64);  total (122/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 196 total aligned reads displayed.

CACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATT  >  NC_000913/2561413‑2561591
                                                                                          |                                                                                        
caccaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCTCTTTTTCCCTTCAGGCATTTTGGTTTAATTCAATCAACGAAAGGCGGt                                                                                        <  2:166285/93‑1 (MQ=255)
caccaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGt                                                                                        >  2:1640490/1‑93 (MQ=255)
 accaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGtt                                                                                       <  2:5425/93‑1 (MQ=255)
 accaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGtt                                                                                       >  2:3813703/1‑93 (MQ=255)
  ccaAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTc                                                                                      >  2:12116/1‑93 (MQ=255)
   caAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTcc                                                                                     <  1:2438002/93‑1 (MQ=255)
       gCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaac                                                                                 >  1:2043284/1‑93 (MQ=255)
         caacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaa                                                                               >  1:1135827/1‑93 (MQ=255)
         caacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaa                                                                               >  1:374698/1‑93 (MQ=255)
          aacaaGATGGGAGTTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaac                                                                              <  1:4714264/93‑1 (MQ=255)
          aacaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaac                                                                              <  1:733494/93‑1 (MQ=255)
           acaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCATTCAACGAAAGGCGGCTCCcaacaaca                                                                             >  1:3964467/1‑93 (MQ=255)
           acaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaaca                                                                             <  2:2283153/93‑1 (MQ=255)
           acaaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCcaacaaca                                                                             >  2:254682/1‑93 (MQ=255)
             aaGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACaaa                                                                           <  2:976790/93‑1 (MQ=255)
             aaGATGGGAGGCGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACaaa                                                                           >  2:733674/1‑93 (MQ=255)
              agagGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACaaaa                                                                          <  1:2479152/93‑1 (MQ=255)
              aGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACaaaa                                                                          <  1:3037539/93‑1 (MQ=255)
                 tGGGAGGTGAATTTTTTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCGCAACAACAAAATaa                                                                       <  2:4184238/93‑1 (MQ=255)
                 tGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCGACGAAAGGCGGTTCCCAACAACAAAATaa                                                                       <  1:893479/93‑1 (MQ=255)
                 tGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaa                                                                       >  1:3109572/1‑93 (MQ=255)
                 tGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaa                                                                       >  2:4524132/1‑93 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                      >  1:2804372/1‑93 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                      <  2:1043676/93‑1 (MQ=255)
                  gggAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATaaa                                                                      >  1:4421797/1‑93 (MQ=255)
                     aGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGAt                                                                   <  2:4746378/93‑1 (MQ=255)
                      ggTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATg                                                                  >  1:3303049/1‑93 (MQ=255)
                           aTTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAgg                                                             >  2:3914231/1‑93 (MQ=255)
                              ttGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaa                                                          >  2:672361/1‑93 (MQ=255)
                               tgtgATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTaaa                                                         <  1:2944966/93‑1 (MQ=255)
                                  gATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAAt                                                      >  2:2899945/1‑93 (MQ=255)
                                     gCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGaa                                                   >  2:4423638/1‑93 (MQ=255)
                                        gATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCCACAACAAACTAAAGATGGGAGGTAAACATGAAACg                                                >  2:2635272/1‑93 (MQ=255)
                                         aTCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACgg                                               <  2:2450679/93‑1 (MQ=255)
                                          tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGt                                              >  1:2347863/1‑93 (MQ=255)
                                          tCGCTTTTTCCCTTCAGGCATTTTTGTTTAATGCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGt                                              <  1:525055/93‑1 (MQ=255)
                                            gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTgg                                            <  1:4190259/93‑1 (MQ=255)
                                            gCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTgg                                            >  1:1021052/1‑93 (MQ=255)
                                             cTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTggg                                           >  1:3895131/1‑93 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGa                                          <  1:1862718/93‑1 (MQ=255)
                                              tttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTGGCAACAACAAAATAAAGATGGGAGGGAAAAATGAAACGGTGGGa                                          >  1:2971893/1‑93 (MQ=255)
                                               ttttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGag                                         >  2:1191729/1‑92 (MQ=255)
                                                tttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAg                                        <  2:1709055/93‑1 (MQ=255)
                                                tttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAg                                        >  2:2718/1‑93 (MQ=255)
                                                 ttCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGt                                       <  2:761213/93‑1 (MQ=255)
                                                  tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTg                                      >  2:648135/1‑93 (MQ=255)
                                                  tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTg                                      <  2:1498259/93‑1 (MQ=255)
                                                  tCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTg                                      >  2:1041984/1‑93 (MQ=255)
                                                   cccTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTgg                                     >  1:1565131/1‑93 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGTAGTGGCt                                   >  2:4340224/1‑93 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCt                                   >  2:4669747/1‑93 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCt                                   >  2:4726755/1‑93 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCt                                   >  2:3428619/1‑93 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCt                                   <  2:3772922/93‑1 (MQ=255)
                                                     cTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCt                                   >  2:862247/1‑93 (MQ=255)
                                                      ttCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCtt                                  >  1:3579960/1‑93 (MQ=255)
                                                        cAGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTa                                <  2:182909/93‑1 (MQ=255)
                                                         aGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTAc                               <  2:3623249/93‑1 (MQ=255)
                                                         aGGCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTAc                               >  1:2402343/1‑93 (MQ=255)
                                                          ggCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACt                              >  2:3280770/1‑93 (MQ=255)
                                                           gCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTg                             <  2:2605517/93‑1 (MQ=255)
                                                           gCATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTg                             >  1:3394562/1‑93 (MQ=255)
                                                            cATTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTgg                            >  2:520502/1‑93 (MQ=255)
                                                             aTTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAGATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                           <  2:2990223/93‑1 (MQ=255)
                                                             aTTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                           >  1:811365/1‑93 (MQ=255)
                                                             aTTTTTGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGt                           >  1:290746/1‑93 (MQ=255)
                                                              tttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                          >  2:1318682/1‑93 (MQ=255)
                                                              tttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTg                          >  2:920379/1‑93 (MQ=255)
                                                              tttttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGGAAAAATGAAACGGGGGGAAGTGGCTTTACTGGTg                          >  1:3131427/1‑93 (MQ=255)
                                                                 ttGTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATg                       <  2:2662854/93‑1 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATgga                     >  1:2174951/1‑92 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                     >  1:3800630/1‑93 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                     >  1:4158892/1‑93 (MQ=255)
                                                                   gTTTAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGt                     >  2:3283140/1‑93 (MQ=255)
                                                                    tttAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTa                    >  2:1559246/1‑93 (MQ=255)
                                                                     ttAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                   <  2:811366/93‑1 (MQ=255)
                                                                     ttAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTaa                   <  1:2697030/93‑1 (MQ=255)
                                                                      tAATTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAAt                  >  2:3217457/1‑93 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                 <  1:1279686/93‑1 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                 >  1:1163963/1‑93 (MQ=255)
                                                                       aaTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATg                 <  2:1804888/93‑1 (MQ=255)
                                                                        aTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                <  2:3468893/93‑1 (MQ=255)
                                                                        aTTCAATCAACGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGt                >  2:1440574/1‑93 (MQ=255)
                                                                        aTTCAATCAACGAAAGGCGCTTCCAAAAAACAAAATAAAGATAGGAGGTAAAAGTGAAACGGTGGGAAGTGGCTTTACTGGTGATAGTAATGt                <  2:3986352/93‑1 (MQ=255)
                                                                         atcaatcaaCGAAAGGCCGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGGGGGAAGTGGCTTTACTGGTGATGGTAATGtt               <  1:572606/92‑1 (MQ=255)
                                                                          tcaatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg              >  1:1064282/1‑93 (MQ=255)
                                                                          tcaatcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTg              >  2:2042427/1‑93 (MQ=255)
                                                                             atcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt           <  2:588096/93‑1 (MQ=255)
                                                                             atcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt           >  2:2159056/1‑93 (MQ=255)
                                                                             atcaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGtt           >  1:4356866/1‑93 (MQ=255)
                                                                              acaaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGttt          <  1:65670/92‑1 (MQ=255)
                                                                               caaCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTg         >  1:1519501/1‑93 (MQ=255)
                                                                                 aCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt       >  2:1956167/1‑93 (MQ=255)
                                                                                 aCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt       >  2:1080314/1‑93 (MQ=255)
                                                                                 aCGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCt       >  2:4050875/1‑93 (MQ=255)
                                                                                  cGAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTg      >  1:1808610/1‑93 (MQ=255)
                                                                                   gAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt     <  2:407018/93‑1 (MQ=255)
                                                                                   gAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGt     >  2:3720014/1‑93 (MQ=255)
                                                                                   gAAAGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGCAATGTTGGTTtgttgt     >  1:3373360/1‑93 (MQ=255)
                                                                                    aaaGGCGGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTa    >  1:1924033/1‑93 (MQ=255)
                                                                                          |                                                                                        
CACCAAAGCCAACAAGATGGGAGGTGAATTTTGTGATGCAGATCGCTTTTTCCCTTCAGGCATTTTTGTTTAATTCAATCAACGAAAGGCAGTTCCCAACAACAAAATAAAGATGGGAGGTAAAAATGAAACGGTGGGAAGTGGCTTTACTGGTGATGGTAATGTTGGTTTGCTGTATT  >  NC_000913/2561413‑2561591

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