Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 847,214 T→A 12.3% intergenic (‑95/+44) glnP ← / ← glnH glutamine transporter subunit/glutamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913847,2140TA12.3% 198.2 / 18.6 90intergenic (‑95/+44)glnP/glnHglutamine transporter subunit/glutamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base A (4/7);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

TCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGC  >  NC_000913/847070‑847361
                                                                                                                                                |                                                                                                                                                   
tCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAgggg                                                                                                                                                <  2:875184/150‑1 (MQ=255)
  aaTCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCCCCCTAATCCTCTCCCCATAGGGgag                                                                                                                                              >  2:729526/1‑150 (MQ=255)
  aaTAAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGgag                                                                                                                                              <  2:985159/150‑1 (MQ=255)
    tCAGAAGCGGAATGGCAGGCCAGATGTCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGagag                                                                                                                                            <  2:1284626/150‑1 (MQ=255)
    tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCTTATGTTGTTCCTGTTACCGGGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTGTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGagag                                                                                                                                            <  2:953099/150‑1 (MQ=255)
    tCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGagag                                                                                                                                            <  1:1221624/150‑1 (MQ=255)
     cAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgg                                                                                                                                           <  1:1133407/150‑1 (MQ=255)
      aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAggg                                                                                                                                          >  1:698563/1‑150 (MQ=255)
      aGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAggg                                                                                                                                          >  2:893809/1‑150 (MQ=255)
      aGAAGCGGAACGGCAGGCCAGATGGCACACCAGTCAAACTGCTTATGTTGTTCCTGTTACCTGGGTGTAAATAATCAAATTCATGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgga                                                                                                                                          <  1:42877/150‑2 (MQ=255)
        aaGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAggatt                                                                                                                                        <  1:931615/150‑4 (MQ=255)
               aTGGCAGGCCAGATGGCACTCCAGTCAACCTGCATAGGTTGTTCCTGTTACCGGGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAg                                                                                                                                 <  1:122554/150‑1 (MQ=255)
               aTGGCAGGCCAGAAGGCACTCCAGTCAAACTGCATATGTGGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAg                                                                                                                                 <  1:603709/150‑1 (MQ=255)
                 ggCATGCCAGATGCCACTCCATTCAAACTGCATATGTTGTTCCTGTTCCCGTGGTGTATATAATCAAAATCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCACCATAGGGGAGAGGGAACTGCCAGTg                                                                                                                               <  1:673898/150‑1 (MQ=255)
                  gCAGGCCAGATGCCACTCCAGTCAAACTGCATATCTTGGTCCTGTCACCGCGGTGCAAATAATCAAATTCAGCGTAATAAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTACTCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGc                                                                                                                              <  1:222168/150‑1 (MQ=255)
                        cAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTAc                                                                                                                       <  2:982926/150‑1 (MQ=255)
                          gATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTCGTGTAAATATTAAAATTCTGGGTAATTAATGGCCCCGGAGCTGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACa                                                                                                                      <  2:214482/150‑1 (MQ=255)
                            tGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGCTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAgg                                                                                                                    <  2:14167/150‑1 (MQ=255)
                              gCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGgtg                                                                                                                  <  2:1133622/150‑1 (MQ=255)
                               cACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTTATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGTAACTGACAGTGCGTTTTACAGgggt                                                                                                                 >  2:226381/1‑150 (MQ=255)
                                    cAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGt                                                                                                            >  2:641266/1‑150 (MQ=255)
                                        cAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGttatt                                                                                                        >  2:124596/1‑150 (MQ=255)
                                               cATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCgg                                                                                                 <  1:967780/150‑1 (MQ=255)
                                                 tatGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGtt                                                                                               >  1:254217/1‑150 (MQ=255)
                                                  atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGtta                                                                                              >  1:565852/1‑149 (MQ=255)
                                                  atGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTc                                                                                              <  1:62087/150‑1 (MQ=255)
                                                   tgttgtTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCa                                                                                             >  1:851780/1‑150 (MQ=255)
                                                   tgttgtTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGAGCCCCTCACCCTAATCCTCTCCCCATAGGGGTGAGGGAACTGCCAGTGCGTTTTTCAGGTGTAGCGTTATTATTTCGGTTCa                                                                                             >  2:1267587/1‑150 (MQ=255)
                                                           cTGTTACCGGGGTGTAAATAATCAAATTCAGGGTACTTAAGGGCCCCGGCGCGGGTGTGCGCCCCTCTCCCTAATCCTCTCCCCAGAGGGGAGAGGGTACTGGCTGGGCGTTTTACAGGTGTAGGGTTTTTATTTCGGTTCAggatcgca                                                                                     >  1:811080/1‑143 (MQ=255)
                                                            tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgg                                                                                                                                           >  2:604036/1‑95 (MQ=255)
                                                            tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCTATGGGGAGAgg                                                                                                                                           <  1:604036/95‑1 (MQ=255)
                                                            tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAc                                                                                    <  1:14942/150‑1 (MQ=255)
                                                            tGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGACCCTAACACTAATCCTCTCACCATAGGGGAGAGGGAACTGCCAGTGCGGTTTACATGTGTAGAGTTATTATTTCGTTTCAGTACCTAAc                                                                                    >  2:1118660/1‑150 (MQ=255)
                                                            tGTTACCGTGGTGTAAATAATCAAAATCAGCGGAATAAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAAc                                                                                    <  1:848268/150‑1 (MQ=255)
                                                                  cGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGACCTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCAtttt                                                                              >  2:113820/1‑150 (MQ=255)
                                                                     ggtgTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTg                                                                           <  1:8570/150‑1 (MQ=255)
                                                                      gtgtAAAAAATCACATTCATGGTAATTAATGGCCCCGGAGCGGGGTTGCGCCACTCACCCTTATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTTCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGt                                                                          <  1:245425/150‑1 (MQ=255)
                                                                      gtgtAAAAAATCAAATTCAGTTTAATTAATGGCCCCGGAGCGGTTTTTCGCACATCATACTAATCCTCTCCCCATTGGGGAGAGGGAACTGCCATTTCGTTTTACATCTGTAGCGTTATTATTTCGGTTCAGTACCTAACCATTTTTTGt                                                                          <  1:931253/150‑1 (MQ=255)
                                                                       tgtAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATGTTTTGTa                                                                         >  1:798772/1‑150 (MQ=255)
                                                                            ataatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCAATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGAttt                                                                    >  2:957878/1‑150 (MQ=255)
                                                                             taatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGAtttg                                                                   >  1:494267/1‑149 (MQ=255)
                                                                             taatCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTc                                                                   <  2:922840/150‑1 (MQ=255)
                                                                               atCAAATTCAGGGTAATTACTGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                 <  2:186766/150‑1 (MQ=255)
                                                                               atCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGt                                                                 <  1:720666/150‑1 (MQ=255)
                                                                                tCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGtt                                                                >  2:175493/1‑150 (MQ=255)
                                                                                   aaTTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTa                                                             <  1:532558/150‑1 (MQ=255)
                                                                                     ttCAGGGTAATTAATTGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAg                                                           <  1:394924/150‑1 (MQ=255)
                                                                                       cAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGtt                                                         >  1:575314/1‑150 (MQ=255)
                                                                                         gggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTcc                                                       <  1:1220709/150‑1 (MQ=255)
                                                                                          ggTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGGTATTATTTCGGTTCAGTACCGCACCATTTTTTGTAGATTTCGTTGTAAGTTCCg                                                      >  2:388396/1‑150 (MQ=255)
                                                                                              attaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtct                                                  >  1:515295/1‑150 (MQ=255)
                                                                                               ttaatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTtctc                                                 <  2:376283/150‑1 (MQ=255)
                                                                                                 aatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGCACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTcgc                                               >  2:1308747/1‑150 (MQ=255)
                                                                                                 aatGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTcgc                                               >  2:621311/1‑150 (MQ=255)
                                                                                                   tGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCgcgc                                             >  2:944576/1‑150 (MQ=255)
                                                                                                    ggCCCCGGAGCGGGTTTGCGCCCCTCCCCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGGTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCa                                            >  2:916087/1‑150 (MQ=255)
                                                                                                    ggCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGGTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCa                                            <  1:93898/150‑1 (MQ=255)
                                                                                                       cccGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAggg                                         >  2:1102398/1‑150 (MQ=255)
                                                                                                                     gcgcCCCTCACCCTAATCCTCACCCCATAGGGGAGAGGGAACTGCCAGGGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCg                           <  1:768650/150‑1 (MQ=255)
                                                                                                                      cgcCCCTCACCGTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGt                          >  1:289730/1‑150 (MQ=255)
                                                                                                                      cgcCCCTCACCCTAATCCTCTCCCCATTGGGGAGAGGGAACTGCCATTGCGTTTTACAGGTGTAGCGATATTATTTCGGTTCAGTACCGAACCATTTTATGTAGCTTTCGTTGTAAGATCCGTTATCGCGCAGGGTTTTCAActctccct                          >  2:1001024/1‑143 (MQ=255)
                                                                                                                       gcCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGtt                         >  1:274934/1‑150 (MQ=255)
                                                                                                                       gcCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGtt                         >  1:925775/1‑150 (MQ=255)
                                                                                                                        ccccTCACCCTAATCCTCTCCCCATAGGGTAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATATTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGttt                        >  1:388340/1‑149 (MQ=255)
                                                                                                                        ccccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                        <  1:1138099/150‑1 (MQ=255)
                                                                                                                        ccccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                        >  1:25661/1‑150 (MQ=255)
                                                                                                                        ccccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGGACCGAACCATTTTTTGTAGATTTCTTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTg                        >  2:769326/1‑150 (MQ=255)
                                                                                                                        ccccTCACACTAATCCTCTCCCCATAGGGGAGAGGGCCCTGCCAGTGCGTGTTACAGGTGTAGCGTTATTATTTCGGGTCAGTACCGATCCATTTTTTGTAGATTTCGTTGTAAGTTCAGTTCTCGCGAAGTGTGTTCTACGCGCCGTTg                        >  1:393840/1‑150 (MQ=255)
                                                                                                                         cccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGa                       >  1:1203462/1‑150 (MQ=255)
                                                                                                                          ccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGAc                      <  2:402314/150‑1 (MQ=255)
                                                                                                                          ccTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCCTTGAc                      >  2:287437/1‑150 (MQ=255)
                                                                                                                           cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                     <  1:601410/150‑1 (MQ=255)
                                                                                                                           cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCAGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                     >  2:493969/1‑150 (MQ=255)
                                                                                                                           cTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACGGCCAGTGCGTTTTACAGGTGGAGCGTTATTATATCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACt                     <  2:121913/150‑1 (MQ=255)
                                                                                                                            acACCCTAATCCTCTCCCGATAGGGGAGAGGGAACTTCCAGTGCGGTTTACTGGTGTAGCGTTATTATTTCGGTTCAGTACCGTACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACtt                    <  2:1022599/149‑1 (MQ=255)
                                                                                                                              ccccTAAACCTCACCCCAAAGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                 <  1:363640/149‑1 (MQ=255)
                                                                                                                              aCCCTAATCCTCTCCCCATCGGGGAGAGGGAACTGCCCGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACACGCCGTTGACTTTg                  <  2:154290/150‑1 (MQ=255)
                                                                                                                               cccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGt                 >  2:1244136/1‑150 (MQ=255)
                                                                                                                                ccTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTc                <  1:335149/150‑1 (MQ=255)
                                                                                                                                   aaTCCTCTCCCCAGAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTCGCGTTATTATTTCGTTGCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTAGAGCAGGTTTTTCAACTCGCCGTTGACTTTGTCACg             <  2:645419/150‑1 (MQ=255)
                                                                                                                                    aTCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc            <  2:569101/150‑1 (MQ=255)
                                                                                                                                    aTCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGc            <  1:400430/150‑1 (MQ=255)
                                                                                                                                     cccTCTCCCCTATGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCa           >  2:1003746/2‑150 (MQ=255)
                                                                                                                                     cccTCTCCCCTATGGGGAGAGAGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCa           >  1:1234504/2‑150 (MQ=255)
                                                                                                                                     cccTCTCCCATAAGGGGAGAGGGAAATACCACTGCGTTTTACAGCTGTACCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTATCGCGAAGGGTTTTAAACGCGCCGTTGACTTTGTCACGCa           >  1:971536/2‑150 (MQ=255)
                                                                                                                                        gctcCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCt        <  1:1269654/149‑1 (MQ=255)
                                                                                                                                         tctcccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtc       <  2:238493/147‑1 (MQ=255)
                                                                                                                                            cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc    <  2:851780/150‑1 (MQ=255)
                                                                                                                                            cccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtc    >  1:874473/1‑150 (MQ=255)
                                                                                                                                             ccATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCtcgtcg   >  1:759718/1‑150 (MQ=255)
                                                                                                                                              cctAGGGGAGAGGGAACTGCCAGTGCGTGTTCCAGGTGTCGCGTGATTAGTTCGGTTCAGGACCGAACCATTTTCTGTAGATTTCGTAGTCAGGACCGTTCTCGCGCAGGGTGTTCAACGCGCCGGTGACTTGGTCACGCAGCTCGTCGc  <  2:839765/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCAATCAGAAGCGGAATGGCAGGCCAGATGGCACTCCAGTCAAACTGCATATGTTGTTCCTGTTACCGTGGTGTAAATAATCAAATTCAGGGTAATTAATGGCCCCGGAGCGGGTTTGCGCCCCTCACCCTAATCCTCTCCCCATAGGGGAGAGGGAACTGCCAGTGCGTTTTACAGGTGTAGCGTTATTATTTCGGTTCAGTACCGAACCATTTTTTGTAGATTTCGTTGTAAGTTCCGTTCTCGCGCAGGGTTTTCAACGCGCCGTTGACTTTGTCACGCAGCTCGTCGC  >  NC_000913/847070‑847361

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