Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,454,581 A→G intergenic (+3/+30) yjfF → / ← fbp galactofuranose ABC transporter putative membrane subunit YjtF/fructose‑1,6‑bisphosphatase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,454,5810AG98.5% 224.5 / ‑4.4 66intergenic (+3/+30)yjfF/fbpgalactofuranose ABC transporter putative membrane subunit YjtF/fructose‑1,6‑bisphosphatase 1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (37/28);  total (37/29)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCC  >  NC_000913/4454442‑4454716
                                                                                                                                           |                                                                                                                                       
cGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGGATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATaaaaa                                                                                                                                         >  2:1264107/1‑140 (MQ=255)
      gTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCctc                                                                                                                                   <  1:1264107/140‑1 (MQ=255)
       ttCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCctct                                                                                                                                  <  1:1097402/140‑1 (MQ=255)
       ttCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCctct                                                                                                                                  <  1:255913/140‑1 (MQ=255)
       ttCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCctct                                                                                                                                  <  2:251705/140‑1 (MQ=255)
       ttCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCctct                                                                                                                                  <  1:1376306/140‑1 (MQ=255)
               ggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTgg                                                                                                                          >  1:948853/1‑140 (MQ=255)
               ggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTgg                                                                                                                          >  1:641576/1‑140 (MQ=255)
               ggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTgg                                                                                                                          >  1:514901/1‑140 (MQ=255)
                   gAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGagag                                                                                                                      <  2:346209/140‑1 (MQ=255)
                    aaaaTCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAgg                                                                                                                     <  1:847524/140‑1 (MQ=255)
                                  tATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTAcg                                                                                                       <  2:1070430/140‑1 (MQ=255)
                                          tgttTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCggg                                                                                               >  2:570258/1‑140 (MQ=255)
                                            tttATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGaa                                                                                             >  2:196515/1‑140 (MQ=255)
                                                           ttACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCg                                                                              >  2:276929/1‑140 (MQ=255)
                                                           ttACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCg                                                                              >  2:183482/1‑140 (MQ=255)
                                                            tACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGt                                                                             >  2:51686/1‑140 (MQ=255)
                                                                gCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGa                                                                         <  1:51686/140‑1 (MQ=255)
                                                                 cGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGAc                                                                        >  2:512737/1‑140 (MQ=255)
                                                                 cGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGAc                                                                        >  2:126270/1‑140 (MQ=255)
                                                                 cGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGAc                                                                        >  2:51376/1‑140 (MQ=255)
                                                                   tGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACAt                                                                      >  1:1155929/1‑140 (MQ=255)
                                                                   tGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACAt                                                                      >  1:106486/1‑140 (MQ=255)
                                                                   tGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACAt                                                                      >  1:612303/1‑140 (MQ=255)
                                                                   tGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACAt                                                                      >  2:539732/1‑140 (MQ=255)
                                                                      aTTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCtt                                                                   >  1:533901/1‑140 (MQ=255)
                                                                      aTTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCtt                                                                   >  2:1063340/1‑140 (MQ=255)
                                                                         aaCGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCaa                                                                <  2:367573/140‑1 (MQ=255)
                                                                         aaCGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCaa                                                                <  2:533901/140‑1 (MQ=255)
                                                                                cTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATg                                                         >  2:920308/1‑140 (MQ=255)
                                                                                       agaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTg                                                  <  2:106486/140‑1 (MQ=255)
                                                                                        gaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTg                                                  >  1:1086008/1‑139 (MQ=255)
                                                                                        gaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTg                                                  <  2:1086008/139‑1 (MQ=255)
                                                                                        gaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGc                                                 >  2:442602/1‑140 (MQ=255)
                                                                                        gaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGc                                                 >  2:1154649/1‑140 (MQ=255)
                                                                                        gaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGc                                                 >  2:1214465/1‑140 (MQ=255)
                                                                                          aTCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCg                                               >  2:425495/1‑140 (MQ=255)
                                                                                          aTCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCg                                               >  1:1060194/1‑140 (MQ=255)
                                                                                            cGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGAc                                             <  1:183482/140‑1 (MQ=255)
                                                                                            cGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGAc                                             <  1:425495/140‑1 (MQ=255)
                                                                                               cAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACaaa                                          >  1:290850/1‑140 (MQ=255)
                                                                                                 gagTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGa                                        <  2:612303/140‑1 (MQ=255)
                                                                                                      cGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGaa                                   >  1:1357023/1‑140 (MQ=255)
                                                                                                      cGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGaa                                   <  1:570258/140‑1 (MQ=255)
                                                                                                         cAGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACgg                                <  2:290850/140‑1 (MQ=255)
                                                                                                          aGTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCC‑‑GTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgcgc                             <  2:659376/140‑1 (MQ=255)
                                                                                                           gTGACAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgcg                              <  1:176215/140‑1 (MQ=255)
                                                                                                              aCAAGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTg                           >  2:217648/1‑140 (MQ=255)
                                                                                                                aaGAGTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGt                         <  1:777315/140‑1 (MQ=255)
                                                                                                                 agagTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTg                        <  2:498725/140‑1 (MQ=255)
                                                                                                                   agTCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCa                      <  1:675030/140‑1 (MQ=255)
                                                                                                                     tCAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAgg                    >  1:191254/1‑140 (MQ=255)
                                                                                                                      cAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAggg                   >  1:1115627/1‑140 (MQ=255)
                                                                                                                      cAACATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAggg                   >  2:376775/1‑140 (MQ=255)
                                                                                                                       aaCATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGt                  >  2:177153/1‑140 (MQ=255)
                                                                                                                       aaCATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGt                  >  2:430952/1‑140 (MQ=255)
                                                                                                                         cATTGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGttt                <  2:641576/140‑1 (MQ=255)
                                                                                                                            tGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCg             >  1:1040726/1‑140 (MQ=255)
                                                                                                                            tGCGCAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCg             <  1:920308/140‑1 (MQ=255)
                                                                                                                             gcgcAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCa                      <  1:671795/130‑1 (MQ=255)
                                                                                                                             gcgcAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCa                      >  2:671795/1‑130 (MQ=255)
                                                                                                                             gcgcAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCgg            >  1:150035/1‑140 (MQ=255)
                                                                                                                              cgcAGCAATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCggg           <  1:177153/140‑1 (MQ=255)
                                                                                                                                cagcaATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat         <  2:1357023/140‑1 (MQ=255)
                                                                                                                                   caATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgatgat      >  1:202993/1‑140 (MQ=255)
                                                                                                                                   caATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgatgat      <  1:1063340/140‑1 (MQ=255)
                                                                                                                                   caATAAAAGCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgatgat      >  1:354472/1‑140 (MQ=255)
                                                                                                                                       aaaagCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATcc  <  2:1060194/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
CGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCC  >  NC_000913/4454442‑4454716

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: