Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,826,728 A→G T507A (ACC→GCC)  recG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,826,7280AG94.4% 256.9 / 2.6 72T507A (ACC→GCC) recGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (38/30);  total (40/32)
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.85e-01

TGACTACGGTCGCTATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTC  >  NC_000913/3826582‑3826874
                                                                                                                                                  |                                                                                                                                                   
tGACTACGGTCGCTATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTg                                                                                                                                                 >  2:611146/1‑151 (MQ=255)
 gACTACGGTCGCTATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTgg                                                                                                                                                >  1:1063051/1‑151 (MQ=255)
      gggtctttATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAag                                                                                                                                            <  1:593713/143‑1 (MQ=255)
          cctttttCCTGTTCCCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTg                                                                                                                                       <  2:1141011/146‑1 (MQ=255)
                 cctaaaccccGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGaaa                                                                                                                                    <  1:377478/140‑1 (MQ=255)
                 ccggATCCCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACt                                                                                                                                  <  2:835702/146‑1 (MQ=255)
                 ccTAATCCCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCt                                                                                                                                                  <  1:544894/133‑1 (MQ=255)
                     agaCCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGaaa                                                                                                                                    <  2:1066747/141‑1 (MQ=255)
                           gccgTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTAcc                                                                                                                         >  2:1010392/1‑148 (MQ=255)
                           gccgTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCaga                                                                                                                      <  1:508137/151‑1 (MQ=255)
                             cgTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGt                                                                                                                    <  1:314670/151‑1 (MQ=255)
                                aCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTg                                                                                                                  >  2:1239751/1‑150 (MQ=255)
                                   gACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACg                                                                                                              >  1:305529/1‑151 (MQ=255)
                                      atcatTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTg                                                                                                           <  2:325977/151‑1 (MQ=255)
                                           tGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCt                                                                                                       <  1:422797/150‑1 (MQ=255)
                                            gACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTgg                                                                                                     <  2:53164/151‑1 (MQ=255)
                                            gACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTg                                                                                                      <  1:586592/150‑1 (MQ=255)
                                             aCCGCGTGCACCACGCCTGCATAACTGAAGGTCGACAGGCATTCTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGTAGCTGCCTGGGAAGAGTTGAAACTGGCGCTAACAGAGTTGAACGTTAGCCTgg                                                                                                     >  2:781664/1‑150 (MQ=255)
                                               cgcgTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTac                                                                                                  >  1:108675/1‑151 (MQ=255)
                                                gcgTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTac                                                                                                  >  2:360605/1‑150 (MQ=255)
                                                  gTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACg                                                                                               >  2:1218779/1‑151 (MQ=255)
                                                   tGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACgg                                                                                              >  1:426787/1‑151 (MQ=255)
                                                    gCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACgg                                                                                              <  2:741817/150‑1 (MQ=255)
                                                     caccacGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCATGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGc                                                                                            >  2:1000182/1‑151 (MQ=255)
                                                      accacGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCg                                                                                           >  1:784603/1‑151 (MQ=255)
                                                         acGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGAt                                                                                        >  2:257942/1‑151 (MQ=255)
                                                         acGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGAt                                                                                        >  2:1136157/1‑151 (MQ=255)
                                                            ccTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGa                                                                                      >  2:951465/1‑150 (MQ=255)
                                                              tGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAAc                                                                                   >  2:77798/1‑151 (MQ=255)
                                                               gCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAAcc                                                                                  <  2:797112/151‑1 (MQ=255)
                                                                   aaCTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGcc                                                                              >  2:1255813/1‑151 (MQ=255)
                                                                       gAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCga                                                                            <  2:189193/149‑1 (MQ=255)
                                                                           gtcgtcAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACa                                                                      <  2:801055/151‑1 (MQ=255)
                                                                             cgtcAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAgg                                                                    >  2:114542/1‑151 (MQ=255)
                                                                              gtcAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGc                                                                   <  1:143127/151‑1 (MQ=255)
                                                                                cAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCg                                                                  <  1:15221/150‑1 (MQ=255)
                                                                                   gCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTg                                                               >  2:854503/1‑150 (MQ=255)
                                                                                    cATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTg                                                               >  1:1235799/1‑149 (MQ=255)
                                                                                      tACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATg                                                            >  1:710326/1‑150 (MQ=255)
                                                                                         tGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATgg                                                           >  1:1130337/1‑148 (MQ=255)
                                                                                          gggTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGTCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGcgtc                                                      <  1:1062932/151‑1 (MQ=255)
                                                                                              ttGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGtt                                                   <  2:7732/150‑1 (MQ=255)
                                                                                               tttACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGttt                                                  >  2:232008/3‑150 (MQ=255)
                                                                                                gTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTa                                                 <  1:504880/150‑1 (MQ=255)
                                                                                                 tACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTaa                                                >  2:556912/1‑150 (MQ=255)
                                                                                                      tgattgaAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGGAGGCGGCGGTTTAACAAgg                                          >  1:981763/1‑151 (MQ=255)
                                                                                                       gattgaAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAgg                                          <  2:1036019/150‑1 (MQ=255)
                                                                                                        attgaAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTg                                        <  2:547068/151‑1 (MQ=255)
                                                                                                         ttgaAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTg                                        >  2:1047494/1‑150 (MQ=255)
                                                                                                         ttgaAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTg                                        >  2:737112/1‑150 (MQ=255)
                                                                                                             agagTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCt                                    >  2:288134/1‑150 (MQ=255)
                                                                                                                gTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTacac                                <  2:117631/151‑1 (MQ=255)
                                                                                                                gTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTaca                                 <  1:863792/150‑1 (MQ=255)
                                                                                                                  cGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACct                              <  2:612910/151‑1 (MQ=255)
                                                                                                                  cGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACct                              <  2:690260/151‑1 (MQ=255)
                                                                                                                   ggAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACctg                             >  1:200306/1‑151 (MQ=255)
                                                                                                                   ggAATTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAAAAAGGTGAGCTaca                                 >  1:562332/1‑147 (MQ=255)
                                                                                                                      aTTGCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACctgctg                          <  2:441052/151‑1 (MQ=255)
                                                                                                                         gCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTTAACAAGGTGAGCTACACCTGCTGGtt                       >  1:377645/1‑151 (MQ=255)
                                                                                                                         gCTGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGAAGGCGTCGGTTAAACAAGGTGAGCTACACCCGCTgg                         >  1:807747/1‑149 (MQ=255)
                                                                                                                           tGGAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGc                     >  1:617955/1‑151 (MQ=255)
                                                                                                                            ggAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGGGGCGTTTAAACAAGGTGAGCTAAACCTGCTGGTTGcc                    >  2:815816/1‑151 (MQ=255)
                                                                                                                            ggAAGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGc                     <  2:11003/150‑1 (MQ=255)
                                                                                                                               aGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCcaca                 <  2:1194245/151‑1 (MQ=255)
                                                                                                                               aGCGCAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAAAAAGGTGAGCTACACCTGCTGGTTGCcac                  >  1:848097/1‑150 (MQ=255)
                                                                                                                                   cAGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCg             <  1:276066/151‑1 (MQ=255)
                                                                                                                                    aGGCGGCGGAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGGTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGGTGGCGACCATTTAAACCAGGGAGCTCAACCTACTGGtgtccccaccc              >  1:841547/1‑139 (MQ=255)
                                                                                                                                        ggcggAAGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGGGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTAtt        >  2:440175/1‑151 (MQ=255)
                                                                                                                                             aaGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGaa     >  1:425321/1‑149 (MQ=255)
                                                                                                                                              aGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTc  >  2:860987/1‑151 (MQ=255)
                                                                                                                                              aGCTGCCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATAAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTc  <  1:674288/151‑1 (MQ=255)
                                                                                                                                               gCTACCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGGTGGCGTCGTTTAAACAAGGCGAGCTACACCTGCTGGTTGCCACAACCGTTAATGAAGt   >  2:1138304/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGACTACGGTCGCTATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGGCGCTACCAGAGTTGAACGTTGGCCTGGTACACGGGCGGATGAAACCTGCCGAGAAACAGGCGGTGATGGCGTCGTTTAAACAAGGTGAGCTACACCTGCTGGTTGCCACAACCGTTATTGAAGTC  >  NC_000913/3826582‑3826874

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