Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,281,132 G→T 5.4% D180Y (GAC→TAC)  narJ → molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,281,1320GT 5.4% 661.6 / 11.0 173D180Y (GAC→TAC) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base G (91/72);  new base T (5/5);  total (96/77)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

AAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTG  >  CP009273/1280986‑1281278
                                                                                                                                                  |                                                                                                                                                  
aaGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCcga                                                                                                                                                   <  2:2196888/148‑1 (MQ=255)
aaGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCcga                                                                                                                                                   <  2:831282/148‑1 (MQ=255)
   gTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCcgacga                                                                                                                                                >  1:2721499/1‑148 (MQ=255)
    tttGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGaa                                                                                                                                               >  1:1737398/1‑148 (MQ=255)
     ttGAAAGATATCGCGCCGATTCTGGCATTGCTTAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATTCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAg                                                                                                                                              <  2:323288/148‑1 (MQ=255)
     ttGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAg                                                                                                                                              <  2:2514431/148‑1 (MQ=255)
     ttGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAg                                                                                                                                              <  1:1620409/148‑1 (MQ=255)
       gAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAgcg                                                                                                                                            >  1:295154/1‑148 (MQ=255)
         aaGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAgcgcg                                                                                                                                          <  1:1571131/148‑1 (MQ=255)
         aaGATATCGCGCCGATTCTGGCAGTGCTGAGCGCGCGTCTTCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGACGACGAcgcgcg                                                                                                                                          >  2:1374015/1‑148 (MQ=255)
           gATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCTACGAAGcgcgcg                                                                                                                                        >  2:516911/1‑148 (MQ=255)
           gATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCTACGAAGcgcgcg                                                                                                                                        >  2:516918/1‑148 (MQ=255)
           gATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGcgcgcg                                                                                                                                        >  2:2286212/1‑148 (MQ=255)
           gATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGcgcgcg                                                                                                                                        >  1:1558974/1‑148 (MQ=255)
            atatCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGACAAAATTGCCGACGAAGCGCGCga                                                                                                                                       >  2:2898564/1‑148 (MQ=255)
            atatCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCga                                                                                                                                       >  2:2805173/1‑148 (MQ=255)
             tatCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgat                                                                                                                                      >  2:1173657/1‑148 (MQ=255)
               tCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatga                                                                                                                                    >  1:1378918/1‑148 (MQ=255)
                cgcgCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatgat                                                                                                                                   >  1:905137/1‑148 (MQ=255)
                 gcgcCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATa                                                                                                                                  <  1:266553/148‑1 (MQ=255)
                   gcCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcg                                                                                                                                <  2:1737398/148‑1 (MQ=255)
                    ccGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgc                                                                                                                               >  1:2151870/1‑148 (MQ=255)
                    ccGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgc                                                                                                                               >  2:1589700/1‑148 (MQ=255)
                      gATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgccg                                                                                                                             >  2:1832256/1‑148 (MQ=255)
                      gATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgccg                                                                                                                             >  2:159058/1‑148 (MQ=255)
                       aTTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgccgc                                                                                                                            <  2:3201264/148‑1 (MQ=255)
                         tCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAg                                                                                                                          >  2:2646274/1‑148 (MQ=255)
                         tCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAg                                                                                                                          >  1:321416/1‑148 (MQ=255)
                              cATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCt                                                                                                                     >  1:2643082/1‑148 (MQ=255)
                                ttGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTgg                                                                                                                   <  2:1643039/148‑1 (MQ=255)
                                 tGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGa                                                                                                                  >  2:3242460/1‑148 (MQ=255)
                                   cTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATg                                                                                                                <  2:1975202/148‑1 (MQ=255)
                                   cTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATg                                                                                                                <  1:629668/148‑1 (MQ=255)
                                   cTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATg                                                                                                                <  2:1241948/148‑1 (MQ=255)
                                    tGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGc                                                                                                               <  1:795030/148‑1 (MQ=255)
                                    tGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGc                                                                                                               <  1:795328/148‑1 (MQ=255)
                                     gAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCt                                                                                                              >  2:2593512/1‑148 (MQ=255)
                                     gAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCt                                                                                                              >  2:463288/1‑148 (MQ=255)
                                      aGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTg                                                                                                             >  1:1010679/1‑148 (MQ=255)
                                         gcgcgTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGTGCGATGATACGCCGCAGGCGCTGGATGCTGTCt                                                                                                          <  2:2176386/148‑1 (MQ=255)
                                               cTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaag                                                                                                    <  1:3303893/148‑1 (MQ=255)
                                                  cAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGaagaag                                                                                                 >  1:837090/1‑148 (MQ=255)
                                                  cAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGaagaag                                                                                                 >  1:1033564/1‑148 (MQ=255)
                                                    aCAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAagag                                                                                               >  1:117993/1‑148 (MQ=255)
                                                     cAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGc                                                                                              >  2:2360550/1‑148 (MQ=255)
                                                        cGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAgg                                                                                           >  1:753579/1‑148 (MQ=255)
                                                         gTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGt                                                                                          <  2:3179019/148‑1 (MQ=255)
                                                            aaaGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTaa                                                                                       <  2:1684813/148‑1 (MQ=255)
                                                            aaaGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTaa                                                                                       >  1:2831408/1‑148 (MQ=255)
                                                               gCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAAtt                                                                                    >  1:673975/1‑148 (MQ=255)
                                                                ccGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTc                                                                                   >  1:908759/1‑148 (MQ=255)
                                                                 cGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCt                                                                                  <  1:2856734/148‑1 (MQ=255)
                                                                 cGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCt                                                                                  <  1:1948258/148‑1 (MQ=255)
                                                                    tATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTg                                                                               <  2:1141066/148‑1 (MQ=255)
                                                                       gCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTg                                                                            <  1:2601263/148‑1 (MQ=255)
                                                                            gCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACaaa                                                                       >  1:2614867/1‑148 (MQ=255)
                                                                            gCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACaaa                                                                       >  1:592310/1‑148 (MQ=255)
                                                                             cTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAg                                                                      <  2:3275187/148‑1 (MQ=255)
                                                                                tttGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTaaa                                                                                      >  1:1472129/1‑129 (MQ=255)
                                                                                tttGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCt                                                                   >  1:778411/1‑148 (MQ=255)
                                                                                tttGATCTGCTGCTGAAACTGGCAAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCt                                                                   >  2:1392389/1‑148 (MQ=255)
                                                                                 ttGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGACGAGCAGGTTAAATTCTTTGCTGACAAAGGCTg                                                                  >  2:487640/1‑148 (MQ=255)
                                                                                    aTCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCgg                                                               >  1:2998889/1‑148 (MQ=255)
                                                                                     tCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGt                                                              >  2:1324788/1‑148 (MQ=255)
                                                                                     tCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGt                                                              <  1:3078766/148‑1 (MQ=255)
                                                                                       tgctgctgAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGa                                                            <  1:234791/148‑1 (MQ=255)
                                                                                          tgctgAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTc                                                         <  2:333537/148‑1 (MQ=255)
                                                                                          tgctgAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTc                                                         <  1:3209780/148‑1 (MQ=255)
                                                                                              gAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCa                                                     <  1:3213875/148‑1 (MQ=255)
                                                                                               aaaCTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCaa                                                    >  2:410013/1‑148 (MQ=255)
                                                                                                aaCTGGCGAAGACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAAt                                                   <  2:3033280/148‑1 (MQ=255)
                                                                                                  cTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCa                                                 <  1:2098114/148‑1 (MQ=255)
                                                                                                    ggCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACt                                               >  1:2532074/1‑148 (MQ=255)
                                                                                                     gCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTg                                              >  2:479206/1‑148 (MQ=255)
                                                                                                        aaTACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTgcgc                                           >  1:3100750/1‑148 (MQ=255)
                                                                                                        aaTACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTgcgc                                           >  2:1153051/1‑148 (MQ=255)
                                                                                                        aaTACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTgcgc                                           <  1:1173645/148‑1 (MQ=255)
                                                                                                         aTACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCa                                          >  1:1163887/1‑148 (MQ=255)
                                                                                                          tACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCTACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTAAGCAATCACTGCGCAt                                         <  1:516918/148‑1 (MQ=255)
                                                                                                           aCCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATc                                        >  2:2698556/1‑148 (MQ=255)
                                                                                                               cTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGTTGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCg                                    <  1:348900/148‑1 (MQ=255)
                                                                                                                tATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGt                                   <  1:2805173/148‑1 (MQ=255)
                                                                                                                 aTCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTc                                  >  2:1880441/1‑148 (MQ=255)
                                                                                                                  tCGACAGCGACACAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCg                                 <  1:2281216/148‑1 (MQ=255)
                                                                                                                    gACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCt                               <  2:429237/148‑1 (MQ=255)
                                                                                                                    gACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCt                               <  2:429223/148‑1 (MQ=255)
                                                                                                                        gCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGc                           >  2:2168218/1‑148 (MQ=255)
                                                                                                                         cGACAAAGTGGCGGAAAAAATTGCCTACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGcc                          >  2:911485/1‑148 (MQ=255)
                                                                                                                         cGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGcc                          >  1:729817/1‑148 (MQ=255)
                                                                                                                         cGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGcc                          >  2:668100/1‑148 (MQ=255)
                                                                                                                                gTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCg                   >  2:1416053/1‑148 (MQ=255)
                                                                                                                                gTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCg                   <  1:2593571/148‑1 (MQ=255)
                                                                                                                                 tGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGt                  >  1:1322946/1‑148 (MQ=255)
                                                                                                                                      gAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCgcgc             <  2:1319123/148‑1 (MQ=255)
                                                                                                                                        aaaaaTTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGcgccg           <  2:2940298/148‑1 (MQ=255)
                                                                                                                                          aaaTTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCa         <  1:2898564/148‑1 (MQ=255)
                                                                                                                                             ttGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAata      <  1:1508915/148‑1 (MQ=255)
                                                                                                                                              tGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAatat     >  2:2248092/1‑148 (MQ=255)
                                                                                                                                              tGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAatat     >  2:1114808/1‑148 (MQ=255)
                                                                                                                                                 cgacgaAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTg  >  2:1698939/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGTTAAATTCTTTGCTGACAAAGGCTGCGGTGATTCAGCAATCACTGCGCATCAGCGTCGCTTTGCCGGTGCCGTCGCGCCGCAATATCTG  >  CP009273/1280986‑1281278

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: