Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,282,099 C→T R1005C (CGC→TGC)  narG → nitrate reductase 1 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,282,0990CT86.5% 687.1 / 68.5 251R1005C (CGC→TGC) narGnitrate reductase 1 alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (14/20);  new base T (117/100);  total (131/120)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

Only 100 of 258 total aligned reads displayed.

CTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCC  >  NC_000913/1282009‑1282190
                                                                                          |                                                                                           
cTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCgtgc                                                                                           >  2:6454957/1‑93 (MQ=255)
cTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCgtgc                                                                                           >  2:2010361/1‑93 (MQ=255)
cTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCgtgc                                                                                           >  2:5325764/1‑93 (MQ=255)
 tGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCa                                                                                          >  2:4434598/1‑93 (MQ=255)
 tGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCa                                                                                          >  2:2542910/1‑93 (MQ=255)
  gAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaa                                                                                         >  2:4977196/1‑93 (MQ=255)
  gAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaa                                                                                         >  1:7199915/1‑93 (MQ=255)
    gCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCaaaa                                                                                       >  2:650895/1‑93 (MQ=255)
    gCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCaaaa                                                                                       >  2:1397986/1‑93 (MQ=255)
     cGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAAt                                                                                      <  2:4694968/93‑1 (MQ=255)
     cGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAAt                                                                                      >  2:2398668/1‑93 (MQ=255)
     agAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAAt                                                                                      <  1:2354127/92‑1 (MQ=255)
       aaTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTa                                                                                    >  1:3070738/1‑93 (MQ=255)
        aTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTAt                                                                                   <  1:1814574/93‑1 (MQ=255)
        aTTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTAt                                                                                   >  2:3592670/1‑93 (MQ=255)
         tttACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTAtc                                                                                  >  2:3657520/1‑93 (MQ=255)
         tttACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAtc                                                                                  <  1:3760594/93‑1 (MQ=255)
         tttACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAtc                                                                                  <  1:1369177/93‑1 (MQ=255)
           tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTAtctc                                                                                <  1:5993604/93‑1 (MQ=255)
           tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTAtctc                                                                                >  2:1077733/1‑93 (MQ=255)
           tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTAtctc                                                                                <  2:3249532/93‑1 (MQ=255)
           tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAtctc                                                                                <  1:3281473/93‑1 (MQ=255)
            aCCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTcc                                                                               <  1:1550266/93‑1 (MQ=255)
             ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCa                                                                              >  2:7349134/1‑93 (MQ=255)
             ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCa                                                                              >  1:5907315/1‑93 (MQ=255)
             ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCa                                                                              <  2:194681/93‑1 (MQ=255)
               ggTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGc                                                                            >  2:4027150/1‑93 (MQ=255)
                gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGcc                                                                           >  2:1256411/1‑93 (MQ=255)
                gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGcc                                                                           >  1:2586723/1‑93 (MQ=255)
                gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGcc                                                                           >  2:6748724/1‑93 (MQ=255)
                gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGcc                                                                           >  1:1682986/1‑93 (MQ=255)
                 tCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGccc                                                                          >  1:2663897/1‑93 (MQ=255)
                   gTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGa                                                                        >  2:177675/1‑93 (MQ=255)
                    tGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGAc                                                                       >  1:530485/1‑93 (MQ=255)
                      aCCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACct                                                                     >  2:7445763/1‑93 (MQ=255)
                       ccATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACctg                                                                    <  1:4632043/93‑1 (MQ=255)
                       ccATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACctg                                                                    <  2:2931864/93‑1 (MQ=255)
                        cATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACctgg                                                                   >  1:6922390/1‑93 (MQ=255)
                         aTACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACctggt                                                                  >  1:4783414/1‑93 (MQ=255)
                             gCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACctggtctgg                                                              >  2:7147446/1‑93 (MQ=255)
                              cATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACctggtctggt                                                             >  1:2211071/1‑93 (MQ=255)
                                 cTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGgtctggtctg                                                          <  2:2321550/93‑1 (MQ=255)
                                  tGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGtctggtctgg                                                         <  2:3420925/93‑1 (MQ=255)
                                  tGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGtctggtctgg                                                         >  2:6222379/1‑93 (MQ=255)
                                    gcgcTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGaa                                                       <  1:5078403/93‑1 (MQ=255)
                                       cTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt                                                    >  2:5669925/1‑93 (MQ=255)
                                       cTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt                                                    <  2:6619645/93‑1 (MQ=255)
                                        tGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATg                                                   <  1:3533302/93‑1 (MQ=255)
                                        tGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATg                                                   <  2:6694773/93‑1 (MQ=255)
                                            tAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAaca                                               <  1:1350523/93‑1 (MQ=255)
                                            aaaaGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCCCCAGCCCGACCTGGTCTGGTCTGGAAGATGAaca                                               <  1:2451854/92‑1 (MQ=255)
                                             aaaGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac                                              <  2:2685636/93‑1 (MQ=255)
                                             aaaGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac                                              <  2:2859947/93‑1 (MQ=255)
                                             aaaGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac                                              <  1:3339988/93‑1 (MQ=255)
                                              aagaagACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGGCCTGGTCTGGTCTGGAAGATGAACACg                                             >  1:3041056/1‑93 (MQ=255)
                                               agaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt                                            >  2:4493058/1‑93 (MQ=255)
                                                 aagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGttt                                          <  1:2875067/93‑1 (MQ=255)
                                                  agaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTc                                         <  2:5904277/93‑1 (MQ=255)
                                                  agaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTc                                         <  1:6990217/93‑1 (MQ=255)
                                                  agaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGCCTGGTCTGGAAGATGAACACGTTTc                                         <  1:4406756/93‑1 (MQ=255)
                                                    aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGCCTGGTCTGGAAGATGAACACGTTTCtt                                       <  2:4259819/93‑1 (MQ=255)
                                                     cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa                                      <  2:792052/93‑1 (MQ=255)
                                                     cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa                                      <  1:6737402/93‑1 (MQ=255)
                                                      gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc                                     <  1:3634527/93‑1 (MQ=255)
                                                      gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc                                     >  1:5833638/1‑93 (MQ=255)
                                                       agaagaTCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                    >  2:1808133/1‑93 (MQ=255)
                                                       agaagaTCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                    >  1:4196768/1‑93 (MQ=255)
                                                        gaagaTCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACaa                                   >  2:4175723/1‑93 (MQ=255)
                                                         aagaTCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAAc                                  <  1:6065841/93‑1 (MQ=255)
                                                         aagaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAAc                                  <  1:1004180/93‑1 (MQ=255)
                                                           gaTCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGc                                >  1:1521807/1‑93 (MQ=255)
                                                           gaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGCCTGGTCTGGAAGATGAACACGTTTCTTACAACGc                                >  2:5610047/1‑93 (MQ=255)
                                                            aTCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc                               >  1:636270/1‑93 (MQ=255)
                                                             tCCGCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                              <  1:1947575/93‑1 (MQ=255)
                                                              ccGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGCCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCgg                             <  2:3379823/93‑1 (MQ=255)
                                                                gCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGtt                           >  2:7118098/1‑93 (MQ=255)
                                                                gCTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGtt                           <  2:1473981/93‑1 (MQ=255)
                                                                 cTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                          <  2:6066428/93‑1 (MQ=255)
                                                                 cTTCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTa                          >  1:5037934/1‑93 (MQ=255)
                                                                  ttCCGCGATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTac                         >  2:2198605/1‑93 (MQ=255)
                                                                        gATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc                   >  2:3835475/1‑93 (MQ=255)
                                                                        gATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc                   <  2:1634846/93‑1 (MQ=255)
                                                                        gATATTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAAc                   >  1:6568664/1‑93 (MQ=255)
                                                                          tatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGCCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGt                 <  1:1062051/93‑1 (MQ=255)
                                                                           atTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                <  2:3815072/93‑1 (MQ=255)
                                                                           atTCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                >  2:6888259/1‑93 (MQ=255)
                                                                            ttCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTc               <  1:4303208/93‑1 (MQ=255)
                                                                             tCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCa              >  1:5575535/1‑93 (MQ=255)
                                                                             tCAGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCa              <  2:302790/93‑1 (MQ=255)
                                                                               aGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACg            >  1:4942730/1‑93 (MQ=255)
                                                                               aGGCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACg            <  1:2524304/93‑1 (MQ=255)
                                                                                ggCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa           <  2:5816888/93‑1 (MQ=255)
                                                                                ggCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa           <  1:7519182/93‑1 (MQ=255)
                                                                                ggCACAGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa           >  1:3376540/1‑93 (MQ=255)
                                                                                  cacaGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc         <  1:5809173/93‑1 (MQ=255)
                                                                                  cacaGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc         <  1:7462843/93‑1 (MQ=255)
                                                                                   acaGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCt        >  1:3859646/1‑93 (MQ=255)
                                                                                   acaGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCt        <  1:950690/93‑1 (MQ=255)
                                                                                   acaGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCt        >  2:1749436/1‑93 (MQ=255)
                                                                                    caGCCGTGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTg       >  2:212541/1‑93 (MQ=255)
                                                                                          |                                                                                           
CTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCC  >  NC_000913/1282009‑1282190

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