Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,083,827 G→C R52P (CGC→CCC)  hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,083,8270GC98.7% 246.7 / NA 76R52P (CGC→CCC) hisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (38/37);  minor base T (0/1);  total (38/38)
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.97e-01

AACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCGTCGTCTGGAT  >  CP009273/2083682‑2083972
                                                                                                                                                 |                                                                                                                                                 
aacaCTCGTTTACGCATAGCTATGCAGAACTCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCccgt                                                                                                                                                <  2:770448/149‑1 (MQ=255)
  cgcccGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCccgtgc                                                                                                                                              <  2:1080862/145‑1 (MQ=255)
   acTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCccgtgcg                                                                                                                                             >  2:823587/1‑149 (MQ=255)
    cTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCccgtgcgt                                                                                                                                            <  2:480933/149‑1 (MQ=255)
        tttACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTgacg                                                                                                                                        >  2:520873/1‑149 (MQ=255)
         tcaCGCAGAGCTATGCACAAATCCGGCCCTTTAAGTGATCACTCACGCGCATTGCTGGCGGCCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTgacga                                                                                                                                       <  2:709798/147‑1 (MQ=255)
            cGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTgacgacga                                                                                                                                    <  2:1403607/149‑1 (MQ=255)
             gCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGCTCGCGATCGCAGACAACATGCCGATTGATATTCGTCCCGTGCGTgacgacgac                                                                                                                                   >  1:424992/1‑149 (MQ=255)
              cATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACa                                                                                                                                  >  2:1032004/1‑149 (MQ=255)
                 aGCTAGGCAGAACTCCGGCCGTTTAAGTGATGACTCACGCGACTTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATACTGCCCGTGCGTGACGACGACATTc                                                                                                                               <  1:728665/149‑1 (MQ=255)
                  gCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAAGTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTcc                                                                                                                              <  2:75650/149‑1 (MQ=255)
                    tATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCg                                                                                                                            <  1:1312401/149‑1 (MQ=255)
                     aTGCAGAAATCCGGCCGTGTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGACGATTGATCTACTGCCCATGCGTgactacaaggttaccga                                                                                                                           >  2:1000506/1‑138 (MQ=255)
                         aGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTg                                                                                                                       <  2:85537/149‑1 (MQ=255)
                         aGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTg                                                                                                                       >  2:1042809/1‑149 (MQ=255)
                         aGAAATCCGGCCGTGTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTg                                                                                                                       <  2:1033521/149‑1 (MQ=255)
                          ggaaTCCGGCCGTTTAAGTGATGACTCACGCGAATTGTTGGCGCGCTGGGGGAGTAAAATTAATCTTCACACCCAGCGTCTGATCGCGATGGCAGAAAACAGGCCGATTGATATTCGGCCCGTGCGTGACGACGACATTCCCGGTCTgg                                                                                                                      <  2:718792/147‑1 (MQ=255)
                          gAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTgg                                                                                                                      >  1:1370109/1‑149 (MQ=255)
                                cGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAatgg                                                                                                                <  1:1003888/149‑1 (MQ=255)
                                 ggCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAatgga                                                                                                               <  2:129908/149‑1 (MQ=255)
                                 ggCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAatgga                                                                                                               >  2:619389/1‑149 (MQ=255)
                                  gCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGATGACATTCCCGGTCTGGTAatggat                                                                                                              <  1:1082488/149‑1 (MQ=255)
                                   ccGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCTGATTGATATTCTGCCAATAGTTTACGACGACATTCCCGGTCTGGTAatggatg                                                                                                             <  2:677319/149‑1 (MQ=255)
                                       ttAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGt                                                                                                         <  2:358405/149‑1 (MQ=255)
                                        tAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTg                                                                                                        <  1:1253297/149‑1 (MQ=255)
                                        tAAGTGATGACTCACGCGAAGTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTg                                                                                                        <  1:1099414/149‑1 (MQ=255)
                                         aaGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTgg                                                                                                       <  2:94271/149‑1 (MQ=255)
                                              atgaCTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGAc                                                                                                  >  1:199256/1‑149 (MQ=255)
                                                    cACGCGAATGGCTGGCGCGCTGTGGCATTAAAATTAACCTTCCCACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTggg                                                                                            <  2:48826/149‑1 (MQ=255)
                                                     aCGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGa                                                                                           >  1:750592/1‑149 (MQ=255)
                                                     aCGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGa                                                                                           >  2:734058/1‑149 (MQ=255)
                                                     aCGCGAATTGCTGGCGCGCTATGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGCCGACGACATTCCCGGTCTGGGAATGGATGGCGTGGTAGACCTTGGGa                                                                                           >  1:553271/1‑149 (MQ=255)
                                                      cgcgAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTggggt                                                                                          >  1:128741/1‑147 (MQ=255)
                                                       gcgAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGAtt                                                                                         >  2:200627/1‑149 (MQ=255)
                                                        cgAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCATGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTa                                                                                        >  2:949815/1‑149 (MQ=255)
                                                        cgAATTGCCGGCACGCTGGGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATAACAAACCAGAACCTGAAGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTa                                                                                        <  2:956987/149‑1 (MQ=255)
                                                          aaTTGCTGGCGCGCGGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATc                                                                                      <  2:631712/149‑1 (MQ=255)
                                                          aaTTGCGCGCGGGCTGTGGCAATAAAATTAAACTTCCGACCCCGCGCCTGGTCCCGATGGCAGAAAACATGCCGATTGATATTCTGCCTGTGCGTGATGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATc                                                                                      <  2:163772/149‑1 (MQ=255)
                                                            ttGCTGGCGCGCTGTGGCATTAAAATTAAACTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGCGTGCTCTTCTGCCCGTGCGTGCCGACGGCGTTCCCGGTCTGGTAATGGATGGCGTGGTAGACTCTGGGATTaccag                                                                                    >  1:54073/1‑145 (MQ=255)
                                                                tGGGGCGCTGTGGCAGTAGAAGTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGaa                                                                                >  2:713853/1‑149 (MQ=255)
                                                                   cgcgcTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAAc                                                                             <  1:1298747/149‑1 (MQ=255)
                                                                     cgcTGGGGCATTAAAACTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCTCGTGCGTGACGACGTCATTCTCTGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGt                                                                           <  1:1000506/149‑1 (MQ=255)
                                                                      gctagggCATTAAAATTAATCTTCCCACCCAGCGCCTGAACGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTg                                                                          <  1:344737/144‑1 (MQ=255)
                                                                       cTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGc                                                                         >  1:71332/1‑149 (MQ=255)
                                                                        tgtgGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCt                                                                        >  1:253807/1‑149 (MQ=255)
                                                                               ttAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGgaa                                                                    >  2:1008440/1‑146 (MQ=255)
                                                                               ttAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGgaa                                                                    <  1:1008440/146‑1 (MQ=255)
                                                                                  aaaTTAATCTTCCGACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAagag                                                              <  2:1338013/149‑1 (MQ=255)
                                                                                    aTTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCGTGGGATTATCGGCGAAAACGTGCggggcggagggct                                                            >  2:790140/1‑139 (MQ=255)
                                                                                     ttAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTggg                                                                                            <  2:112377/116‑1 (MQ=255)
                                                                                     ttAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTggg                                                                                            >  1:112377/1‑116 (MQ=255)
                                                                                        aTCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCtt                                                        >  1:623403/1‑149 (MQ=255)
                                                                                           ttCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTgacgac                                                                                                                                      >  1:1117792/1‑68 (MQ=255)
                                                                                           ttCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTgacgac                                                                                                                                      <  2:1117792/68‑1 (MQ=255)
                                                                                           ttCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAACCGTGCTGGAAGGAGAGCTGCGTAAc                                                     >  1:49813/1‑149 (MQ=255)
                                                                                                acCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAAccgccg                                                >  1:953641/1‑149 (MQ=255)
                                                                                                  ccAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCcgcg                                              >  1:317967/1‑149 (MQ=255)
                                                                                                  ccAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCcgcg                                              >  2:1140577/1‑149 (MQ=255)
                                                                                                      cgcCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTACAACCTCGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCa                                          <  1:1001882/149‑1 (MQ=255)
                                                                                                        ccTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAgg                                        >  1:122524/1‑149 (MQ=255)
                                                                                                         cTGATCGCGATGGCAGAAAACCTGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAAGGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGccagcaa                                       <  2:54073/149‑6 (MQ=255)
                                                                                                           gATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGCTTCTCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTg                                     <  1:775257/149‑1 (MQ=255)
                                                                                                           gATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTg                                     >  1:197468/1‑149 (MQ=255)
                                                                                                                 gATGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTACTGGGTAGAGTGCTAGACCTTGTGATTATCGGCGAGAACGTGCTGGAAGAAGAGCTGCTTAACAGACCCGGCCAGGCTGAAGATc                               >  2:1241116/1‑149 (MQ=255)
                                                                                                                   tGGCAGAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCa                             <  2:128741/149‑1 (MQ=255)
                                                                                                                       aGAAAACATGCCGATTGATATTCTGCCCGTGCGTGCCGACGACATTCCCGGTCTGGTCATGGATGGCGTGGTAGACCTTGGGACTATCGGCGAGAACGTGCTGGCCGACGTGCTGCCTAACCGCCGCGCCCAGGGTGTAGATCCACGca                         >  2:204948/1‑148 (MQ=255)
                                                                                                                        gAAAACATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTa                        >  1:1306834/1‑149 (MQ=255)
                                                                                                                         ataaCATGCCGATGGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTAc                       <  2:284569/147‑1 (MQ=255)
                                                                                                                           aaCATGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATAGATGGCGTGGTAGACCTTGGGCTTATCGGCGAAAACGTGCTGGAAGAAGAGCGGCTTAACCGCCGCGCCCAGGGGGAAGATCCCCGCTACtt                     >  2:187685/1‑149 (MQ=255)
                                                                                                                              aTGCCGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTAc                  >  2:1236591/1‑149 (MQ=255)
                                                                                                                                 ccGATTGATATTCTGCCCGTGCGTGACGACGACATTCCCTGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGCAGATCCACGCTACTTTACCCt               <  1:1029528/149‑1 (MQ=255)
                                                                                                                                      tGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGcgtc          <  1:200627/149‑1 (MQ=255)
                                                                                                                                      tGATATTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCCAGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGcgtc          >  1:341274/1‑149 (MQ=255)
                                                                                                                                         tatTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCgtcgtc       <  2:199256/149‑1 (MQ=255)
                                                                                                                                         tatTCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCgtcgtc       >  2:620189/1‑149 (MQ=255)
                                                                                                                                           ttCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCGTCGTCTg     <  2:953641/149‑1 (MQ=255)
                                                                                                                                            tCTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGGAGAGCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGCGAAGATCCACGCTCCTTTAACCTGCGTCGTCTgg    >  2:7693/1‑149 (MQ=255)
                                                                                                                                             cTGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCGTCGTCTGGa   <  2:623403/149‑1 (MQ=255)
                                                                                                                                              tGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGACAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTTCGTCGTCTgggt  >  2:24524/1‑147 (MQ=255)
                                                                                                                                              tGCCCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTCAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCGTCGTCTGGAt  <  2:605293/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGCGCGCTGTGGCATTAAAATTAATCTTCACACCCAGCGCCTGATCGCGATGGCAGAAAACATGCCGATTGATATTCTGCGCGTGCGTGACGACGACATTCCCGGTCTGGTAATGGATGGCGTGGTAGACCTTGGGATTATCGGCGAAAACGTGCTGGAAGAAGAGCTGCTTAACCGCCGCGCCCAGGGTGAAGATCCACGCTACTTTACCCTGCGTCGTCTGGAT  >  CP009273/2083682‑2083972

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