Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 439,372 C→T 6.7% G374G (GGC→GGT)  guaB → IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730439,3720CT 6.7% 243.2 / 2.4 81G374G (GGC→GGT) guaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base T (2/4);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAAGCACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCGAAG  >  CP000730/439232‑439511
                                                                                                                                            |                                                                                                                                           
aCTGAAGCACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGgc                                                                                                                                             >  1:693496/1‑141 (MQ=255)
aCTGAAGCACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGgc                                                                                                                                             >  2:312918/1‑141 (MQ=255)
  tGAAGCACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGgcgc                                                                                                                                           <  1:503160/141‑1 (MQ=255)
   gAAGCACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCa                                                                                                                                          <  2:979822/141‑1 (MQ=255)
       cACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGCCATGCGGCTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAg                                                                                                                                      >  2:425414/1‑141 (MQ=255)
        aCGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGa                                                                                                                                     >  2:579642/1‑141 (MQ=255)
        aCGCAAACATGGTAAAGCTATCATTGCTGATGGTGATAATAAATGCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCATTTACTGAAGAAATCCCAGGCNCATCAGa                                                                                                                                     >  2:813043/1‑141 (MQ=255)
         cGCAAACATGGTAACGCTATCATTGCTTATGCTGTTATTAACTGCTCAGGAGATATCATTAAAGCATTAGCTGTGGGTGTACATATGGTTATTTTAGGTAGATGATTAGCAGGAACTGAAGAAAGCCCAGGCGCAAcgcat                                                                                                                                    <  2:647019/141‑5 (MQ=255)
         cGCAAACATGGTAAAGCTATCATTGCTGATGGTGCTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGaa                                                                                                                                    <  1:1005376/141‑1 (MQ=255)
         cGCAAACATGGTAAAGATATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGaa                                                                                                                                    >  1:681230/1‑141 (MQ=255)
             aaCATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAAttt                                                                                                                                <  2:482534/141‑1 (MQ=255)
             aaCATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAAttt                                                                                                                                >  2:95719/1‑141 (MQ=255)
               cATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTc                                                                                                                              >  2:732895/1‑141 (MQ=255)
                   gTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAg                                                                                                                          >  2:328998/1‑141 (MQ=255)
                        gCTATCATTGCTGATGGGGGTATTAAATTCTGAGTAGATATCATTAAAGCATTAGCTGCTTGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGa                                                                                                                     <  1:762510/141‑1 (MQ=255)
                         cTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGTCCAGGTGCTACAGAAATTTTCCAAGGTAGAc                                                                                                                    <  2:356561/141‑1 (MQ=255)
                         cTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGAc                                                                                                                    <  2:1060464/141‑1 (MQ=255)
                           aTCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACaa                                                                                                                  <  2:307764/141‑1 (MQ=255)
                              aTTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAatat                                                                                                               >  1:356941/1‑141 (MQ=255)
                                tGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATaa                                                                                                             <  1:865260/141‑1 (MQ=255)
                                        gtggtATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATAcc                                                                                                     >  1:760511/1‑141 (MQ=255)
                                         tggtATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCg                                                                                                    <  2:1052176/141‑1 (MQ=255)
                                         tggtATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCACAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCg                                                                                                    <  2:396921/141‑1 (MQ=255)
                                           gtATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTg                                                                                                  <  1:321564/141‑1 (MQ=255)
                                             aTTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGt                                                                                                >  2:842398/1‑141 (MQ=255)
                                              ttAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGTCCAGGTGCTACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTa                                                                                               <  2:392480/141‑1 (MQ=255)
                                               tAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTAt                                                                                              >  2:741562/1‑141 (MQ=255)
                                                aaaTTCTCGGGGGATATCATTAGGGCATTGGGTGCTGGTGGACATGGGGTTATGTTAGGTGGCTTATTAGCAGGTACTGAAGAAAGCCGAGGTGCAACAGACATTTGCCAAGGTGGACAATATAAGGTATACCGTGGGATg                                                                                             <  2:242927/141‑1 (MQ=255)
                                                aaaTTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATg                                                                                             >  2:606655/1‑141 (MQ=255)
                                                 aaTTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGCACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGAGGTATgg                                                                                            >  2:665835/1‑141 (MQ=255)
                                                  aTTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATggg                                                                                           >  1:682889/1‑141 (MQ=255)
                                                      tcAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGctct                                                                                       >  1:191308/1‑141 (MQ=255)
                                                       cAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCtt                                                                                      <  2:443487/141‑1 (MQ=255)
                                                        aGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCttt                                                                                     <  2:895827/141‑1 (MQ=255)
                                                          gagaTATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTTTGGGCTCTTTAg                                                                                   >  2:952216/1‑141 (MQ=255)
                                                               atCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCg                                                                              <  1:136195/141‑1 (MQ=255)
                                                                tCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGGTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGa                                                                             <  1:851003/141‑1 (MQ=255)
                                                                 gattaaagtcgttgctgctGGTTGTCATGCGGTTTTGTTAGGTAGCTTATTAGGAGGTACTGAAGAAAGCGCAGGTGCAACAGAAATTTTCCCGGGTAGACAATATCATGTAGACCGTGGTATGGGCTCTTTAGGTGAGAt                                                                            <  2:993792/128‑1 (MQ=255)
                                                                 cATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGAt                                                                            <  1:261555/141‑1 (MQ=255)
                                                                   ttAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGGGATgg                                                                          >  2:673110/1‑141 (MQ=255)
                                                                   ttAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATgg                                                                          <  2:370015/141‑1 (MQ=255)
                                                                      aaGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGaaa                                                                       >  2:737318/1‑141 (MQ=255)
                                                                        gCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGaaaaa                                                                     <  1:451128/141‑1 (MQ=255)
                                                                        gCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGaaaaa                                                                     >  2:1020005/1‑141 (MQ=255)
                                                                          aTTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGTGATGGAAAAAgg                                                                   >  2:1071879/1‑141 (MQ=255)
                                                                          aTTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAAGATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAgg                                                                   >  1:860366/1‑141 (MQ=255)
                                                                              gctgctGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGTCCAGGTGCTACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCa                                                               >  2:209176/1‑141 (MQ=255)
                                                                               ctgctgGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGTCCAGGTGCTACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCaa                                                              >  1:415229/1‑141 (MQ=255)
                                                                                 gctgGTGGACATGCGGTTATGTTAGGTAGTTTATTAGCAGGTACTGAAGAAAGCCTAGGCGCAACAGAAATTTTCCGAGTTAGACAATATAAAGTATTCCGTGGTATGGGCTCTTTAGGTGCGATAGAAAAAGGTTCAAAc                                                            <  2:291931/141‑1 (MQ=255)
                                                                                 gctgGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAAc                                                            >  1:697720/1‑141 (MQ=255)
                                                                                    ggtggACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAGCGAc                                                         <  2:63151/141‑1 (MQ=255)
                                                                                     gtggACATGCGGTTATGTTAGGTAGTTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGAcc                                                        <  2:558819/141‑1 (MQ=255)
                                                                                       ggACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGt                                                      <  2:1086767/141‑1 (MQ=255)
                                                                                       ggACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGt                                                      <  1:328998/141‑1 (MQ=255)
                                                                                         aCATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTa                                                    <  2:465701/141‑1 (MQ=255)
                                                                                             gCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTc                                                >  2:1427/1‑141 (MQ=255)
                                                                                              cGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGACATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTcc                                               >  1:566055/1‑141 (MQ=255)
                                                                                              cGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTcc                                               >  2:494436/1‑141 (MQ=255)
                                                                                              cGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTcc                                               >  1:936345/1‑141 (MQ=255)
                                                                                              cGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTCCGATGGAAAAAGGTTCAAACGACCGTTACTTcc                                               <  1:563806/141‑1 (MQ=255)
                                                                                                  tATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAaga                                           >  2:796370/1‑141 (MQ=255)
                                                                                                      ttAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGAc                                       >  2:413142/1‑141 (MQ=255)
                                                                                                        aGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACaa                                     <  1:755058/141‑1 (MQ=255)
                                                                                                        aGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACaa                                     >  2:277861/1‑141 (MQ=255)
                                                                                                            agaatattaGCAGGTACTGAAGAAAGCCCAGGCGAAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAgcg                                 <  2:776852/137‑1 (MQ=255)
                                                                                                              cTTATTAGCAGGTACTGAAGAAAGCCGAGGCGCAGCAGAAATTTTACAAGGTAGACAATATAAAGTATACCGTGGTATGGTCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGcc                               >  2:107278/1‑141 (MQ=255)
                                                                                                              cTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGcc                               >  2:913150/1‑141 (MQ=255)
                                                                                                                  gtaGCAGCGACTGAAGAAAGCCCAGGCGCAACGGAAGTTTTCCAAGGTATACAAGAGAAAGTATACCGTGGTATGGGCTCTTGAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAg                           <  2:167142/140‑1 (MQ=255)
                                                                                                                   taGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGa                          <  1:1035660/141‑1 (MQ=255)
                                                                                                                   taGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGa                          <  1:296718/141‑1 (MQ=255)
                                                                                                                     gCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACATCCAAGAGGACAAAGCGCCTAAGaaa                        >  2:860575/1‑141 (MQ=255)
                                                                                                                     gCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGAGATGAAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGaaa                        <  2:681230/141‑1 (MQ=255)
                                                                                                                      cAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAAt                       >  1:513428/1‑141 (MQ=255)
                                                                                                                        ggTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAAttt                     <  1:532626/141‑1 (MQ=255)
                                                                                                                            cTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTc                 >  2:294277/1‑141 (MQ=255)
                                                                                                                              gaagaaAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCt               >  2:924370/1‑141 (MQ=255)
                                                                                                                                agaaAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGa             >  2:14231/1‑141 (MQ=255)
                                                                                                                                agaaAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTAGTTCCTGa             >  2:856155/1‑141 (MQ=255)
                                                                                                                                  acaGCCCAGTCGCAACCGATATTGTCCATGGTAGACAATATACAGTATACCGTTGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGGAATTTGTTCCTGAAg           <  2:486642/139‑1 (MQ=255)
                                                                                                                                    aGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGt         >  2:986507/1‑141 (MQ=255)
                                                                                                                                    aGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGt         <  1:930111/141‑1 (MQ=255)
                                                                                                                                        cAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCg     >  1:127988/1‑141 (MQ=255)
                                                                                                                                        cAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCg     >  2:299313/1‑141 (MQ=255)
                                                                                                                                          ggCGCAACAAAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAACCGCCTAAGAAATTTGTTCCTGAAGGTATCgca   >  2:872976/1‑139 (MQ=255)
                                                                                                                                           gcgcAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCGAAg  <  1:338879/141‑1 (MQ=255)
                                                                                                                                           gcgcAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCGAAg  >  2:528565/1‑141 (MQ=255)
                                                                                                                                           gcgcAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCGAAg  >  2:369505/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                           
ACTGAAGCACGCAAACATGGTAAAGCTATCATTGCTGATGGTGGTATTAAATTCTCAGGAGATATCATTAAAGCATTAGCTGCTGGTGGACATGCGGTTATGTTAGGTAGCTTATTAGCAGGTACTGAAGAAAGCCCAGGCGCAACAGAAATTTTCCAAGGTAGACAATATAAAGTATACCGTGGTATGGGCTCTTTAGGTGCGATGGAAAAAGGTTCAAACGACCGTTACTTCCAAGAGGACAAAGCGCCTAAGAAATTTGTTCCTGAAGGTATCGAAG  >  CP000730/439232‑439511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑