Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 34,120 (T)9→8 intergenic (+82/‑180) carB → / → caiF carbamoyl‑phosphate synthase large subunit/cai operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091334,1120T.98.6% 328.2 / ‑0.4 73intergenic (+74/‑188)carB/caiFcarbamoyl‑phosphate synthase large subunit/cai operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (40/31);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

CCACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCAAAAC  >  NC_000913/33962‑34258
                                                                                                                                                      |                                                                                                                                                  
ccACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGt                                                                                                                                                    <  1:459707/151‑1 (MQ=255)
 cACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtt                                                                                                                                                   >  1:699716/1‑151 (MQ=255)
    cgcgATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGtt                                                                                                                                                   <  2:256517/148‑1 (MQ=255)
      cgATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttt                                                                                                                                              >  1:997930/1‑151 (MQ=255)
      cgATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttt                                                                                                                                              >  1:1114183/1‑151 (MQ=255)
          ggCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGttttttttgc                                                                                                                                           <  2:43023/150‑3 (MQ=255)
                gAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTa                                                                                                                                    >  1:1007602/1‑150 (MQ=255)
                 aaTGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAg                                                                                                                                  <  2:649445/151‑1 (MQ=255)
                    gCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTa                                                                                                                                >  1:981109/1‑150 (MQ=255)
                    gCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTa                                                                                                                                >  2:536598/1‑150 (MQ=255)
                      cGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAAt                                                                                                                             >  1:343680/1‑151 (MQ=255)
                          gCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCttt                                                                                                                         >  2:560307/1‑151 (MQ=255)
                             aCTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTgg                                                                                                                      >  1:111466/1‑151 (MQ=255)
                               tGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTCGCCCATTAAGTAAATCTTTTgggg                                                                                                                    >  1:1137001/1‑151 (MQ=255)
                                gAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTgggg                                                                                                                    >  1:603336/1‑150 (MQ=255)
                                gAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTgggg                                                                                                                    <  1:36295/150‑1 (MQ=255)
                                 aaaaaGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTTATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTggg                                                                                                                     >  2:61124/1‑148 (MQ=255)
                                  aaaaGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAAt                                                                                                                 >  1:805527/1‑151 (MQ=255)
                                   aaaGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAAt                                                                                                                 >  1:661767/1‑150 (MQ=255)
                                   aaaGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATc                                                                                                                <  1:913797/151‑1 (MQ=255)
                                       tAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGat                                                                                                             >  2:149242/1‑150 (MQ=255)
                                        aaTTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGata                                                                                                            <  2:121631/150‑1 (MQ=255)
                                         aTTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAtt                                                                                                          >  2:596759/1‑151 (MQ=255)
                                         aTTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAtt                                                                                                          >  1:1204801/1‑151 (MQ=255)
                                         aTTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAtt                                                                                                          <  2:197120/151‑1 (MQ=255)
                                           ttCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAtttt                                                                                                        >  1:60945/1‑151 (MQ=255)
                                           ttCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttt                                                                                                         <  1:310977/150‑1 (MQ=255)
                                           ttCCGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAtttt                                                                                                        >  2:221918/1‑151 (MQ=255)
                                             cGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATAttttt                                                                                                       <  2:367529/150‑1 (MQ=255)
                                                 gCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatg                                                                                                   <  1:394101/150‑1 (MQ=255)
                                                  cAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTtgatg                                                                                                   <  2:676028/149‑1 (MQ=255)
                                                    ggAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACAt                                                                                               >  1:452167/1‑151 (MQ=255)
                                                    ggAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACAt                                                                                               >  2:616167/1‑151 (MQ=255)
                                                     gAAATGCTCGCACAGATCAAATAATAGCGTGTCGTGGCAGAGATTTTGCTTCCGCTTATTTGCTAGAATAAGTTATTCGGTTCTCTGGTGAGGGCCCTTTTTTTTGCCCATTCAGTAAATTTTTGGGGGCATCGATAGTTTtggtaata                                                                                                <  2:634522/149‑7 (MQ=255)
                                                       aaTGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATa                                                                                              >  1:390195/1‑149 (MQ=255)
                                                          gcacgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAg                                                                                            >  2:907916/1‑148 (MQ=255)
                                                           cacgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAgg                                                                                        >  1:242148/1‑151 (MQ=255)
                                                            acgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAgg                                                                                        >  1:3854/1‑150 (MQ=255)
                                                             cgcacAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAt                                                                                      >  1:14520/1‑151 (MQ=255)
                                                               cgcaGATAAAATAATAGCGTGTCATGGCATATATTTTGCATCCGCTAATTTGATCGAATAACTAATCCCGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAtt                                                                                     <  2:293718/148‑1 (MQ=255)
                                                                acaGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGAtt                                                                                     <  2:985906/149‑1 (MQ=255)
                                                                    aTCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTc                                                                               >  1:1028956/1‑151 (MQ=255)
                                                                            ataGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTAt                                                                       <  1:1055418/151‑1 (MQ=255)
                                                                               gCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGa                                                                    <  2:1114933/151‑1 (MQ=255)
                                                                                   gtCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGt                                                                <  2:1072441/151‑1 (MQ=255)
                                                                                      aTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGa                                                              >  2:743276/1‑150 (MQ=255)
                                                                                      aTGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGa                                                              >  1:598801/1‑150 (MQ=255)
                                                                                        ggCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTCCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAg                                                            <  2:201947/150‑1 (MQ=255)
                                                                                          cAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTg                                                         >  1:1171896/1‑151 (MQ=255)
                                                                                          cAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTg                                                         >  1:245069/1‑151 (MQ=255)
                                                                                          cAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTg                                                         >  1:482426/1‑151 (MQ=255)
                                                                                             atatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCAt                                                      <  1:894466/151‑1 (MQ=255)
                                                                                              tatTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCAt                                                      >  2:741982/1‑150 (MQ=255)
                                                                                                 ttttCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACt                                                   >  2:348145/1‑150 (MQ=255)
                                                                                                     cATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCg                                               <  1:1112137/150‑1 (MQ=255)
                                                                                                     cATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCg                                               <  2:1027649/150‑1 (MQ=255)
                                                                                                     cATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCg                                               <  1:40947/150‑1 (MQ=255)
                                                                                                        ccGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGatat                                           <  1:100594/151‑1 (MQ=255)
                                                                                                        ccGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGata                                            <  1:1041088/150‑1 (MQ=255)
                                                                                                             aaTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTAGCGATATATCca                                      <  1:1022998/151‑1 (MQ=255)
                                                                                                              aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGGATACTATCGATATATCcac                                     <  1:536620/151‑1 (MQ=255)
                                                                                                              aTTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCcac                                     >  1:220524/1‑151 (MQ=255)
                                                                                                                  gATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAAtt                                 <  2:780680/151‑1 (MQ=255)
                                                                                                                  gATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAAt                                  <  2:732022/150‑1 (MQ=255)
                                                                                                                  gATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAAt                                  <  2:1114203/150‑1 (MQ=255)
                                                                                                                     cGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAAtttt                               >  2:941859/1‑150 (MQ=255)
                                                                                                                      gAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTa                              >  1:1109180/1‑150 (MQ=255)
                                                                                                                         tAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAatat                          >  2:596178/1‑151 (MQ=255)
                                                                                                                            cTAATACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGc                       >  2:935273/1‑151 (MQ=255)
                                                                                                                              aaTACGGTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCt                     <  2:393278/151‑1 (MQ=255)
                                                                                                                                 aCGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTg                   <  1:1079550/151‑1 (MQ=255)
                                                                                                                                 aCGGGTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTAATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGt                  >  1:800408/1‑151 (MQ=255)
                                                                                                                                    gTTCTCTGATGAGGACCGTTTTTTTTGCCCATTAAGTAATTCTTTTGGGGAATCGATATTTTTGATTACCTAAGCTGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGttt                >  2:569529/1‑150 (MQ=255)
                                                                                                                                       ctctGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAAtt            >  1:1160293/1‑151 (MQ=255)
                                                                                                                                         ctGATGAGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAAtt            <  2:588031/149‑1 (MQ=255)
                                                                                                                                            atgaGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTc       <  2:1046914/151‑1 (MQ=255)
                                                                                                                                             tgaGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTc       >  2:1091780/1‑150 (MQ=255)
                                                                                                                                              gaGGACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCa      <  1:15430/150‑1 (MQ=255)
                                                                                                                                                 gACCGTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCAAAAc  <  1:560328/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCACCGCGATGGCGCTGAATGCCGATGCGACTGAAAAAGTAATTTCGGTGCAGGAAATGCACGCACAGATCAAATAATAGCGTGTCATGGCAGATATTTTTCATCCGCTAATTTGATCGAATAACTAATACGGTTCTCTGATGAGGACCGTTTTTTTTTGCCCATTAAGTAAATCTTTTGGGGAATCGATATTTTTGATGACATAAGCAGGATTTAGCTCACACTTATCGACGGTGAAGTTGCATACTATCGATATATCCACAATTTTAATATGGCCTTGTTTAATTGCTTCAAAAC  >  NC_000913/33962‑34258

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