Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 984,870 T→A 8.2% N315Y (AAT→TAT)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273984,8700TA 8.2% 899.5 / 51.6 257N315Y (AAT→TAT) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (104/132);  new base A (15/6);  total (119/138)
Fisher's exact test for biased strand distribution p-value = 2.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 261 total aligned reads displayed.

TTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGA  >  CP009273/984721‑985019
                                                                                                                                                     |                                                                                                                                                     
 ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTg                                                                                                                                                      >  2:1629554/1‑150 (MQ=255)
  tACCGTTACACTCTACCAACTTAATGACAATATTCGGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGt                                                                                                                                                     >  2:1890930/1‑150 (MQ=255)
    ccGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCa                                                                                                                                                   <  1:4108647/150‑1 (MQ=255)
     cGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCag                                                                                                                                                  <  2:3080754/150‑1 (MQ=255)
      gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga                                                                                                                                                 <  1:1848543/150‑1 (MQ=255)
       ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCagag                                                                                                                                                >  1:3341122/1‑150 (MQ=255)
       ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCagag                                                                                                                                                >  1:2800580/1‑150 (MQ=255)
        tACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAgaga                                                                                                                                               >  1:1430275/1‑150 (MQ=255)
         acacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGaa                                                                                                                                              >  1:1589858/1‑150 (MQ=255)
          cacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGaaa                                                                                                                                             <  1:85898/150‑1 (MQ=255)
           acTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGaaaa                                                                                                                                            <  2:2089184/150‑1 (MQ=255)
            ctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAAc                                                                                                                                           >  1:4358743/1‑150 (MQ=255)
              ctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCa                                                                                                                                         <  2:1961436/150‑1 (MQ=255)
              ctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCa                                                                                                                                         >  1:2356931/1‑150 (MQ=255)
                aCCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGc                                                                                                                                       >  1:1730802/1‑150 (MQ=255)
                aCCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGc                                                                                                                                       >  1:4071815/1‑150 (MQ=255)
                   aaCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGtt                                                                                                                                    <  2:3676219/150‑1 (MQ=255)
                    aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGttt                                                                                                                                   <  1:84963/150‑1 (MQ=255)
                    aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGttt                                                                                                                                   <  1:1247073/150‑1 (MQ=255)
                       tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTAc                                                                                                                                >  2:3924045/1‑150 (MQ=255)
                       tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTAc                                                                                                                                >  2:2552258/1‑150 (MQ=255)
                         aTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTActc                                                                                                                              >  1:2905043/1‑150 (MQ=255)
                              aaTATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTgg                                                                                                                         >  1:1700392/1‑150 (MQ=255)
                              aaTATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTgg                                                                                                                         <  1:3341689/150‑1 (MQ=255)
                               atatTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGa                                                                                                                        <  2:4861971/150‑1 (MQ=255)
                                 atTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATc                                                                                                                      >  1:3215830/1‑150 (MQ=255)
                                   tCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCaa                                                                                                                    <  2:1704414/150‑1 (MQ=255)
                                    cAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAAt                                                                                                                   >  2:2425162/1‑150 (MQ=255)
                                      ggggTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTc                                                                                                                 >  2:4684887/1‑150 (MQ=255)
                                      ggggTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTc                                                                                                                 >  2:50632/1‑150 (MQ=255)
                                       gggTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTcc                                                                                                                <  1:3090200/150‑1 (MQ=255)
                                       gggTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTcc                                                                                                                >  1:1858249/1‑150 (MQ=255)
                                       gggTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTcc                                                                                                                <  1:639276/150‑1 (MQ=255)
                                        ggTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTccc                                                                                                               <  1:2234579/150‑1 (MQ=255)
                                          ttttACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAg                                                                                                             <  1:2988297/150‑1 (MQ=255)
                                          ttttACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAg                                                                                                             <  1:555629/150‑1 (MQ=255)
                                          ttttACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAg                                                                                                             <  2:379509/150‑1 (MQ=255)
                                            ttACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCt                                                                                                           <  2:3439307/150‑1 (MQ=255)
                                              aCCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCttt                                                                                                         <  2:3117557/150‑1 (MQ=255)
                                              aCCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCttt                                                                                                         >  2:1254856/1‑150 (MQ=255)
                                               ccTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAAGTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCtttt                                                                                                        <  1:1162831/150‑1 (MQ=255)
                                                 tGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCa                                                                                                      >  1:2329137/1‑150 (MQ=255)
                                                  gggggATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCAt                                                                                                     >  2:3385018/1‑150 (MQ=255)
                                                  gggggATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCAt                                                                                                     <  1:664493/150‑1 (MQ=255)
                                                    gggATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATaa                                                                                                   >  2:4317230/1‑150 (MQ=255)
                                                      gATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAAtg                                                                                                 >  2:3753266/1‑150 (MQ=255)
                                                           cGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCaa                                                                                            >  1:364122/1‑150 (MQ=255)
                                                              aGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATgg                                                                                         >  2:148803/1‑150 (MQ=255)
                                                              aGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATgg                                                                                         <  2:2354752/150‑1 (MQ=255)
                                                                    ggCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCttttt                                                                                   <  2:171550/150‑1 (MQ=255)
                                                                    ggCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCttttt                                                                                   <  2:3341122/150‑1 (MQ=255)
                                                                     gCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTc                                                                                  <  2:703667/150‑1 (MQ=255)
                                                                       gAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCAc                                                                                <  1:4825725/150‑1 (MQ=255)
                                                                          tCCCAATACAAAAACTTAATTGCACGCTGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCATATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCAcccc                                                                             <  1:3807448/150‑1 (MQ=255)
                                                                              aaTACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTc                                                                         >  1:2429527/1‑150 (MQ=255)
                                                                               aTACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCa                                                                        <  2:4090390/150‑1 (MQ=255)
                                                                                 aCCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAAc                                                                      >  2:3252874/1‑150 (MQ=255)
                                                                                  ccAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAAcc                                                                     <  1:2935179/150‑1 (MQ=255)
                                                                                   cAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCg                                                                    >  1:4390080/1‑150 (MQ=255)
                                                                                   cAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCg                                                                    >  1:4321225/1‑150 (MQ=255)
                                                                                       aCTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTc                                                                <  2:4574073/150‑1 (MQ=255)
                                                                                        cTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCg                                                               >  2:1675662/1‑150 (MQ=255)
                                                                                        cTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCg                                                               >  2:1674768/1‑150 (MQ=255)
                                                                                         ttAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGc                                                              >  1:4598272/1‑150 (MQ=255)
                                                                                          tAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGcc                                                             <  2:132210/150‑1 (MQ=255)
                                                                                           aaTTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCa                                                            <  1:948347/150‑1 (MQ=255)
                                                                                              tGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCcgag                                                         >  2:582289/1‑150 (MQ=255)
                                                                                                cACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTc                                                       <  1:3893527/150‑1 (MQ=255)
                                                                                                cACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTc                                                       >  2:4582638/1‑150 (MQ=255)
                                                                                                cACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTc                                                       <  2:1065793/150‑1 (MQ=255)
                                                                                                 aCGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCa                                                      <  1:2162307/150‑1 (MQ=255)
                                                                                                  cgcgGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCAt                                                     <  2:2796317/150‑1 (MQ=255)
                                                                                                        aaGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCa                                               >  1:1312597/1‑150 (MQ=255)
                                                                                                         agagAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAg                                              <  1:4425404/150‑1 (MQ=255)
                                                                                                             aaGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGccc                                          >  1:4232791/1‑150 (MQ=255)
                                                                                                             aaGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGccc                                          >  1:4232563/1‑150 (MQ=255)
                                                                                                             aaGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGccc                                          >  1:900399/1‑150 (MQ=255)
                                                                                                              aGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                                         <  1:2411602/150‑1 (MQ=255)
                                                                                                              aGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                                         >  2:2294869/1‑150 (MQ=255)
                                                                                                              aGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCt                                         >  2:4746442/1‑150 (MQ=255)
                                                                                                               gTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTg                                        <  2:1430275/150‑1 (MQ=255)
                                                                                                                 ggCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGAt                                      <  2:1609168/150‑1 (MQ=255)
                                                                                                                 ggCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGAt                                      <  1:1317617/150‑1 (MQ=255)
                                                                                                                 ggCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGAt                                      <  2:4821358/150‑1 (MQ=255)
                                                                                                                    ggTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATAcc                                   <  2:1799075/150‑1 (MQ=255)
                                                                                                                    ggTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATAcc                                   <  1:1033295/150‑1 (MQ=255)
                                                                                                                        taGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACt                               >  2:4267086/1‑150 (MQ=255)
                                                                                                                        taGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACt                               <  1:3105459/150‑1 (MQ=255)
                                                                                                                         aGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTa                              <  2:3962896/150‑1 (MQ=255)
                                                                                                                         aGCTCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTa                              >  2:2848650/1‑150 (MQ=255)
                                                                                                                            tCAACCGCTTTGCGTAAATCCAGATAGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCg                           >  2:1478029/1‑150 (MQ=255)
                                                                                                                              aaCCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGaa                         >  1:467493/1‑150 (MQ=255)
                                                                                                                                ccGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACt                       >  1:4810973/1‑150 (MQ=255)
                                                                                                                                  gCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATCATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCg                     <  1:3275756/150‑1 (MQ=255)
                                                                                                                                  gCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCg                     >  1:1669028/1‑150 (MQ=255)
                                                                                                                                  gCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCg                     <  1:2561788/150‑1 (MQ=255)
                                                                                                                                   cTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGa                    <  2:2023491/150‑1 (MQ=255)
                                                                                                                                     ttGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGaca                  <  2:1893826/150‑1 (MQ=255)
                                                                                                                                            aaTCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAAt           >  2:3177790/1‑150 (MQ=255)
                                                                                                                                                 aGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGc      >  1:2871288/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAGAAAACCAGCGTTTTACTCTGTGGATCAATTCCCAGCTTTTCATAATGTGCAATGGCTTTTTCACCCCATTCAACCGGGTCGCCAGAGTCATGACGCAGGCCCTGATACCGACTAGCGAACTCGACACCGAAATCACGCAGGA  >  CP009273/984721‑985019

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