Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,337,755 T→C 5.8% L642S (TTA→TCA)  acnA → aconitate hydratase 1; aconitase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,337,7550TC 5.8% 204.0 / 3.4 69L642S (TTA→TCA) acnAaconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base C (2/2);  total (37/32)
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.14e-01

ATTGCCCGTGCGGTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCAGCGGGTCGAT  >  NC_000913/1337610‑1337904
                                                                                                                                                 |                                                                                                                                                     
aTTGCCCGTGCGGTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCg                                                                                                                                                   <  1:1156148/150‑1 (MQ=255)
aTTGCCCGTGCGGTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCAAACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCg                                                                                                                                                   >  1:780486/1‑150 (MQ=255)
  tGCCCGTGCGGTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGcc                                                                                                                                                 >  1:686580/1‑150 (MQ=255)
          cGGTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTtg                                                                                                                                         >  2:621446/1‑150 (MQ=255)
          cGGTAGAACAAGTCTCCACAGAAATGTACCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACATACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATAGCctttctttt                                                                                                                                         >  1:861800/1‑149 (MQ=255)
            ggaGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTtgat                                                                                                                                       <  1:574765/148‑1 (MQ=255)
            gTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTtgat                                                                                                                                       >  2:570996/1‑150 (MQ=255)
            gTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTtgat                                                                                                                                       <  1:381557/150‑1 (MQ=255)
            gTAGAACAAGTCTCCACAGAAATGTTCCCCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTtgat                                                                                                                                       <  1:850959/150‑1 (MQ=255)
              aGAACAAGTCTCAACAGAAATGTTCCGCAAAGAGTACGCAGCAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGAACCGATACCTACGGTTTGCAGGAGGACTCAACATATATTCGCTTATCGCCTTTCTTtgatga                                                                                                                                     >  1:540635/1‑150 (MQ=255)
               gAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGATGGGAATTAACGTCACACGATCCGATACCTACGGTTTGCAGGAGGACTCAACCTATATTCGCTTATCGCCGTTCTTTGATGaa                                                                                                                                    >  1:421024/1‑150 (MQ=255)
                          cacaGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAac                                                                                                                         <  1:134947/150‑1 (MQ=255)
                            caGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAacac                                                                                                                       <  2:220286/150‑1 (MQ=255)
                              gAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAAcacca                                                                                                                     <  2:742417/150‑1 (MQ=255)
                                   gTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACATGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAAcaccagcacc                                                                                                                >  1:385800/1‑150 (MQ=255)
                                       cGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTg                                                                                                            <  1:669996/150‑1 (MQ=255)
                                       cGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACAAGTg                                                                                                            >  2:666753/1‑150 (MQ=255)
                                        gCAAAGAGTACGCAGAAGTTTTTGAATGCACAGCAGACGGGAAGGGAATTAACGGCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTgg                                                                                                           <  2:874547/150‑1 (MQ=255)
                                        gCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTgg                                                                                                           >  2:971296/1‑150 (MQ=255)
                                         cAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGa                                                                                                          <  1:1199947/150‑1 (MQ=255)
                                                  cGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCa                                                                                                 >  2:1141843/1‑150 (MQ=255)
                                                    cAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACTGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCaag                                                                                               >  1:199826/1‑148 (MQ=255)
                                                    cAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACg                                                                                               >  2:1015416/1‑150 (MQ=255)
                                                    cAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTCATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACg                                                                                               >  2:879813/1‑150 (MQ=255)
                                                    cAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTCATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACg                                                                                               >  2:875331/1‑150 (MQ=255)
                                                          tttttGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTgcgc                                                                                         >  1:589111/1‑150 (MQ=255)
                                                          tttttGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTgcgc                                                                                         >  1:776593/1‑150 (MQ=255)
                                                          tttttGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTgcgc                                                                                         <  2:712683/150‑1 (MQ=255)
                                                               gAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGGTGGCAGGAGGACTCACCCTATATTCGCAGATCGCCGTTCTGTGATGAAATGCAGCCAACACCAGCACCAGTGCAAGATCTTCACGGTGCGCGGCTc                                                                                    >  1:226137/1‑150 (MQ=255)
                                                               gAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGAGCCGATCCCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAATTGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATc                                                                                    <  1:895651/150‑1 (MQ=255)
                                                                aaGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATcc                                                                                   <  1:141375/150‑1 (MQ=255)
                                                                  ggCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTc                                                                                 >  1:1102381/1‑150 (MQ=255)
                                                                  ggCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTc                                                                                 >  2:76449/1‑150 (MQ=255)
                                                                  ggCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTc                                                                                 >  2:1179736/1‑150 (MQ=255)
                                                                            agTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTgg                                                                       >  2:932712/1‑150 (MQ=255)
                                                                             gTGGAAGGGAATTAACGTCACACGATCCGATAACTACGGTTGTCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTggg                                                                      <  2:328176/150‑1 (MQ=255)
                                                                               ggAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTggggg                                                                    <  2:623436/150‑1 (MQ=255)
                                                                                 aaGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGAt                                                                  >  2:1047723/1‑150 (MQ=255)
                                                                                 aaGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGAt                                                                  >  2:512860/1‑150 (MQ=255)
                                                                                    ggAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATAATCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCa                                                               >  1:1017584/1‑150 (MQ=255)
                                                                                      aaTTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTCATCGCCTTTCTTTGATGAAATGGAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGt                                                             <  1:879813/150‑1 (MQ=255)
                                                                                      aaTTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTCATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGt                                                             <  1:875331/150‑1 (MQ=255)
                                                                                        ttAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTca                                                           >  1:434662/1‑150 (MQ=255)
                                                                                               cacaCGATCCGATACCTACGGTTGGCAGGATGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGATATGCAGGCAACAGCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGa                                                    <  1:652839/150‑1 (MQ=255)
                                                                                               cacaCGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGa                                                    >  2:265570/1‑150 (MQ=255)
                                                                                                 cacGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGAAATGCTGGGGGATTCAGTCACCACTGAcc                                                  <  2:1175960/150‑1 (MQ=255)
                                                                                                  acGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCa                                                 >  2:245848/1‑150 (MQ=255)
                                                                                                          aTACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAAAGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATAtctctc                                         <  2:11921/150‑1 (MQ=255)
                                                                                                           tACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGAGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTcc                                        <  2:194531/150‑1 (MQ=255)
                                                                                                             ccTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCCCGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCCCCACTGACCATATCTCTCcgg                                      >  1:1215708/1‑150 (MQ=255)
                                                                                                             ccTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGCCCATATCTCTCcgg                                      >  1:321709/1‑150 (MQ=255)
                                                                                                             ccTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCcgg                                      >  1:788077/1‑150 (MQ=255)
                                                                                                             ccTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCcgg                                      <  2:686580/150‑1 (MQ=255)
                                                                                                                  ggTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTGCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGc                                 <  1:669198/150‑1 (MQ=255)
                                                                                                                       gCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTAt                            <  2:1225822/150‑1 (MQ=255)
                                                                                                                         aggaggACTCAACCTATATTCGCTTATCGCCTTTCTTAGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTa                          <  2:199826/150‑1 (MQ=255)
                                                                                                                          ggaggaCTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTaa                         <  2:540635/150‑1 (MQ=255)
                                                                                                                            aggaCTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCTGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGc                       <  2:489116/150‑1 (MQ=255)
                                                                                                                               aCTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCg                    >  2:743087/1‑150 (MQ=255)
                                                                                                                                 tCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGAc                  <  2:538495/150‑1 (MQ=255)
                                                                                                                                  cAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCGCCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACa                 >  2:252940/1‑150 (MQ=255)
                                                                                                                                  cAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCACGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACa                 >  2:127648/1‑150 (MQ=255)
                                                                                                                                   aaCCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAg                <  1:512860/150‑1 (MQ=255)
                                                                                                                                    aCCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGc               <  1:1141843/150‑1 (MQ=255)
                                                                                                                                         tatTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCTCTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCAGc          >  2:1024187/1‑150 (MQ=255)
                                                                                                                                          atTCGCTTATCGGCTTTATTTGATTAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGTTCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTTTTAAGCCCGACAGCCCAGCg         <  1:399306/150‑1 (MQ=255)
                                                                                                                                             cGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCCCCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGCCAGCCCAGCGGGt      >  2:1187640/1‑150 (MQ=255)
                                                                                                                                             cGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCAGCGGGt      <  1:583754/150‑1 (MQ=255)
                                                                                                                                             cGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACAATAACTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCACCGGGt      >  1:168016/1‑150 (MQ=255)
                                                                                                                                                 tATCGCCTTTCTTTGATGAAATGCAGGCAACCCCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGCCCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCAGCGGGTCgct  >  1:788051/1‑148 (MQ=255)
                                                                                                                                                 tATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCAGCGGGTCGat  >  2:111460/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ATTGCCCGTGCGGTAGAACAAGTCTCCACAGAAATGTTCCGCAAAGAGTACGCAGAAGTTTTTGAAGGCACAGCAGAGTGGAAGGGAATTAACGTCACACGATCCGATACCTACGGTTGGCAGGAGGACTCAACCTATATTCGCTTATCGCCTTTCTTTGATGAAATGCAGGCAACACCAGCACCAGTGGAAGATATTCACGGTGCGCGGATCCTCGCAATGCTGGGGGATTCAGTCACCACTGACCATATCTCTCCGGCGGGCAGTATTAAGCCCGACAGCCCAGCGGGTCGAT  >  NC_000913/1337610‑1337904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: