Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,070,707 C→T 63.4% G180S (GGT→AGT)  fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,070,7070CT63.4% 134.1 / 148.8 142G180S (GGT→AGT) fbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base C (21/31);  new base T (40/50);  total (61/81)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 146 total aligned reads displayed.

CCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCG  >  NC_000913/3070562‑3070852
                                                                                                                                                 |                                                                                                                                                 
ccGTGTACGTTACCGAAGGACGTTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGt                                                                                                                                                <  3:820464/149‑1 (MQ=255)
ccGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGt                                                                                                                                                <  1:1109574/149‑1 (MQ=255)
        gTTACCGAAGGGCGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGATGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCa                                                                                                                                        >  3:976116/1‑149 (MQ=255)
         ttACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAg                                                                                                                                       >  4:1387794/1‑149 (MQ=255)
         ttACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAg                                                                                                                                       >  3:1324532/1‑149 (MQ=255)
                aGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGAtt                                                                                                                                <  3:1374348/149‑1 (MQ=255)
                aGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGAtt                                                                                                                                <  1:622460/149‑1 (MQ=255)
                     gCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCag                                                                                                                           >  1:513793/1‑149 (MQ=255)
                     gCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCag                                                                                                                           >  4:228795/1‑149 (MQ=255)
                      cTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCaga                                                                                                                          <  2:1247349/149‑1 (MQ=255)
                      cTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCaga                                                                                                                          <  4:1200264/149‑1 (MQ=255)
                      cTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCaga                                                                                                                          <  1:1059872/149‑1 (MQ=255)
                          gATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCa                                                                                                                      >  2:167568/1‑149 (MQ=255)
                            tGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATg                                                                                                                    <  2:28714/149‑1 (MQ=255)
                               tGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCg                                                                                                                 >  2:376583/1‑149 (MQ=255)
                                    cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAtttt                                                                                                            <  2:1013774/149‑1 (MQ=255)
                                    cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAtttt                                                                                                            <  4:392954/149‑1 (MQ=255)
                                    cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAtttt                                                                                                            <  3:566750/149‑1 (MQ=255)
                                    cgcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAtttt                                                                                                            <  4:1254488/149‑1 (MQ=255)
                                     gcgGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTg                                                                                                           <  1:942321/149‑1 (MQ=255)
                                       gggCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCa                                                                                                                            >  3:973468/1‑130 (MQ=255)
                                       gggCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCa                                                                                                                            <  4:973468/130‑1 (MQ=255)
                                         gCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACa                                                                                                       >  1:1469965/1‑149 (MQ=255)
                                          cTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACAt                                                                                                      <  1:1115457/149‑1 (MQ=255)
                                           tGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATg                                                                                                     >  4:1107179/1‑149 (MQ=255)
                                            gATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATgc                                                                                                    <  1:167568/149‑1 (MQ=255)
                                            gATTTTGCTCAGTTCGGTGAATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATgc                                                                                                    <  4:1010974/149‑1 (MQ=255)
                                             aTTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATgcg                                                                                                   >  4:1134221/1‑149 (MQ=255)
                                              ttttGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATgcgc                                                                                                  >  1:967218/1‑149 (MQ=255)
                                              ttttGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATgcgc                                                                                                  <  1:1348102/149‑1 (MQ=255)
                                              ttttGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATgcgc                                                                                                  >  4:640681/1‑149 (MQ=255)
                                                   cTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGCCATGCCGATTTTgg                                                                                                          >  3:1360698/1‑136 (MQ=255)
                                                   cTCAGTTCGGTGTATGCGTAATCAACGTCTGCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAg                                                                                             <  4:137289/149‑1 (MQ=255)
                                                    tCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAgg                                                                                            >  2:1080879/1‑149 (MQ=255)
                                                      aGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTcc                                                                                                                             <  2:1437437/114‑1 (MQ=255)
                                                      aGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTcc                                                                                                                             >  1:1437437/1‑114 (MQ=255)
                                                      aGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTa                                                                                          <  1:546999/149‑1 (MQ=255)
                                                        ttCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTAtt                                                                                        <  1:651660/149‑1 (MQ=255)
                                                        ttCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCGCCAGGTAtt                                                                                        >  2:1013221/1‑149 (MQ=255)
                                                         tCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTAttt                                                                                       <  1:1494973/149‑1 (MQ=255)
                                                                 tGCGTAATCAACGTGTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGa                                                                               >  1:1134807/1‑149 (MQ=255)
                                                                 tGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGa                                                                               >  3:828714/1‑149 (MQ=255)
                                                                 tGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGa                                                                               <  2:719540/149‑1 (MQ=255)
                                                                       aTCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGAttt                                                                                                             >  4:98729/1‑113 (MQ=255)
                                                                       aTCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGa                                                                         >  1:1081517/1‑149 (MQ=255)
                                                                        tCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGAt                                                                        <  3:1089500/149‑1 (MQ=255)
                                                                           aCGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGtt                                                                     >  3:908962/1‑149 (MQ=255)
                                                                           aCGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGtt                                                                     >  3:1425358/1‑149 (MQ=255)
                                                                             gTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTtct                                                                   >  4:1002230/1‑149 (MQ=255)
                                                                             gTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCa                                                                                              <  3:478592/122‑1 (MQ=255)
                                                                             gTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCa                                                                                              >  4:478592/1‑122 (MQ=255)
                                                                                ttCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCtt                                                                >  4:762837/1‑149 (MQ=255)
                                                                                ttCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCtt                                                                <  1:94771/149‑1 (MQ=255)
                                                                                ttCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCtt                                                                >  3:1217113/1‑149 (MQ=255)
                                                                                ttCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCtt                                                                >  1:371862/1‑149 (MQ=255)
                                                                                ttCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCtt                                                                >  1:777935/1‑149 (MQ=255)
                                                                                ttCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCtt                                                                <  2:538433/149‑1 (MQ=255)
                                                                                 tCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTg                                                               <  3:228795/149‑1 (MQ=255)
                                                                                  ccGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGc                                                              <  1:263666/149‑1 (MQ=255)
                                                                                   cGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCa                                                             <  3:1421811/149‑1 (MQ=255)
                                                                                   cGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCa                                                             <  2:1134807/149‑1 (MQ=255)
                                                                                   cGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCa                                                             <  3:287043/149‑1 (MQ=255)
                                                                                    ggCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTcc                                                                                                                             >  4:1033104/1‑84 (MQ=255)
                                                                                    ggCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTcc                                                                                                                             <  3:1033104/84‑1 (MQ=255)
                                                                                     gCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCaga                                                           <  2:371862/149‑1 (MQ=255)
                                                                                      cTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCagag                                                          >  3:86703/1‑149 (MQ=255)
                                                                                       tGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAgaga                                                         >  3:330210/1‑149 (MQ=255)
                                                                                            gtACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGAttctt                                                    <  3:870826/149‑1 (MQ=255)
                                                                                            gtACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGAttctt                                                    <  3:164590/149‑1 (MQ=255)
                                                                                            gtACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGAttctt                                                    <  1:1299440/149‑1 (MQ=255)
                                                                                              aCAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGACTCTTCa                                                  <  1:1464967/149‑1 (MQ=255)
                                                                                                aGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGa                                                <  2:400658/149‑1 (MQ=255)
                                                                                                   gCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAg                                             >  1:222317/1‑149 (MQ=255)
                                                                                                     aGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGt                                           <  3:640681/149‑1 (MQ=255)
                                                                                                      gAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTc                                          <  2:431187/149‑1 (MQ=255)
                                                                                                         gCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTTGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGAt                                       >  3:1143384/1‑149 (MQ=255)
                                                                                                         gCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGAt                                       <  3:170135/149‑1 (MQ=255)
                                                                                                          cGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATc                                      >  2:1173535/1‑149 (MQ=255)
                                                                                                           gTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCa                                     <  3:1134221/149‑1 (MQ=255)
                                                                                                            tCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCAt                                    <  3:1320664/149‑1 (MQ=255)
                                                                                                                tgtgGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATCTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTga                                <  3:887929/149‑1 (MQ=255)
                                                                                                                  tgGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTgaga                              <  1:1246598/149‑1 (MQ=255)
                                                                                                                  tgGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTga                                <  2:75854/147‑1 (MQ=255)
                                                                                                                       gttgtCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTgagaagaga                         <  4:188539/149‑1 (MQ=255)
                                                                                                                        ttgtCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAATCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTgagaagagaa                        >  2:1294133/1‑149 (MQ=255)
                                                                                                                          gtCCACGCCGTCTTCTTCACCACTGGTGAAACCCAGTTCGATTTCCAGAGTCATGCCGATGTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCCCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACa                      <  4:167523/149‑1 (MQ=255)
                                                                                                                          gtCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACa                      <  3:95657/149‑1 (MQ=255)
                                                                                                                           tCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTgagaagag                          >  1:57352/1‑144 (MQ=255)
                                                                                                                           tCCACGCCGTCTTCTTCACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTgagaagag                          <  2:57352/144‑1 (MQ=255)
                                                                                                                                   gTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTAc             <  1:775042/149‑1 (MQ=255)
                                                                                                                                      ttcttcACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTtc                                                                    <  4:840132/91‑1 (MQ=255)
                                                                                                                                      ttcttcACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCgg          >  3:457187/1‑149 (MQ=255)
                                                                                                                                      ttcttcACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCgg          >  2:900059/1‑149 (MQ=255)
                                                                                                                                      ttcttcACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCgg          >  1:1510476/1‑149 (MQ=255)
                                                                                                                                      ttcttcACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCgg          <  4:828714/149‑1 (MQ=255)
                                                                                                                                        cttcACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTa        <  4:86703/149‑1 (MQ=255)
                                                                                                                                         ttcACCACTGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAg       >  2:1223059/1‑149 (MQ=255)
                                                                                                                                            accaccGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTg    <  4:1324532/149‑1 (MQ=255)
                                                                                                                                             ccaccGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGc   >  2:880641/1‑149 (MQ=255)
                                                                                                                                              cactGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCg  >  4:628537/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCGTGTACGTTACCGAAGGACGCTGCGATGGTGAAACGCGGGCTGATTTTGCTCAGTTCGGTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCG  >  NC_000913/3070562‑3070852

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