Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,237,247 +C intergenic (+60/‑337) ygjT → / → ygjU putative transport protein/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,237,2451.C99.7% 931.1 / ‑4.2 299intergenic (+58/‑339)ygjT/ygjUputative transport protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (157/141);  total (157/142)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 346 total aligned reads displayed.

AATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGG  >  NC_000913/3237155‑3237332
                                                                                            |                                                                                       
 aTTATCGGCATGATAAGCAGCGTGGTGGATAATTTGTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                                       <  2:5861947/93‑1 (MQ=255)
 aTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                                       <  2:3809685/93‑1 (MQ=255)
 aTTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcc                                                                                       <  1:4834073/93‑1 (MQ=255)
  ttATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccc                                                                                      <  2:3679705/93‑2 (MQ=255)
   tATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccct                                                                                     <  2:169206/93‑3 (MQ=255)
    aTCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctccctt                                                                                    >  1:7531376/1‑90 (MQ=255)
     tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcccttc                                                                                   >  1:1882748/1‑89 (MQ=255)
     tCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTctccctctcccttc                                                                                   <  2:6792034/93‑5 (MQ=255)
       ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                  >  2:5637997/1‑93 (MQ=255)
       ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                  <  2:3550646/93‑1 (MQ=255)
       ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                  <  2:2366199/93‑1 (MQ=255)
       ggCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                  >  1:7502375/1‑93 (MQ=255)
       ggCATGATAAGCAGCGGGGTAGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCa                                                                                  <  2:3159441/93‑1 (MQ=255)
        gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                 <  2:7929238/93‑1 (MQ=255)
        gCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAg                                                                                 >  1:4252135/1‑93 (MQ=255)
         cATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAgg                                                                                <  2:2826623/93‑1 (MQ=255)
            gATAATCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                             <  1:3101511/93‑1 (MQ=255)
            gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                             >  1:3721925/1‑93 (MQ=255)
            gATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgag                                                                             >  2:4648278/1‑93 (MQ=255)
              tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaggg                                                                           >  2:7302178/1‑91 (MQ=255)
              tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaggg                                                                           <  2:4970974/93‑3 (MQ=255)
              tAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGgaggg                                                                           <  1:1079218/93‑3 (MQ=255)
               aaGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAgggg                                                                          <  1:7059131/93‑1 (MQ=255)
                 gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGt                                                                        >  2:471431/1‑93 (MQ=255)
                 gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGt                                                                        <  1:6016930/93‑1 (MQ=255)
                 gcagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGt                                                                        <  1:5328114/93‑1 (MQ=255)
                  cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTc                                                                       <  2:5284295/93‑1 (MQ=255)
                  cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTc                                                                       <  2:7596122/93‑1 (MQ=255)
                  cagcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTc                                                                       <  1:6339181/93‑1 (MQ=255)
                   agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCg                                                                      <  2:1871693/93‑1 (MQ=255)
                   agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGGCg                                                                      >  2:3148419/1‑93 (MQ=255)
                     cGGGGTGGATAATTTTTAATCTGACTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCggg                                                                    <  2:4325242/93‑1 (MQ=255)
                      ggggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCgggg                                                                   >  2:7795071/1‑93 (MQ=255)
                      ggggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCgggg                                                                   <  2:3871391/93‑1 (MQ=255)
                       gggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGt                                                                  <  2:6441481/93‑1 (MQ=255)
                       gagtgtATAATTGTTCATCTGCCTAAGCCGTGTACTCTGTCATTAACATGAGCTCCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGt                                                                  <  1:3833858/91‑1 (MQ=255)
                        ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTg                                                                 >  2:200856/1‑93 (MQ=255)
                        ggTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTg                                                                 >  1:4051178/1‑93 (MQ=255)
                         gTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGa                                                                >  1:5094985/1‑93 (MQ=255)
                          tGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAg                                                               <  1:7791367/93‑1 (MQ=255)
                           tgATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTTATTAACATGAGCGCCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAgg                                                              <  1:6731205/92‑1 (MQ=255)
                            gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAggg                                                             >  1:2614317/1‑93 (MQ=255)
                            gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAggg                                                             >  2:1535016/1‑93 (MQ=255)
                            gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAggg                                                             >  2:7171678/1‑93 (MQ=255)
                            gATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAggg                                                             >  2:6191550/1‑93 (MQ=255)
                               aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTaa                                                          >  1:158725/1‑93 (MQ=38)
                                aTTTTTAATCTGCGTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAt                                                         <  1:4640998/93‑1 (MQ=38)
                                aTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAt                                                         <  1:3928297/93‑1 (MQ=255)
                                 tttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAtt                                                        >  1:5883218/1‑93 (MQ=38)
                                  ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAttt                                                       <  1:6250048/93‑1 (MQ=38)
                                   tttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAtttt                                                      >  1:7363947/1‑93 (MQ=38)
                                    ttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAAttttt                                                     >  1:971613/1‑93 (MQ=38)
                                     tAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTc                                                    <  2:6472746/93‑1 (MQ=38)
                                       aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGc                                                  >  1:3548255/1‑93 (MQ=38)
                                       aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGc                                                  >  1:5455139/1‑93 (MQ=38)
                                        tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCa                                                 >  2:4880606/1‑93 (MQ=38)
                                        tCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCa                                                 >  1:1992333/1‑93 (MQ=38)
                                          tGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCAcc                                               >  1:5596001/1‑93 (MQ=38)
                                           gCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCg                                              <  2:929933/93‑1 (MQ=38)
                                            ccTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGa                                             >  1:4171754/1‑93 (MQ=38)
                                             cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                            <  2:1789961/93‑1 (MQ=38)
                                             cTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGAt                                            >  1:6653596/1‑93 (MQ=38)
                                              tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                                           >  2:7384960/1‑93 (MQ=38)
                                              tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                                           >  1:4563663/1‑93 (MQ=38)
                                              tAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATg                                           >  2:5156835/1‑93 (MQ=38)
                                                aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCt                                         <  1:7179488/93‑1 (MQ=38)
                                                aGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCt                                         <  1:7113252/93‑1 (MQ=38)
                                                   cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc                                      <  2:6754552/93‑1 (MQ=38)
                                                   cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGc                                      >  1:4360793/1‑93 (MQ=37)
                                                   cGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGCGAGGGGAATTTTTCGCACCGATGCTGGc                                      >  2:2682821/1‑93 (MQ=37)
                                                    gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                                     >  1:1027574/1‑93 (MQ=37)
                                                    gtgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGcc                                     <  1:3134098/93‑1 (MQ=38)
                                                     tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCt                                    >  2:4024046/1‑93 (MQ=37)
                                                     tgtACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCt                                    >  2:3471239/1‑93 (MQ=37)
                                                       tACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                                  >  1:804229/1‑93 (MQ=37)
                                                       tACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGt                                  <  1:7464845/93‑1 (MQ=37)
                                                        aCCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                                 >  1:7709780/1‑93 (MQ=37)
                                                         cccTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTc                                >  1:4450884/1‑93 (MQ=37)
                                                            tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTtccc                             >  1:2948554/1‑93 (MQ=37)
                                                            tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                             >  1:883915/1‑93 (MQ=37)
                                                            tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                             <  2:5384047/93‑1 (MQ=37)
                                                            tGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                             <  1:3341867/93‑1 (MQ=38)
                                                               cATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCa                          >  2:4137566/1‑93 (MQ=25)
                                                                aTTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAc                         >  1:5174186/1‑93 (MQ=25)
                                                                    aCATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTa                     <  1:1373388/93‑1 (MQ=21)
                                                                     cATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCCACCTaa                    >  1:2443152/1‑93 (MQ=21)
                                                                         aGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCct                >  2:1179429/1‑93 (MQ=21)
                                                                            aCCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc             <  1:2235652/93‑1 (MQ=18)
                                                                            aCCGTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCtctc             <  2:2936656/93‑1 (MQ=18)
                                                                               gTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTtccc          >  2:4041590/1‑93 (MQ=18)
                                                                               gTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTcccc          <  1:839167/93‑1 (MQ=18)
                                                                               gTTTTCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCCAACCCTCTcccc          >  2:6296694/1‑93 (MQ=18)
                                                                                ttttCTCCCTCTCCCTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTTTTCTCACCCCCACCCTccccccc         >  2:1428825/1‑92 (MQ=12)
                                                                                   tctccctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc      >  2:2563493/5‑93 (MQ=18)
                                                                                   tctccctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAc      <  1:5818288/89‑1 (MQ=18)
                                                                                     tccctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg    <  2:6626554/91‑1 (MQ=18)
                                                                                     tccctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg    >  2:5063398/3‑93 (MQ=18)
                                                                                     tccctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg    >  1:1217994/3‑93 (MQ=18)
                                                                                     cccctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgg    <  2:6965174/91‑1 (MQ=18)
                                                                                       cctctcccTTCCAGGGAGGGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACgggg  >  1:2744138/1‑93 (MQ=18)
                                                                                            |                                                                                       
AATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCTCTCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGG  >  NC_000913/3237155‑3237332

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